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.G91.7% 145.9 / 6.7 48intergenic (‑150/+218)fabI/ycjDenoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (26/18);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTT  >  NC_000913/1351041‑1351324
                                                                                                                                                     |                                                                                                                                           
gCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑gg                                                                                                                                          >  1:365079/1‑151 (MQ=255)
  tttAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑gggg                                                                                                                                       >  2:309604/1‑151 (MQ=255)
           tATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTAGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAACAGGCTGCGGCA‑‑ggggggggcgc                                                                                                                                >  1:309430/1‑145 (MQ=255)
            aTTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCttt                                                                                                                              <  1:784476/151‑1 (MQ=255)
            aTTGTTGAGGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCttt                                                                                                                              >  1:542191/1‑151 (MQ=255)
                  gATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCt                                                                                                                        >  2:250706/1‑151 (MQ=255)
                           tgtgCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCAccc                                                                                                                   >  1:464279/1‑147 (MQ=255)
                              gCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAAccc                                                                                                             <  1:580923/150‑1 (MQ=255)
                                 tGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCCCCCTAACCCtctc                                                                                                         >  1:461829/1‑151 (MQ=255)
                                       tCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCCCACCCTAACCCCCTaccc                                                                                                      >  2:507547/1‑148 (MQ=255)
                                           gCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAgggg                                                                                               <  1:741427/151‑1 (MQ=255)
                                            cGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCAAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTTAACCCCTTCCCCGAgggg                                                                                               >  2:731834/1‑150 (MQ=255)
                                               ttCGTCGTTTTAGTAAACAGTACGAACAGATCAACGGTTATTATATTCAACCTGGCTGTGAGTATTTATAGTTACTAGGTCCGCCCGGAGCAGGCTGCGGCC‑‑GGGAGGGGGGTTTTTCTCCTCCCCCTAAATCTCTCCCCaaag                                                                                                  >  1:558695/1‑144 (MQ=255)
                                                    cgtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAgggg                                                                                        >  1:614031/1‑149 (MQ=255)
                                                     gtTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCCCCCCCTCACCCCCCCCCCAGAGGGGGGAgggg                                                                                        >  2:696914/1‑148 (MQ=255)
                                                        ttAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTTCCCTCACCCCAACCCTCTCCCCCGAGGGGTGAGGGGAcc                                                                                     >  2:729709/1‑148 (MQ=255)
                                                          aGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCCGAGGGGCGAGGGGACCGTAtt                                                                                >  1:26553/1‑151 (MQ=255)
                                                           gTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTAtt                                                                                >  1:339016/1‑150 (MQ=255)
                                                           gTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATtgt                                                                              <  2:328808/151‑1 (MQ=255)
                                                                  gTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAat                                                                       <  1:588635/151‑1 (MQ=255)
                                                                    aCGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAata                                                                      <  1:490444/150‑1 (MQ=255)
                                                                       aaCAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGtt                                                                 <  2:134407/151‑1 (MQ=255)
                                                                          aGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGGCCGGATTGTGCAAATATTGtt                                                                 >  1:47221/1‑149 (MQ=255)
                                                                           gATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTAcc                                                              <  2:461865/151‑1 (MQ=255)
                                                                              aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCa                                                           <  2:47243/151‑1 (MQ=255)
                                                                              aaaCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCCAGGGGGCCGCTTTGTGCAAATATTTTTAcccc                                                            >  1:442258/1‑150 (MQ=38)
                                                                                aCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTTTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAg                                                          >  2:248500/1‑150 (MQ=255)
                                                                                aCGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCAACCTTACCCTCCCCCCAAGGGGGGGAGGGGGACGTATTGTGCAAAAAATTTTTAccc                                                             >  1:855557/1‑148 (MQ=255)
                                                                                     tattatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCCCTCCCCAGAGGGGGGAGGGGGCCGGATTTTGCAAATATTTTTTcccccacaaaa                                                     >  2:454861/1‑143 (MQ=255)
                                                                                      attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGcaacaa                                                    <  2:154722/150‑1 (MQ=255)
                                                                                      attatAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGAAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCCAACCCCCCCCCCAGAAGGGCGGGGGGGCCGTCTCGGGCAAATCTTTTTACCCCAGcaa                                                       >  1:781288/1‑147 (MQ=255)
                                                                                         ataATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCaacaaaca                                                 <  2:594644/150‑1 (MQ=255)
                                                                                           aaTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAgg                                               <  2:84058/150‑1 (MQ=255)
                                                                                            aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTTCCCCTAACCCTCTCCCCCGAGGGGGGGGGGGCACGTATTTTGCAAATATTTTTAACCCCGCaacaaaaa                                                 >  1:433361/1‑145 (MQ=255)
                                                                                            aTCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCCAACCCTCCCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAgg                                               >  1:396290/1‑149 (MQ=255)
                                                                                               aaCCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCAt                                          <  2:723457/151‑1 (MQ=255)
                                                                                                   tGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCCCCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCCAATATTGTTACCCCAGCAAAAAACAGGCTCATACa                                      >  2:105654/1‑150 (MQ=255)
                                                                                                      atgtgAGTAGCTTTAGTTGCGAGGTCCGACCGGAGCAGGCTGCATCAGGGGGGGGCGCTTTTCCCCTCCTCCTAACCCTTTCCCCACTGGGGGGGGGGGACCCTATTGTGCAAATATTTTTACCCCAGCAACAAAAAGGCTTATACAGcc                                   >  1:36056/2‑150 (MQ=255)
                                                                                                        gtgAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCt                                <  1:154666/151‑1 (MQ=255)
                                                                                                          gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCCGAGGGGCGGGGGGACCGTATTGTGCCAATATTTTTACCCCCGCAACAAACAGGCTCATACAGcccc                                 >  2:686571/1‑148 (MQ=38)
                                                                                                          gAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAAACCTCTCCCCCGAGGGGCGGGGGGGCCGGATTGTGCAAATATTGTTAACCCAGCAACAAACAGGCTCATACAGccccca                               >  2:837666/1‑148 (MQ=255)
                                                                                                           aGTAGCTATAGTAGCCAGGTCTGACCGGAGCAGGCTGTGGCA‑‑GGGGGGGCTCTTTTCTCCTCACCCTAACCCTCTCCCCAGACGGGCGAGGGGACCGTATTGTGTAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCGCCTaa                              <  1:344357/149‑1 (MQ=38)
                                                                                                             tAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGGGGGGACCGTATTGTGCCAATATTGTTACCCCCGCAAAAAACAGGCTCATACAACCCCTAAccc                           >  2:90300/1‑151 (MQ=38)
                                                                                                               gCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTTAACCTCTCCCCCGAGGGGCGGGGGGACCGGATTGTGCAAATATTTTTACCCCCGCAACAAACAAGGTCATAACGGCCCCAAAcctt                          >  2:148244/1‑149 (MQ=255)
                                                                                                                  ataGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCtt                         <  2:485146/149‑1 (MQ=255)
                                                                                                                    aGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATg                    <  2:570053/151‑1 (MQ=255)
                                                                                                                         ccAGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCCGAGGGGCGAGGGGGCCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGa                >  1:388322/1‑150 (MQ=255)
                                                                                                                           aGGTCCGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATgg             <  1:300063/151‑1 (MQ=255)
                                                                                                                                cGACCGGAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTggg        <  1:681319/151‑1 (MQ=255)
                                                                                                                                     ggAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGt   <  1:674120/151‑1 (MQ=255)
                                                                                                                                      gAGCAGGCTGCGGCAGGGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGtt  <  2:725459/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GCTTTAATCCTTATTGTTGATGCTTGTTGTGCCTGAAAATCAGGCGATTCGTCGTTTTAGTAAACAGTACGAACAGATAAACGGTTATTATAATCAACCTGGCTGTGAGTAGCTATAGTTGCCAGGTCCGACCGGAGCAGGCTGCGGCA‑‑GGGGGGGCGCTTTTCCCCTCACCCTAACCCTCTCCCCAGAGGGGCGAGGGGACCGTATTGTGCAAATATTGTTACCCCAGCAACAAACAGGCTCATACAGCCCCTAACCCTTTCATGGCGATGGCTGGGACGGTT  >  NC_000913/1351041‑1351324

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