Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,701,211 C→A 8.7% T108T (ACG→ACT)  PP_3266 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,701,2110CA 8.7% 131.5 / 2.7 46T108T (ACG→ACT) PP_3266hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/18);  new base A (2/2);  total (26/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 = 7.15e-01

GAACTGCCACAGGAACGGCCCGAGTTTTTCGCCCAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCTTTGGGGTAGAAGTCCTTGCGCCAGGGGC  >  NC_002947/3701064‑3701359
                                                                                                                                                   |                                                                                                                                                    
gAACTGCCACAGGAACGGCCCGAGTTTTTCGCCCAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGt                                                                                                                                                    <  2:240342/150‑1 (MQ=255)
   cTGCCACAGGAACGGCCCGAGTTTTTCGCCCAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGc                                                                                                                                                 >  1:715360/1‑150 (MQ=255)
        acaGGAACGGCCCGAGTTTTTCGCCCAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCAc                                                                                                                                            <  2:715360/150‑1 (MQ=255)
          aGGAACGGCCCGAGTTTTTCGCCCAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGc                                                                                                                                          <  2:521122/150‑1 (MQ=255)
           ggAACGGCCCGAGTTTTTCGCCCAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGcc                                                                                                                                         <  1:770411/150‑1 (MQ=255)
             aaCGGCCCGAGTTTTTCGCCCAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATAGAGCGGGGCCCCTTGACTGCGAACACGAAGGCTTCGCGCGTTGGTTCACGCCAg                                                                                                                                       >  1:813050/1‑150 (MQ=255)
              aCGGCCCGAGTTTTTCGCCCAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGc                                                                                                                                      >  1:722751/1‑150 (MQ=255)
                   ccGAGTTTTTCGCCCAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGa                                                                                                                                 <  2:722751/150‑1 (MQ=255)
                    cGAGTTTTTCGCCCAACAGCAGCGGGCCGGAAGCGAAAAAATTTGCCACGTGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGAt                                                                                                                                >  1:578325/1‑150 (MQ=255)
                               cccAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGc                                                                                                                     <  1:796024/150‑1 (MQ=255)
                               cccAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGc                                                                                                                     >  1:547239/1‑150 (MQ=255)
                                ccAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGGTTGCGCACGGGCg                                                                                                                    >  1:40292/1‑150 (MQ=255)
                                ccAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCg                                                                                                                    >  1:39934/1‑150 (MQ=255)
                                  aaCAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCTCTCGGGCGTc                                                                                                                  >  2:748201/1‑150 (MQ=255)
                                                ggAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGaa                                                                                                    <  1:249669/150‑1 (MQ=255)
                                                   aGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGa                                                                                                 <  2:547239/150‑1 (MQ=255)
                                                        aaaaaTTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCAATCAAGCGTTAGTTCGGGCGATCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCAt                                                                                            >  1:829431/1‑150 (MQ=255)
                                                               ggCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGGTATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGAtctc                                                                                     <  2:614751/150‑1 (MQ=255)
                                                                           ttCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCAc                                                                         >  2:111103/1‑150 (MQ=255)
                                                                                     ttCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGccc                                                                                                                                    >  1:104039/1‑81 (MQ=255)
                                                                                     ttCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGccc                                                                                                                                    <  2:104039/81‑1 (MQ=255)
                                                                                          ccgccgGACATGGGTGATATAGCGCGTCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCGGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGcaga                                                          >  1:49746/1‑147 (MQ=255)
                                                                                                   aTGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCAtt                                                                                           <  1:46911/108‑1 (MQ=255)
                                                                                                   aTGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCAtt                                                                                           >  2:46911/1‑108 (MQ=255)
                                                                                                                   ggCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGaa                                                                                                    >  1:76011/1‑83 (MQ=255)
                                                                                                                   ggCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGaa                                                                                                    <  2:76011/83‑1 (MQ=255)
                                                                                                                     ccccTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGCTAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGCTTCAGccc                               >  2:827116/1‑150 (MQ=255)
                                                                                                                      cccTTGACTGCGAACGCGAAGCGTTCCGGCGTTGCTTCACGCCAGCCGAGATAGCGCTCGGGCGTCTGCCGGCTGTAGAACGACCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCt                              <  2:40292/150‑1 (MQ=255)
                                                                                                                      cccTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCt                              <  2:39934/150‑1 (MQ=255)
                                                                                                                       ccTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCtt                             >  2:462265/1‑150 (MQ=255)
                                                                                                                       ccTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCtt                             >  2:499029/1‑150 (MQ=255)
                                                                                                                       ccTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCtt                             >  2:111985/1‑150 (MQ=255)
                                                                                                                         ttGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCTTTg                           <  1:462265/150‑1 (MQ=255)
                                                                                                                               gCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCGTTGGGGtcg                     >  1:132126/1‑148 (MQ=255)
                                                                                                                                cGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTg                                                                                          <  2:286137/80‑1 (MQ=255)
                                                                                                                                cGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTg                                                                                          >  1:286137/1‑80 (MQ=255)
                                                                                                                                  aaCACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCTTTGGGGTAGAAg                  >  1:134209/1‑150 (MQ=255)
                                                                                                                                     acGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCTTTGGGGTAGAAGTcc               >  2:808764/1‑150 (MQ=255)
                                                                                                                                       gAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCTTTGGGGTAGAAGTCCtt             >  2:594538/1‑150 (MQ=255)
                                                                                                                                         aGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCTTTGGGGTAGAAGTCCTTgc           >  1:682569/1‑150 (MQ=255)
                                                                                                                                           ccTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCTTTGGGGTAGAAGTCCTTgcgc         <  2:829431/150‑1 (MQ=255)
                                                                                                                                            cTTCCGGAGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGccc                                                                    >  2:122729/1‑90 (MQ=255)
                                                                                                                                            cTTCCGGAGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGccc                                                                    >  2:122685/1‑90 (MQ=255)
                                                                                                                                            cTTCCGGAGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGccc                                                                    <  1:122729/90‑1 (MQ=255)
                                                                                                                                            cTTCCGGAGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGccc                                                                    <  1:122685/90‑1 (MQ=255)
                                                                                                                                                  gCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCTTTGGGGTAGAAGTCCTTGCGCCAGGGGc  <  1:111985/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GAACTGCCACAGGAACGGCCCGAGTTTTTCGCCCAACAGCAGCGGGCCGGAAGCGAAAAAATTGGCCACGGGCTCTTCGATTTCCTTCAACCGCCGGACATGGGTGATATAGCGCGGCCCCTTGACTGCGAACACGAAGCCTTCCGGCGTTGCTTCACGCCAGCCCAGATAGCGCTCGGGCGTCTGCAGGCTGTAGAACGAGCCATTGATCTCGATGCTGTTCACCGCCCGTGAAGCGAATGCCAGTTCATCGTCCTGTTTCAGGCCTTTGGGGTAGAAGTCCTTGCGCCAGGGGC  >  NC_002947/3701064‑3701359

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