Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,668,183 T→A 7.5% T336S (ACG→TCG) ‡ yhjA ← putative cytochrome C peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,668,1830TA 7.5% 183.6 / 4.9 53T336S (ACG→TCG) ‡yhjAputative cytochrome C peroxidase
Reads supporting (aligned to +/- strand):  ref base T (25/24);  new base A (2/2);  total (27/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.88e-01

AAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTT  >  NC_000913/3668067‑3668321
                                                                                                                    |                                                                                                                                          
aaaGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCAt                                                                                                                >  1:1179394/1‑145 (MQ=255)
aaaGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCAt                                                                                                                <  2:1179394/145‑1 (MQ=255)
    tCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCaaa                                                                                                        >  1:451475/1‑149 (MQ=255)
      ttttttCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCAtttt                                                                                                                                 >  4:725396/1‑122 (MQ=255)
      ttttttCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCAtttt                                                                                                                                 <  3:725396/122‑1 (MQ=255)
       tttttCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATgg                                                                                                     <  2:451475/149‑1 (MQ=255)
            cAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCAtttt                                                                                                                                 >  3:116890/1‑116 (MQ=255)
            cAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCAtttt                                                                                                                                 <  4:116890/116‑1 (MQ=255)
                cccAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAgc                                                                                                                                        >  2:790018/1‑105 (MQ=255)
                cccAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAgc                                                                                                                                        <  1:790018/105‑1 (MQ=255)
                 ccAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGc                                                                                           >  3:462125/1‑149 (MQ=255)
                               gAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATggg                                                                                                    <  3:171767/126‑1 (MQ=255)
                               gAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATggg                                                                                                    >  4:171767/1‑126 (MQ=255)
                               gAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGtt                                                                               >  1:1025140/1‑147 (MQ=255)
                               gAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGtt                                                                               <  2:1025140/147‑1 (MQ=255)
                                  cGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtctc                                                                          <  4:462125/149‑1 (MQ=255)
                                   gTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTa                                                                                                             <  4:1407973/113‑1 (MQ=255)
                                   gTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTa                                                                                                             >  3:1407973/1‑113 (MQ=255)
                                   gTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCa                                                                         >  2:318668/1‑149 (MQ=255)
                                             aataatACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCaa                                                               <  1:318668/149‑1 (MQ=255)
                                                 ataCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtc                                                                            >  1:29294/1‑132 (MQ=255)
                                                 ataCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtc                                                                            <  2:29294/132‑1 (MQ=255)
                                                     caccaTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtct                                                                           >  3:315121/1‑129 (MQ=255)
                                                     caccaTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtct                                                                           <  4:315121/129‑1 (MQ=255)
                                                        caTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCcgc                                                                                                                        <  3:102571/81‑1 (MQ=255)
                                                        caTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCcgc                                                                                                                        >  4:102571/1‑81 (MQ=255)
                                                            aCAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCcgc                                                                                                                        <  3:172017/77‑1 (MQ=255)
                                                            aCAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCcgc                                                                                                                        >  4:172017/1‑77 (MQ=255)
                                                              aaGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCaaa                                                                       <  1:1180225/124‑1 (MQ=255)
                                                              aaGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCaaa                                                                       >  2:1180225/1‑124 (MQ=255)
                                                                 ttGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCAtt                                                                                                                                   >  4:772470/1‑61 (MQ=255)
                                                                 ttGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCAtt                                                                                                                                   <  3:772470/61‑1 (MQ=255)
                                                                            ttatCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGtaattaat                                                                                  <  4:298684/99‑1 (MQ=255)
                                                                            ttatCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGtaattaat                                                                                  >  3:298684/1‑99 (MQ=255)
                                                                                 tttAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGtaattaa                                                                                   <  1:1029312/93‑1 (MQ=255)
                                                                                 tttAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGtaattaa                                                                                   >  2:1029312/1‑93 (MQ=255)
                                                                                      ataatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCaaa                                                                                                        <  1:1009174/67‑1 (MQ=255)
                                                                                      ataatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCaaa                                                                                                        >  2:1009174/1‑67 (MQ=255)
                                                                                      ataatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGc                            >  2:1396374/1‑143 (MQ=255)
                                                                                      ataatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGc                            <  1:1396374/143‑1 (MQ=255)
                                                                                         atTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGAt                        <  2:214482/144‑1 (MQ=255)
                                                                                         atTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGAt                        >  1:214482/1‑144 (MQ=255)
                                                                                                 ccTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGttt                                                                              >  2:987468/1‑82 (MQ=255)
                                                                                                 ccTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGttt                                                                              <  1:987468/82‑1 (MQ=255)
                                                                                                  cTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAAc          >  3:1339572/1‑149 (MQ=255)
                                                                                                   ttttttCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtc                                                                            >  2:811156/1‑82 (MQ=255)
                                                                                                   ttttttCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtc                                                                            <  1:811156/82‑1 (MQ=255)
                                                                                                   ttttttCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtc                                                                            <  2:1306862/82‑1 (MQ=255)
                                                                                                   ttttttCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtc                                                                            >  1:1306862/1‑82 (MQ=255)
                                                                                                          aaatGCGCTGACGGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGtttt  >  3:1370186/4‑149 (MQ=255)
                                                                                                           aatGCGCTGACGGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTc                                            >  1:24505/3‑106 (MQ=255)
                                                                                                           aatGCGCTGACGGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTc                                            <  2:24505/104‑1 (MQ=255)
                                                                                                           aatGCGCTGACGGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGaaa                                  <  2:644159/114‑1 (MQ=255)
                                                                                                                    tCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATggg                                                                                                    >  3:509694/1‑41 (MQ=255)
                                                                                                                    tCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATggg                                                                                                    <  4:509694/41‑1 (MQ=255)
                                                                                                                    |                                                                                                                                          
AAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTT  >  NC_000913/3668067‑3668321

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: