Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,411 C→A 8.3% M291I (ATG→ATT)  feaR ← transcriptional activator for tynA and feaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,4110CA 8.3% 163.5 / 5.9 48M291I (ATG→ATT) feaRtranscriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base C (27/17);  new base A (2/2);  total (29/19)
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.37e-02

CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCG  >  NC_000913/1446263‑1446559
                                                                                                                                                    |                                                                                                                                                    
cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                                                                                                                                      >  1:142465/1‑149 (MQ=255)
    aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATcc                                                                                                                                                  >  2:404480/1‑149 (MQ=255)
        aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaaa                                                                                                                                              <  4:499035/149‑1 (MQ=255)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGc                                                                                                                                            >  1:138577/1‑149 (MQ=255)
          gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGc                                                                                                                                            >  1:533491/1‑149 (MQ=255)
             aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                                                                                                                         >  1:9547/1‑149 (MQ=255)
              atatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTg                                                                                                                                        >  3:264186/1‑149 (MQ=255)
                           ccTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                                                                                                                           <  4:378283/149‑1 (MQ=255)
                               tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTAATCCCAAAGCGtt                                                                                                                                         >  4:362233/1‑131 (MQ=255)
                               tCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTAATCCCAAAGCGtt                                                                                                                                         <  3:362233/131‑1 (MQ=255)
                                 aaaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa                                                                                                                     >  4:72336/1‑149 (MQ=255)
                                         ttttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTg                                                                                                             <  3:7623/149‑1 (MQ=255)
                                                cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaa                                                                                                      >  2:493352/1‑149 (MQ=255)
                                                 tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaaa                                                                                                     >  2:282551/1‑149 (MQ=255)
                                                     aaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATcc                                                                                                 <  2:9547/149‑1 (MQ=255)
                                                     aaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTAATCCCAAAGCGtt                                                                                                                                         >  3:179448/1‑109 (MQ=255)
                                                     aaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTAATCCCAAAGCGtt                                                                                                                                         <  4:179448/109‑1 (MQ=255)
                                                                   atatCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGa                                                                                   >  2:277938/1‑149 (MQ=255)
                                                                          ttttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAtgcctgc                                                                            >  2:455369/1‑149 (MQ=255)
                                                                            ttAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCa                                                                          <  2:514551/149‑1 (MQ=255)
                                                                                 tAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGtttt                                                                     >  1:516606/1‑149 (MQ=255)
                                                                                   atcggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTtc                                                                   >  4:294355/4‑149 (MQ=255)
                                                                                      ggCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTtcatc                                                                >  2:129073/1‑149 (MQ=255)
                                                                                            gCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGc                                                          <  3:496442/149‑1 (MQ=255)
                                                                                             cAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCg                                                         <  3:138610/149‑1 (MQ=255)
                                                                                               gTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGc                                                       >  2:283366/1‑149 (MQ=255)
                                                                                               gTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGc                                                       >  1:498738/1‑149 (MQ=255)
                                                                                                   aTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATgg                                                   >  2:487288/1‑149 (MQ=255)
                                                                                                   aTGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATgg                                                   >  1:270022/1‑149 (MQ=255)
                                                                                                    tGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGc                                                  <  3:83773/149‑1 (MQ=255)
                                                                                                     gCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAtgcct                                                                              <  1:25961/120‑1 (MQ=255)
                                                                                                     gCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAtgcct                                                                              >  2:25961/1‑120 (MQ=255)
                                                                                                       aaCGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGaa                                               >  3:518061/1‑149 (MQ=255)
                                                                                                       aaCGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGaa                                               >  1:203566/1‑149 (MQ=255)
                                                                                                         cGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATc                                             <  2:433794/149‑1 (MQ=255)
                                                                                                         cGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATc                                             <  2:306741/149‑1 (MQ=255)
                                                                                                          gCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCg                                            >  3:399607/1‑149 (MQ=255)
                                                                                                               ttagttagCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCt                                       <  1:493352/149‑1 (MQ=255)
                                                                                                                tagttagCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTg                                      >  2:370960/1‑149 (MQ=255)
                                                                                                                   ttagCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGcac                                   <  2:68962/149‑1 (MQ=255)
                                                                                                                     agCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACaa                                 >  2:231172/1‑149 (MQ=255)
                                                                                                                      gCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACaaa                                >  4:188697/1‑149 (MQ=255)
                                                                                                                         gAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAAt                             >  4:492037/1‑149 (MQ=255)
                                                                                                                                 gTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAgacga                     >  2:71302/1‑149 (MQ=255)
                                                                                                                                   cGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGacgacg                   <  4:264186/149‑1 (MQ=255)
                                                                                                                                     aTACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAAACGCATCTGCACAAAAATCGAGACGACGGt                 >  1:496409/1‑149 (MQ=255)
                                                                                                                                        cTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTAc              <  3:416212/149‑1 (MQ=255)
                                                                                                                                            ccTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAat          <  2:533491/149‑1 (MQ=255)
                                                                                                                                                   tCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTgc   <  3:79163/149‑1 (MQ=255)
                                                                                                                                                    cATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTgcg  >  3:314531/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCAGTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGGCATGGCGAATCGCATCTGCACAAAAATCGAGACGACGGTTACGAATATATTGCG  >  NC_000913/1446263‑1446559

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