Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,769,912 G→C 14.4% A89G (GCC→GGC)  PP_2422 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,769,9120GC14.4% 84.6 / 5.1 35A89G (GCC→GGC) PP_2422hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (12/16);  new base C (3/2);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AGCCGCTCAGTCTGGGTGGCGCGGCGCCGTCGCGCAAGGGCGGCGACGCCGTTTGTGCTCAATCGTTGGCCTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAGGAGATCAGCTCGCGGGTCTTGGCATCCAGGTGGTTGGTCTTGGCACCT  >  NC_002947/2769766‑2770058
                                                                                                                                                  |                                                                                                                                                  
aGCCGCTCAGTCTGGGTGGCGCGGCGCCGTCGCGCAAGGGCGGCGACGCCGTTTGTGCTCAATCGTTGGCCTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCt                                                                                                                                                 >  2:592430/1‑150 (MQ=255)
            tGGGTGGCGCGGCGCCGTCGCGCAAGGGCGGCGACGCCGTTTGTGCTCAATCGTTGGCCTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCt                                                                                                                                     >  1:233345/1‑150 (MQ=255)
                    gcggcgcCGTCGCGCAAGGGCGGCGACGCCGTTTGTGCTCAATCGTTGGCCTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTg                                                                                                                             <  2:600605/150‑1 (MQ=255)
                    gcggcgcCGTCGCGCAAGGGCGGCGACGCCGTTTGTGCTCAATCGTTGGCCTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTg                                                                                                                             <  1:592430/150‑1 (MQ=255)
                              cgcgCAAGGGCGGCGACGCCGTTTGTGCTCAATCGTTGGCCTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCtt                                                                                                                   <  2:371394/150‑1 (MQ=255)
                                 gcAAGGGCGGCGACGCCGTTTGTGCTCAATCGTTGGCCTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGa                                                                                                                <  1:328094/150‑1 (MQ=255)
                                   aaGGGCGGCGACGCCGTTTGTGCTCAATCGTTGGCCTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACg                                                                                                              <  2:433043/150‑1 (MQ=255)
                                                   tttGTGCTCAATCGTTGGCCTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCa                                                                                              <  2:599673/150‑1 (MQ=255)
                                                    ttGTGCACAGGCGTTGGCCGTGCCGAGCGCGTCCATTGCACGGCCGCAGTACCCCAGCGCGGCACCGGAGTTGATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCAc                                                                                             <  1:168009/150‑1 (MQ=255)
                                                          tCAATCGTTGGCCTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGAt                                                                                       >  2:293834/1‑150 (MQ=255)
                                                          tCAATCGTTGGCCTTGCCCACCGCGTCCATTGCACGAGCGGAATACACAAGCGCGGCACCGGCGTTCATGGCCACGGCAACGCCGAGGCCCTCGGCAAGTTCCTCGCGGCTCGCCCCCGCTTTGACCGCCTGCTACGCACGCTCGGCgca                                                                                       >  2:98100/1‑148 (MQ=255)
                                                            aaTCGTTGGCCTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGACTGTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGc                                                                                     >  2:391857/1‑150 (MQ=255)
                                                                     ccTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGc                                                                            >  1:292650/1‑150 (MQ=255)
                                                                           ccACCGCGTCCATTGCACGGGCGGAGTACACTAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCgag                                                                      >  1:308188/1‑148 (MQ=255)
                                                                                  gTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGggtg                                                                    >  2:10429/1‑145 (MQ=255)
                                                                                  gTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGggtg                                                                    <  1:10429/145‑1 (MQ=255)
                                                                                   tCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCAGGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCCGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTcacg                                                              <  1:98100/150‑2 (MQ=255)
                                                                                   tCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTcaca                                                              >  1:92963/1‑150 (MQ=255)
                                                                                   tCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTcaca                                                              >  1:470802/1‑150 (MQ=255)
                                                                                        tGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGACGAGGTCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCAc                                                         <  1:409554/150‑1 (MQ=255)
                                                                                               gCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAg                                                  >  1:362438/1‑150 (MQ=255)
                                                                                               gCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAg                                                  <  2:233345/150‑1 (MQ=255)
                                                                                                       cacCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAGGAGATCAg                                          <  2:292650/150‑1 (MQ=255)
                                                                                                        acCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGAAACTCGCCTATCTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCCGCGGGTGCTCACAGCCACGGCCAAGgagaaccgc                                         >  1:108631/1‑145 (MQ=255)
                                                                                                                   aCCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAGGAGATCAGCTCGCGGGTCtt                              >  2:482994/1‑150 (MQ=255)
                                                                                                                    ccGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAGGAGATCAGCTCGCGGGTCTTg                             <  2:129032/150‑1 (MQ=255)
                                                                                                                               ggccacggccacGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAGGAGATCAGCTCGCGGGTCTTGGCATCCAggtg                  <  2:92963/150‑1 (MQ=255)
                                                                                                                               ggccacggccacGCCGAGGGCCTAGCCAGGTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAGGAGATCAGCTCGCGGGTCTTGGCATCCAggtg                  >  1:478971/1‑150 (MQ=255)
                                                                                                                                         ttgCCGAGGCCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTcacc                                                              <  2:375718/94‑2 (MQ=255)
                                                                                                                                         ttgCCGAGGCCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTcacc                                                              >  1:375718/3‑95 (MQ=255)
                                                                                                                                         ttgCCGAGGCCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAGGAGATCAGCTCGCGGGTCTTg                             >  2:127141/3‑129 (MQ=255)
                                                                                                                                         ttgCCGAGGCCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAGGAGATCAGCTCGCGGGTCTTg                             <  1:127141/127‑1 (MQ=255)
                                                                                                                                           gCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAGGAGATCAGCTCGCGGGTCTTGGCATCCAGGTGGTTGGTCTTGGc      <  1:482994/150‑1 (MQ=255)
                                                                                                                                             ggAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCGCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGGCGCAGCGGGTGGTCACAGCCACGGCCAAGGAGATCAGCTCGCGGGTCTTGGCATCCAGGTGGTTGGTCTTGGCAc    <  2:108631/149‑1 (MQ=255)
                                                                                                                                               aGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAGGAGATCAGCTCGCGGGTCTTGGCATCCAGGTGGTTGGTCTTGGCACCt  <  2:308188/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGCCGCTCAGTCTGGGTGGCGCGGCGCCGTCGCGCAAGGGCGGCGACGCCGTTTGTGCTCAATCGTTGGCCTTGCCCACCGCGTCCATTGCACGGGCGGAGTACACCAGCGCGGCACCGGCGTTCATGGCCACGGCCACGCCGAGGGCCTCGGCAATTTCCTCGCGGCTGGCCCCCGCTTTGACGGCCTGCTGCGAATGCACGGCGATGCAGCCGTCGCAGCGGGTGGTCACAGCCACGGCCAAGGAGATCAGCTCGCGGGTCTTGGCATCCAGGTGGTTGGTCTTGGCACCT  >  NC_002947/2769766‑2770058

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