Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,250,989 G→C 13.8% intergenic (+35/‑35) yqjF → / → yqjG putative quinol oxidase subunit/putative S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,250,9890GC13.8% 106.7 / 7.6 36intergenic (+35/‑35)yqjF/yqjGputative quinol oxidase subunit/putative S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (14/17);  new base C (2/3);  total (16/20)
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.91e-01

TGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAAC  >  NC_000913/3250854‑3251119
                                                                                                                                       |                                                                                                                                  
tGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGtt                                                                                                                        <  1:455372/148‑1 (MQ=255)
        aaaaaCCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGtt                                                                                                                        <  1:160012/140‑1 (MQ=255)
                  cAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGgaggag                                                                                                        <  3:337055/146‑1 (MQ=255)
                  cAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGgaggag                                                                                                        >  4:337055/1‑146 (MQ=255)
                    aTTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGtt                                                                                                                        <  4:799671/128‑1 (MQ=255)
                    aTTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGtt                                                                                                                        >  3:799671/1‑128 (MQ=255)
                                  ccTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTg                                                                                      <  2:264323/148‑1 (MQ=255)
                                          tgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAAt                                                                                                <  1:316211/130‑1 (MQ=255)
                                          tgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAAt                                                                                                >  2:316211/1‑130 (MQ=255)
                                                        ccGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGa                                                                                 <  1:101296/131‑1 (MQ=255)
                                                        ccGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGa                                                                                 >  2:101296/1‑131 (MQ=255)
                                                            gcgcGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGa                                                                                     <  3:466985/123‑1 (MQ=255)
                                                            gcgcGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGa                                                                                     >  4:466985/1‑123 (MQ=255)
                                                                       aTCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGa                                                                                     >  2:835367/1‑112 (MQ=255)
                                                                       aTCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGa                                                                                     <  1:835367/112‑1 (MQ=255)
                                                                                  gctgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATggg                                                                                             <  3:879779/93‑1 (MQ=255)
                                                                                  gctgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATggg                                                                                             >  4:879779/1‑93 (MQ=255)
                                                                                     gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAt                                                                                                                >  3:738040/1‑71 (MQ=255)
                                                                                     gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAt                                                                                                                <  4:738040/71‑1 (MQ=255)
                                                                                     gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATAc                                                   >  1:1053152/1‑132 (MQ=255)
                                                                                     gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATAc                                                   <  2:1053152/132‑1 (MQ=255)
                                                                                      gtcaaccagttgacccatttttctctcctccaatagcAAAAACGAGGAGGAAGCACCTCGtt                                                                                                                        <  2:463285‑M1/25‑1 (MQ=255)
                                                                                          atcagttgacccatttttctctcctccaatagcAAAAACGAGGAGGAAGCTCCTCg                                                                                                                          >  3:159468‑M1/34‑56 (MQ=255)
                                                                                          atcagttgacccatttttctctcctccaatagcAAAAACGAGGAGGAAGCTCCTCg                                                                                                                          <  4:159468‑M1/23‑1 (MQ=255)
                                                                                           aaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGctcct                                                                                                                            >  2:830219/1‑53 (MQ=255)
                                                                                           aaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGctcct                                                                                                                            <  1:830219/53‑1 (MQ=255)
                                                                                               ttgacccatttttctctcctccaatagcAAAAACGAGGAGGAAGCTCCTCg                                                                                                                          <  3:680174‑M1/23‑1 (MQ=255)
                                                                                               ttgacccatttttctctcctccaatagcAAAAACGAGGAGGAAGCTCCTCg                                                                                                                          >  4:680174‑M1/29‑51 (MQ=255)
                                                                                                   aaGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTAcc                                         <  2:629032/128‑1 (MQ=255)
                                                                                                   aaGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTAcc                                         >  1:629032/1‑128 (MQ=255)
                                                                                                    aGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTgga                                                                                                            <  2:439589/60‑1 (MQ=255)
                                                                                                    aGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTgga                                                                                                            >  1:439589/1‑60 (MQ=255)
                                                                                                       aCGCGCTATACTTAACGATAAAAAACGAGGAGCATCCTCCTCg                                                                                                                          >  2:178218/1‑43 (MQ=255)
                                                                                                                   tAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGt     >  2:288598/1‑148 (MQ=255)
                                                                                                                      cGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAAc  <  1:288598/148‑1 (MQ=255)
                                                                                                                              AACGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGACAACTGATTGAc                                                                                <  1:463285‑M1/62‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
TGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAAC  >  NC_000913/3250854‑3251119

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