Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,762,985 A→T 16.7% N62I (AAC→ATC)  PP_3326 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,762,9850AT16.7% 95.8 / 14.3 42N62I (AAC→ATC) PP_3326hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/18);  new base T (2/5);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

CCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTC  >  NC_002947/3762836‑3763133
                                                                                                                                                      |                                                                                                                                                    
ccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGaa                                                                                                                                                      <  2:764443/150‑1 (MQ=255)
       gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGGGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTAAAAGTCGTAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc                                                                                                                                               >  1:692271/1‑150 (MQ=255)
          aGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTCCGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTCGGCGTTCGA‑CGTGACCGGGAACATGagctcg                                                                                                                                            >  1:47497/1‑145 (MQ=255)
          aGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTccc                                                                                                                                             <  1:407991/150‑1 (MQ=255)
           gCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCgg                                                                                                                                           <  2:249055/150‑1 (MQ=255)
              tCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCC‑AGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCa                                                                                                                                        <  2:86615/150‑1 (MQ=38)
                      ggCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCACgtcgaa                                                                                                                                 <  2:6298/150‑6 (MQ=38)
                         gcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAgcgc                                                                                                                             >  2:35246/1‑150 (MQ=255)
                           ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGATCATg                                                                                                                                                  >  1:367685/1‑127 (MQ=255)
                           ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGATCATg                                                                                                                                                  <  2:367685/127‑1 (MQ=255)
                               ggCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCAGGGCGTTCGA‑CGTGACCGGGAACACGATCCCGGTCACGGTCGAGCGCCTGGCg                                                                                                                       >  1:814943/1‑150 (MQ=255)
                                        aCCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGcc                                                                                                              >  1:774651/1‑150 (MQ=255)
                                          cGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACgg                                                                                                                                     <  2:634102/125‑1 (MQ=255)
                                          cGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACgg                                                                                                                                     >  1:634102/1‑125 (MQ=255)
                                             ggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtgct                                                                                                         <  2:774651/150‑1 (MQ=255)
                                                 cGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCaa                                                                                                     >  1:211477/1‑150 (MQ=255)
                                                    cgcggCAATATCAAGACATCGTAAGGCGATGGCGCCGAGGATAACGCCTACAATGCGGAAAAGATCAATCGCGCCTGGGCGGTCGA‑CGTGACCGATAACATGATCCCGGGCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCa                                                                                                  <  2:47497/147‑1 (MQ=255)
                                                     gTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCAc                                                                                                 >  1:636604/1‑150 (MQ=255)
                                                            tatCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAcccc                                                                                          <  2:636604/150‑1 (MQ=255)
                                                                   gCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCg                                                                                   <  1:35246/150‑1 (MQ=255)
                                                                            ggCCACGGCGCCGAGGATAACGCCTAAAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGcc                                                                          <  1:226130/150‑1 (MQ=255)
                                                                              ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGt                                                                        <  2:591271/150‑1 (MQ=255)
                                                                               cACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtg                                                                       >  1:41124/1‑150 (MQ=255)
                                                                               cACGGCGCCGAGGATAACGCCTACAAGGAGGAAAAGATCAAGCGCGCATGTGCGTGCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtg                                                                       >  1:322984/1‑150 (MQ=255)
                                                                                            aTAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACCAGCAGTTCTGCGGGACAACg                                                          >  2:794490/1‑150 (MQ=255)
                                                                                               aCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCACACACCTACCCTCCCGTCGGCACTTGCCGTGCTGCTGGACAACGGGc                                                       >  1:615570/1‑150 (MQ=255)
                                                                                               aCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGc                                                       >  2:102874/1‑150 (MQ=255)
                                                                                                     aCAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGt                                                 <  2:41124/150‑1 (MQ=255)
                                                                                                       aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGGCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTca                                               <  2:814943/150‑1 (MQ=255)
                                                                                                         ggCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                             <  1:102874/150‑1 (MQ=255)
                                                                                                            ggAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgcg                                          >  2:66444/1‑150 (MQ=255)
                                                                                                             gAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCaa                                                                                                     <  2:217857/90‑1 (MQ=255)
                                                                                                             gAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCaa                                                                                                     >  1:217857/1‑90 (MQ=255)
                                                                                                                  gATCAAGCGCGCCTGTGTGTTCTA‑CGTCACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGGTCGCCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGa                                    <  1:270551/150‑1 (MQ=255)
                                                                                                                      aaGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAAc                                                                                                    <  2:561404/82‑1 (MQ=255)
                                                                                                                      aaGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAAc                                                                                                    >  1:561404/1‑82 (MQ=255)
                                                                                                                          gcgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAAt                            <  2:211477/150‑1 (MQ=255)
                                                                                                                             ccTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATggg                         <  1:148406/150‑1 (MQ=255)
                                                                                                                              cTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt                        <  1:66444/150‑1 (MQ=255)
                                                                                                                                     ctcgaccGTGACCGGGATCATGTTCCCGGTCAC‑GTCGAACGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAAt                    <  2:390602/141‑1 (MQ=38)
                                                                                                                                     ctcgaccGTGACCGGGATCATGTTCCCGGTCAC‑GTCGAACGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAAt                    >  1:390602/7‑147 (MQ=38)
                                                                                                                                        gA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt              <  2:322984/150‑1 (MQ=255)
                                                                                                                                                 cGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGAGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCa      >  2:695807/1‑150 (MQ=255)
                                                                                                                                                     aaCATGATCCCGGTCACGGTCGCGCGCCTGGCGGATCATGCCGTGCTCCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGTTCAATCAGGATCGCTCCCATGTc  <  2:692271/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTC  >  NC_002947/3762836‑3763133

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: ‑