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.G97.6% 137.0 / ‑2.6 42intergenic (‑150/+218)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (22/19);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

ATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTG  >  NC_000913/1351047‑1351333
                                                                                                                                               |                                                                                                                                                 
aTCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAgggggggg                                                                                                                                           <  2:515875/151‑2 (MQ=255)
  ccTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAggggggggc                                                                                                                                          <  2:771455/150‑3 (MQ=255)
    ttATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAggggggggcgct                                                                                                                                       >  1:496132/1‑146 (MQ=255)
     tATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCtt                                                                                                                                       >  1:561882/1‑151 (MQ=255)
       ttgttgATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCtttt                                                                                                                                     >  1:707733/1‑151 (MQ=255)
       ttgttgATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCtttt                                                                                                                                     >  1:192308/1‑151 (MQ=255)
         gttgATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTc                                                                                                                                    <  2:61621/150‑1 (MQ=255)
           tgATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTcccc                                                                                                                                 >  1:369768/1‑151 (MQ=255)
               gCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCa                                                                                                                              <  1:522027/150‑1 (MQ=255)
                 ttgttgTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCAccc                                                                                                                          >  1:590235/1‑151 (MQ=255)
                 ttgttgTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTACCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCCCAccc                                                                                                                          >  1:669159/1‑151 (MQ=255)
                            gAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCctct                                                                                                                 <  1:670575/149‑1 (MQ=255)
                            gAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTcc                                                                                                               <  2:455707/151‑1 (MQ=255)
                               aaTCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCa                                                                                                            <  2:192918/151‑1 (MQ=255)
                                    ggCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAggg                                                                                                       <  2:440232/151‑1 (MQ=255)
                                      cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGc                                                                                                     <  2:707807/151‑1 (MQ=255)
                                      cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGc                                                                                                     <  2:382960/151‑1 (MQ=255)
                                      cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGc                                                                                                     <  2:192341/151‑1 (MQ=255)
                                      cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCCAACCCTCTCCCCAGAgggg                                                                                                      >  1:429281/1‑150 (MQ=255)
                                        aTTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCg                                                                                                    >  2:118243/1‑150 (MQ=255)
                                         ttCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTCCCCTCACCCTAAACCTCTCCCCCGAGGGGCGa                                                                                                   >  2:582319/1‑149 (MQ=255)
                                          tcgtcgTTTTAGTAAACAGTACGAACAAATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGGTGCCAGGGCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGa                                                                                                   <  1:752451/149‑1 (MQ=255)
                                            gtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgggg                                                                                               >  2:142370/1‑151 (MQ=255)
                                                 tttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCCCCCCCCAACCCTCTCCCCCGCGGGGCGAGGGGACCg                                                                                           >  1:109870/1‑150 (MQ=255)
                                                   tAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTa                                                                                         <  2:586056/150‑1 (MQ=255)
                                                       aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCCCACCCTAACCCTCTCCCCCGAGGGGCGAGGGGACCGTATtgt                                                                                     >  1:504818/1‑150 (MQ=255)
                                                              aCGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTTTGCAAatat                                                                             >  1:425176/1‑151 (MQ=255)
                                                               cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATAtt                                                                            <  1:186863/151‑1 (MQ=255)
                                                                 aaCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGt                                                                          <  1:708631/151‑1 (MQ=255)
                                                                   cAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTa                                                                        <  2:655865/151‑1 (MQ=255)
                                                                   cAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGATTTTCCCCTTACCCCAACCCTCTCCCCCGAGGGGCGAGGGGACCGGATTGTGGAAATATTGt                                                                          >  1:218675/1‑149 (MQ=255)
                                                                        aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCCCACCCCAACCCCCCCCCCCAAGGGGCGGGGGGGCCGGATTTTGCAAAAATTTTTTcccc                                                                    >  2:733500/1‑150 (MQ=255)
                                                                               tattatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCCCATCCCACCCCCCAAGGGGGGGGGGGGCCCACTTCTGCCAAAATTGGTAcccacacaa                                                               >  2:211326/1‑142 (MQ=255)
                                                                               tattatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGGCCGGATTTTGCAAATATTGTTACCCCAGcaa                                                               >  2:142053/1‑148 (MQ=255)
                                                                               tattatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGcaaca                                                             <  1:733426/150‑1 (MQ=255)
                                                                                   ataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAg                                                        <  2:104278/151‑1 (MQ=255)
                                                                                     aaTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGc                                                      <  2:561954/150‑1 (MQ=255)
                                                                                        cAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCCCCCCCTAACCCCCTCCCCAGAGGGGGGAGGGGAACCGATTTTGCAAAAATTTTTAACCCCGCCACAAAAAAGCTc                                                    >  2:691151/1‑150 (MQ=255)
                                                                                              ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGGATTGTGCAAATATTTTTACCCCCACAACAAACCAGCTCAAACa                                               >  1:657834/1‑149 (MQ=37)
                                                                                                              aGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCAt                              <  2:525780/150‑1 (MQ=255)
                                                                                                              aGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCAt                              <  1:24156/150‑1 (MQ=255)
                                                                                                               gTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCCCCCCAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATgg                            >  1:179435/1‑151 (MQ=255)
                                                                                                                    cAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCCCCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCCACAAAAAGGGTCCTACAGCCCCTAACCCTTTCATGGCGAt                        >  2:247310/1‑150 (MQ=38)
                                                                                                                              cGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCCGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCCTGGCGATGGCTGGGACgg             >  1:289268/1‑151 (MQ=255)
                                                                                                                                         gCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTg  <  2:212264/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
ATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTG  >  NC_000913/1351047‑1351333

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