Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,624,901 C→T 11.7% L210L (CTG→CTA) ‡ PP_3195 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,624,9010CT11.7% 101.9 / 6.0 43L210L (CTG→CTA) ‡PP_3195hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/19);  new base T (3/2);  total (22/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 = 8.58e-01

TCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGCCTTCATACGGGATCTCTAC  >  NC_002947/3624764‑3625047
                                                                                                                                         |                                                                                                                                                   
tCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCT‑GCCTGGCTGGCcaggcgaggc                                                                                                                                        <  2:1040643/150‑6 (MQ=255)
  aGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggcaggcag                                                                                                                                       <  1:422815/150‑1 (MQ=255)
    cACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCaggcaggcaggc                                                                                                                                     >  2:783604/1‑150 (MQ=255)
    cACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCaggcaggcaggc                                                                                                                                     <  2:604566/150‑1 (MQ=255)
               agccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTat                                                                                                                          >  1:556876/1‑150 (MQ=255)
                gccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatt                                                                                                                         >  2:139734/1‑150 (MQ=255)
                        aCTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTgtgtcagtgcc                                                                                                                 >  2:812700/1‑145 (MQ=255)
                          tATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCtt                                                                                                               >  2:1420788/1‑150 (MQ=255)
                           aTGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTg                                                                                                              <  1:1420788/150‑1 (MQ=255)
                              cGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGtt                                                                                                           <  2:1162332/150‑1 (MQ=255)
                                  tGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCTAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAAt                                                                                                       <  2:416454/150‑1 (MQ=255)
                                     cccAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGACTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCa                                                                                                    <  1:812700/150‑1 (MQ=255)
                                            ttCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGtt                                                                                             <  2:158741/150‑1 (MQ=255)
                                            ttCGCTGTCATGGCGTCCGATCTGGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGtt                                                                                             <  1:1142410/150‑1 (MQ=255)
                                                 tGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACtt                                                                                        <  1:1004091/150‑1 (MQ=255)
                                                   tCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcagg                                                                                                                                              >  1:656533/1‑94 (MQ=255)
                                                   tCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcagg                                                                                                                                              <  2:656533/94‑1 (MQ=255)
                                                     aTGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCTGGCcaggcga                                                                                                                                           >  2:339429/1‑93 (MQ=255)
                                                     aTGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCTGGCcaggcga                                                                                                                                           <  1:339429/95‑3 (MQ=255)
                                                             ccGTCATGCCGTGGAGGCGAAACGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCg                                                                           <  2:1056516/150‑1 (MQ=255)
                                                               gTCATGCCGTGGAGGCAAAGCGCTTCGCCGGTGTCGGGATGGTCGACAATCAGTGTCGACACGCCTCGCCTGGTCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGTCTTG‑CACCAGCAGCGGt                                                                         <  2:870612/150‑1 (MQ=255)
                                                                tCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGAcgcg                                                                                                                              >  1:829936/1‑97 (MQ=255)
                                                                tCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGAcgcg                                                                                                                              <  2:829936/97‑1 (MQ=255)
                                                                tCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTg                                                                        <  1:783604/150‑1 (MQ=255)
                                                                                              ttccGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACca                                          <  1:139734/148‑1 (MQ=255)
                                                                                                 cGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACcagca                                       >  1:316368/1‑150 (MQ=255)
                                                                                                 cGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACcagca                                       >  1:1392955/1‑150 (MQ=255)
                                                                                                      tGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTccc                                  >  1:1338457/1‑150 (MQ=255)
                                                                                                        gTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCga                                >  2:742059/1‑150 (MQ=255)
                                                                                                        gTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCga                                <  1:1093895/150‑1 (MQ=255)
                                                                                                         tCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGGATCTAGCCCGTTGACTTGCCATGATAAACCAT‑CCGGCCCACTTACGCCTTCTGCGCGCACCAGCCGTCCCgag                               >  1:319203/1‑150 (MQ=255)
                                                                                                           gACAATCAGCGACGACACGCCCCAGCTGGCCACCCAAGCACGCAGGCCGATGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGagag                             >  2:728800/1‑150 (MQ=255)
                                                                                                             cAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGt                           >  1:93045/1‑150 (MQ=255)
                                                                                                                 cAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCtt                       >  2:1243421/1‑150 (MQ=255)
                                                                                                                   gCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGc                     >  1:973977/1‑150 (MQ=255)
                                                                                                                   gCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGc                     <  1:1312394/150‑1 (MQ=255)
                                                                                                                        gacACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGCCTTCa                >  2:1248421/1‑150 (MQ=255)
                                                                                                                             gcctcgcctGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGCCTTCATACgg           >  2:263579/1‑150 (MQ=255)
                                                                                                                             gcctcgcctGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGCCTTCATACgg           <  2:93045/150‑1 (MQ=255)
                                                                                                                               cctgcctGGCTGGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCCGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGCCTTCATCCGGGa         >  2:1294300/4‑150 (MQ=255)
                                                                                                                                ctgcctGGCTGGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGCCTTCATACGGGAt        >  1:1352522/3‑150 (MQ=255)
                                                                                                                                  gcctGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGCCTTCATACGGGAtct      >  2:524633/1‑150 (MQ=255)
                                                                                                                                      ggCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTCCTTGCACCAGCAGTCCCGAGAGGTGCTTGCCTTCATAGGGGATGTCTAc  <  2:556876/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
TCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTG‑CACCAGCAGCGGTGCCGGGCCACTCACGCCTTCTGCGCGCACCAGCAGTCCCGAGAGGTGCTTGCCTTCATACGGGATCTCTAC  >  NC_002947/3624764‑3625047

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