Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,351,196 (G)7→8 intergenic (‑157/+211) fabI ← / ← ycjD enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,351,1891.G94.2% 176.9 / 3.0 52intergenic (‑150/+218)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (24/25);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

GATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGA  >  NC_000913/1351059‑1351319
                                                                                                                                   |                                                                                                                                   
gATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑gggggggg                                                                                                                            >  1:949852/1‑138 (MQ=255)
 aTGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑GGGGGGgcg                                                                                                                           >  1:848929/1‑139 (MQ=255)
 aTGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑GGGGGGgcg                                                                                                                           >  2:784773/1‑139 (MQ=255)
   gCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATTAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑ggggggggcgc                                                                                                                         >  1:1680028/1‑135 (MQ=255)
    cTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑GGGGGGGCGCtt                                                                                                                        <  1:784773/139‑1 (MQ=255)
        ttgtGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑gggggggg                                                                                                                            >  2:441731/1‑130 (MQ=255)
        ttgtGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑gggggggg                                                                                                                            <  1:441731/131‑2 (MQ=255)
                          cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCctct                                                                                                   <  1:586933/139‑1 (MQ=255)
                           gATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCct                                                                                                     <  2:244251/136‑1 (MQ=255)
                           gATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCct                                                                                                     >  1:244251/1‑136 (MQ=255)
                                       tAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCg                                                                                      >  2:1440330/1‑139 (MQ=255)
                                              cAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAAccc                                                                                                      >  1:346840/1‑116 (MQ=255)
                                              cAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAAccc                                                                                                      <  2:346840/116‑1 (MQ=255)
                                                 tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCt                                                                                                                 >  2:716631/1‑102 (MQ=255)
                                                 tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCt                                                                                                                 >  1:1672930/1‑102 (MQ=255)
                                                 tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCt                                                                                                                 <  1:716631/102‑1 (MQ=255)
                                                 tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCt                                                                                                                 <  2:1672930/102‑1 (MQ=255)
                                                 tACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGAGAGGGGACCGt                                                                            >  2:269457/1‑139 (MQ=255)
                                                  aCGAACAGATAAACGGTTATTATAATTAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTa                                                                           <  2:1680028/139‑1 (MQ=255)
                                                  aCGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTa                                                                           <  2:1281223/139‑1 (MQ=255)
                                                   cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCct                                                                                                     <  2:736846/112‑1 (MQ=255)
                                                   cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCct                                                                                                     >  1:736846/1‑112 (MQ=255)
                                                   cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAt                                                                          >  1:816539/1‑139 (MQ=255)
                                                                  ttattaTAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGc                                                                                       >  2:456262/1‑111 (MQ=255)
                                                                  ttattaTAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGc                                                                                       <  1:456262/111‑1 (MQ=255)
                                                                    attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTcccc                                                                                                                  >  2:1007428/1‑82 (MQ=255)
                                                                    attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTcccc                                                                                                                  <  1:1007428/82‑1 (MQ=255)
                                                                    attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                                                                     >  1:1697731/1‑111 (MQ=255)
                                                                    attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                                                                     <  2:1697731/111‑1 (MQ=255)
                                                                    attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAc                                                         >  1:793705/1‑139 (MQ=255)
                                                                           tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCt                                                                                                                 <  1:1640380/76‑1 (MQ=255)
                                                                           tCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCt                                                                                                                 >  2:1640380/1‑76 (MQ=255)
                                                                            cAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGcaa                                                 <  2:793705/139‑1 (MQ=255)
                                                                               ccTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTa                                                                           >  1:662953/1‑110 (MQ=255)
                                                                               ccTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTa                                                                           <  2:662953/110‑1 (MQ=255)
                                                                               ccTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCCCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGcaacaa                                              >  1:245838/1‑139 (MQ=255)
                                                                                cTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaa                                             <  2:797507/139‑1 (MQ=255)
                                                                                  ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACGCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAccc                                                       >  2:1083909/1‑127 (MQ=255)
                                                                                  ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAccc                                                       <  1:1083909/127‑1 (MQ=255)
                                                                                   gCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCt                                                                                                                 <  1:1232058/68‑1 (MQ=255)
                                                                                   gCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCt                                                                                                                 >  2:1232058/1‑68 (MQ=255)
                                                                                   gCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCt                                                                                                                 >  2:355425/1‑68 (MQ=255)
                                                                                   gCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCt                                                                                                                 <  1:355425/68‑1 (MQ=255)
                                                                                    cTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAgg                                         <  1:535415/139‑1 (MQ=255)
                                                                                                   gTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaa                                             >  2:311946/1‑120 (MQ=255)
                                                                                                   gTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaa                                             <  1:311946/120‑1 (MQ=255)
                                                                                                        cAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAAcc                                                                                                       <  2:610097/57‑1 (MQ=38)
                                                                                                        cAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAAcc                                                                                                       >  1:610097/1‑57 (MQ=38)
                                                                                                                aCCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGc             >  2:427133/1‑139 (MQ=255)
                                                                                                                     aGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCaga                                                                                            <  2:667112/55‑1 (MQ=21)
                                                                                                                     aGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCaga                                                                                            >  1:667112/1‑55 (MQ=21)
                                                                                                                      gCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGc       <  1:1444625/139‑1 (MQ=255)
                                                                                                                          gCTGCGGCAGGGGGGGGGCGCTTTTCCCCTCACCCTAAcc                                                                                                       >  1:1244139/1‑40 (MQ=18)
                                                                                                                          gCTGCGGCAGGGGGGGGGCGCTTTTCCCCTCACCCTAAcc                                                                                                       <  2:1244139/40‑1 (MQ=18)
                                                                                                                          gCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTggg   >  2:871422/1‑139 (MQ=255)
                                                                                                                           cTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGa  <  2:949852/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                   
GATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGA  >  NC_000913/1351059‑1351319

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