Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,800,475 G→A *110* (TAG→TAA)  ygaM → putative membrane‑anchored DUF883 family ribosome‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,800,4750GA96.3% 159.9 / ‑1.6 54*110* (TAG→TAA) ygaMputative membrane‑anchored DUF883 family ribosome‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (26/26);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

CCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGC  >  NC_000913/2800346‑2800606
                                                                                                                                 |                                                                                                                                    
ccAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAAt                                                                                                                             <  1:420791/139‑1 (MQ=255)
                    ccGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGCTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATACTAAATACGC‑GTCttt                                                                                                        <  1:1198500/139‑1 (MQ=255)
                                                                aGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATc                                                            >  2:1584019/1‑139 (MQ=255)
                                                                 gCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTgg                                                                  >  2:129242/1‑132 (MQ=255)
                                                                 gCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTgg                                                                  <  1:129242/132‑1 (MQ=255)
                                                                  cGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                          >  1:1625843/1‑139 (MQ=255)
                                                                  cGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                          >  2:1283165/1‑139 (MQ=255)
                                                                      ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGata                                                                                                                         >  1:1012303/1‑73 (MQ=255)
                                                                      ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGata                                                                                                                         <  2:1012303/73‑1 (MQ=255)
                                                                          cagcagCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcc                                                  <  2:66774/139‑1 (MQ=255)
                                                                           agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAg                                                                                              >  1:1225046/1‑94 (MQ=255)
                                                                           agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAg                                                                                              <  2:1225046/94‑1 (MQ=255)
                                                                           agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATccc                                                 >  1:877215/1‑139 (MQ=255)
                                                                                 tGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAcc                                           >  2:1193053/1‑139 (MQ=255)
                                                                                  gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGGaaa                                                                                                                                <  2:1221252/54‑1 (MQ=39)
                                                                                  gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCaaa                                                                                                                                >  1:1221252/1‑54 (MQ=255)
                                                                                       aggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                     <  1:1554934/139‑1 (MQ=255)
                                                                                         gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                   <  2:901261/139‑1 (MQ=255)
                                                                                             ttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCa                               >  1:1750458/1‑139 (MQ=255)
                                                                                                  tCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCaaaa                                                                                                                               >  2:725597/1‑39 (MQ=255)
                                                                                                  tCGGCGCACTGTTAAGCATGCGCAAATCGTAATGCaaaa                                                                                                                               <  1:725597/39‑1 (MQ=255)
                                                                                                    ggCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtt                                                                           <  1:212260/88‑1 (MQ=255)
                                                                                                    ggCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtt                                                                           >  2:212260/1‑88 (MQ=255)
                                                                                                    ggCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCt                                 <  1:106246/130‑1 (MQ=255)
                                                                                                    ggCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCt                                 >  2:106246/1‑130 (MQ=255)
                                                                                                    ggCGCACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                        >  2:653865/1‑139 (MQ=255)
                                                                                                      cgcACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                     <  1:463606/124‑1 (MQ=255)
                                                                                                      cgcACTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                     >  2:463606/1‑124 (MQ=255)
                                                                                                         aCTGTTAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTaa                   >  1:1152900/1‑139 (MQ=255)
                                                                                                             ttAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcc                                                  <  1:72712/104‑1 (MQ=255)
                                                                                                             ttAAGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcc                                                  >  2:72712/1‑104 (MQ=255)
                                                                                                                aGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCtt                                                                                                         >  1:13312/1‑46 (MQ=255)
                                                                                                                aGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCtt                                                                                                         <  2:13312/46‑1 (MQ=255)
                                                                                                                aGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc            <  1:1193053/139‑1 (MQ=255)
                                                                                                                aGCATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc            <  2:1625843/139‑1 (MQ=255)
                                                                                                                  cATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                          <  2:159403/91‑1 (MQ=255)
                                                                                                                  cATGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcac                                                          >  1:159403/1‑91 (MQ=255)
                                                                                                                    tggccAAATCCTTATGGAAAAATGATTATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACt        <  1:1523839/135‑1 (MQ=255)
                                                                                                                    tGCGCAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACt        >  1:90122/1‑139 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATAcgc                                                                                                               <  2:302283/36‑1 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATAcgc                                                                                                               >  1:302283/1‑36 (MQ=38)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTg                                                                       >  1:1372727/1‑75 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTg                                                                       <  2:1372727/75‑1 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                   >  2:425297/1‑111 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                   <  1:425297/111‑1 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTa       >  2:1481652/1‑139 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTa       >  1:683882/1‑139 (MQ=255)
                                                                                                                     gcgcAAATCGTAATGCAAAAATGATAAT‑AATAAGCGGCCTTTAACCCGAAAGCCTGTCTTGGGGGTTTTTTTTTCCTTGTGGAATAAAAAAATTCCCCCTCCCCCTTCACCCTCTTATCAAGGGTATTTTGGACAAcat       >  2:1523839/1‑137 (MQ=255)
                                                                                                                      cgcAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt      <  1:1584019/139‑1 (MQ=255)
                                                                                                                      cgcAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt      <  2:683882/139‑1 (MQ=255)
                                                                                                                      cgcAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt      <  2:1152900/139‑1 (MQ=255)
                                                                                                                      cgcAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGAGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt      <  1:1283165/139‑1 (MQ=255)
                                                                                                                       gcAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGa                                                                                                >  2:289552/1‑48 (MQ=255)
                                                                                                                       gcAAATCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGa                                                                                                <  1:289552/48‑1 (MQ=255)
                                                                                                                          aaTCGTAATGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGc  <  2:1750458/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                    
CCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGC‑GTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGC  >  NC_000913/2800346‑2800606

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑