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

ATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTG  >  NC_000913/1351047‑1351333
                                                                                                                                               |                                                                                                                                                 
aTCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑gggggggg                                                                                                                                           <  2:515875/151‑2 (MQ=255)
  ccTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑ggggggggc                                                                                                                                          <  2:771455/150‑3 (MQ=255)
    ttATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑ggggggggcgct                                                                                                                                       >  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)
                 ttgttgTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑ccc                                                                                                                          >  1:590235/1‑151 (MQ=255)
                 ttgttgTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTACCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCCCA‑ccc                                                                                                                          >  1:669159/1‑151 (MQ=255)
                            gAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCctct                                                                                                                 <  1:670575/149‑1 (MQ=255)
                            gAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTcc                                                                                                               <  2:455707/151‑1 (MQ=255)
                               aaTCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCa                                                                                                            <  2:192918/151‑1 (MQ=255)
                                    ggCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAggg                                                                                                       <  2:440232/151‑1 (MQ=255)
                                      cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGc                                                                                                     <  2:707807/151‑1 (MQ=255)
                                      cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGc                                                                                                     <  2:382960/151‑1 (MQ=255)
                                      cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGc                                                                                                     <  2:192341/151‑1 (MQ=255)
                                      cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCCAACCCTCTCCCCAGAgggg                                                                                                      >  1:429281/1‑150 (MQ=255)
                                        aTTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCg                                                                                                    >  2:118243/1‑150 (MQ=255)
                                         ttCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGC‑TTTCCCCTCA‑CCCTAAACCTCTCCCCCGAGGGGCGa                                                                                                   >  2:582319/1‑149 (MQ=255)
                                          tcgtcgTTTTAGTAAACAGTACGAACAAATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGGTGCCAGGGCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGa                                                                                                   <  1:752451/149‑1 (MQ=255)
                                            gtcgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAgggg                                                                                               >  2:142370/1‑151 (MQ=255)
                                                 tttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCCCC‑CCCCAACCCTCTCCCCCGCGGGGCGAGGGGACCg                                                                                           >  1:109870/1‑150 (MQ=255)
                                                   tAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTa                                                                                         <  2:586056/150‑1 (MQ=255)
                                                       aaaCAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCCCA‑CCCTAACCCTCTCCCCCGAGGGGCGAGGGGACCGTATtgt                                                                                     >  1:504818/1‑150 (MQ=255)
                                                              aCGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTTTGCAAatat                                                                             >  1:425176/1‑151 (MQ=255)
                                                               cGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATAtt                                                                            <  1:186863/151‑1 (MQ=255)
                                                                 aaCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGt                                                                          <  1:708631/151‑1 (MQ=255)
                                                                   cAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTa                                                                        <  2:655865/151‑1 (MQ=255)
                                                                   cAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGATTTTCCCCTTA‑CCCCAACCCTCTCCCCCGAGGGGCGAGGGGACCGGATTGTGGAAATATTGt                                                                          >  1:218675/1‑149 (MQ=255)
                                                                        aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCCCA‑CCCCAACCCCCCCCCCCAAGGGGCGGGGGGGCCGGATTTTGCAAAAATTTTTTcccc                                                                    >  2:733500/1‑150 (MQ=255)
                                                                               tattatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCCCATCCCACCCCCCA‑AGGGGGGGGGGGGCCCACTTCTGCCAAAATTGGTAcccacacaa                                                               >  2:211326/1‑142 (MQ=255)
                                                                               tattatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGGCCGGATTTTGCAAATATTGTTACCCCAGcaa                                                               >  2:142053/1‑148 (MQ=255)
                                                                               tattatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGcaaca                                                             <  1:733426/150‑1 (MQ=255)
                                                                                   ataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAg                                                        <  2:104278/151‑1 (MQ=255)
                                                                                     aaTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGc                                                      <  2:561954/150‑1 (MQ=255)
                                                                                        cAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCCCC‑CCCTAACCCCCTCCCCAGAGGGGGGAGGGGAACCGATTTTGCAAAAATTTTTAACCCCGCCACAAAAAAGCTc                                                    >  2:691151/1‑150 (MQ=255)
                                                                                              ggCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGGATTGTGCAAATATTTTTACCCCCACAACAAACCAGCTCAAACa                                               >  1:657834/1‑149 (MQ=37)
                                                                                                              aGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCAt                              <  2:525780/150‑1 (MQ=255)
                                                                                                              aGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCAt                              <  1:24156/150‑1 (MQ=255)
                                                                                                               gTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCC‑CCCCAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATgg                            >  1:179435/1‑151 (MQ=255)
                                                                                                                    cAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCC‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCCACAAAAAGGGTCCTACAGCCCCTAACCCTTTCATGGCGAt                        >  2:247310/1‑150 (MQ=38)
                                                                                                                              cGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCCGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCCTGGCGATGGCTGGGACgg             >  1:289268/1‑151 (MQ=255)
                                                                                                                                         gCGGCAGGGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTg  <  2:212264/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
ATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑GGGGGGGCGCTTTTCCCCTCA‑CCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTTCAGACCTTG  >  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: ‑