Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,414,748 C→A 10.8% intergenic (+290/‑55) PP_2114 → / → PP_2115 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,414,7480CA10.8% 82.6 / 4.6 37intergenic (+290/‑55)PP_2114/PP_2115transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (14/19);  new base A (2/2);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

ATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAG  >  NC_002947/2414600‑2414892
                                                                                                                                                    |                                                                                                                                                
aTTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTccc                                                                                                                                                  >  7:73031/1‑149 (MQ=38)
      gATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAg                                                                                                                                             <  1:140245/148‑1 (MQ=34)
      gATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAg                                                                                                                                             >  3:215458/1‑148 (MQ=34)
      gATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCAGCGCTGCGCACCGGGGTCCCCTCACAACCGCCTTGCTCCCGGCAg                                                                                                                                             >  2:140245/1‑148 (MQ=11)
          gCGACACTGACTCTCCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACcg                                                                                                                                        >  6:167268/1‑149 (MQ=21)
                   aCTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAgg                                                                                                                               >  8:181662/1‑149 (MQ=32)
                   aCTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAg                                                                                                                                <  1:300889/148‑1 (MQ=32)
                    cTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGTGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGc                                                                                                                              <  1:13074/149‑1 (MQ=17)
                       aCCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGcccc                                                                                                                           <  8:271023/149‑1 (MQ=32)
                       aCCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGcc                                                                                                                             <  7:44548/147‑1 (MQ=32)
                         cTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCa                                                                                                                                              <  5:119298/128‑1 (MQ=32)
                         cTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCa                                                                                                                                              >  6:119298/1‑128 (MQ=32)
                             tCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCt                                                                                                                      <  8:162605/148‑1 (MQ=31)
                                    ctattgctgtagctattgctgttgctgttgctgttgatcttgatcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCt                                                                                                               >  7:303448‑M2/110‑148 (MQ=255)
                                                                  ctgttgatcttgatcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAcc                                                                                <  8:303448‑M2/70‑1 (MQ=255)
                                                                               tcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcGAAGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAc                                                                                 <  2:211419‑M2/68‑1 (MQ=255)
                                                                               tcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcGAAGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAc                                                                                 <  8:211186‑M2/68‑1 (MQ=255)
                                                                               tcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcGAAGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAc                                                                                 >  1:211419‑M2/67‑134 (MQ=255)
                                                                               tcttgatcttgatcttgatcttgatcttgatcttcgcgacttcaggaggccgaacgcaggccttgcGAAGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAc                                                                                 >  7:211186‑M2/67‑134 (MQ=255)
                                                                                aaaCGCTGCGCTTCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGc                                                                  <  7:176980/149‑1 (MQ=17)
                                                                                              catcatCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGcctc                                                    <  5:65522/149‑1 (MQ=32)
                                                                                                        ccccTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTcgc                                          <  8:73031/149‑1 (MQ=32)
                                                                                                         cccTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCGCAGCGCTTGACCGGCCTCCAGGCCCGTCCCCTGCCGTCGCCAGGCCGGCGTTCGGCCGCGGGAAGTcgcg                                         >  7:240339/1‑149 (MQ=11)
                                                                                                                     ctCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatc                             <  5:167268/149‑1 (MQ=32)
                                                                                                                        cGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaag                          <  5:163660/149‑1 (MQ=32)
                                                                                                                        cGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaa                           <  8:271203/148‑1 (MQ=32)
                                                                                                                          gggTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTg                                                 <  6:312479/124‑1 (MQ=32)
                                                                                                                          gggTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTg                                                 >  5:312479/1‑124 (MQ=32)
                                                                                                                             tCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGTGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaa                                 <  5:101962/137‑1 (MQ=17)
                                                                                                                               cccTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatc                       >  5:6410/1‑145 (MQ=32)
                                                                                                                               cccTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatc                       >  5:137733/1‑145 (MQ=32)
                                                                                                                                 cTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatc                 <  3:224241/149‑1 (MQ=32)
                                                                                                                                        cGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGcctcc                                                   >  7:21132/1‑108 (MQ=32)
                                                                                                                                        cGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGcctcc                                                   <  8:21132/108‑1 (MQ=32)
                                                                                                                                         cgttctgcTCCCGGCAGGACCGCGCCTAAGGCCCCATCCTGGGGCCTAAGCGCTTGACGGGCATCCATGCCCGTCACTTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTGGCTatgatcaagatcaagatcaagatcaagatcaa         <  1:19671‑M2/141‑1 (MQ=255)
                                                                                                                                          gccttgcTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCGGAAGTCGCGaagatcaagatcaagatcaagatcaagatc           >  1:131481‑M2/8‑146 (MQ=255)
                                                                                                                                              tgcTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCAGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatcaagatcaagatcaagatc     >  6:94318‑M2/4‑148 (MQ=255)
                                                                                                                                                cTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAACtcgcgaagatcaagatcaagatcaagatcaagatcaagagaacg  >  2:19671‑M2/2‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ATTCGCGATGGCGACACTGACTCACCTTTTCGCCCTTACGGCGAGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGCGCTCCCATCATCCGGCCCCTGCGCTGCGCTCCGGGGTCCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTTCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAGATCAAGATCAAGATCAAG  >  NC_002947/2414600‑2414892

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