Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,838,124 C→T 7.7% A559V (GCC→GTC)  PP_3388 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,838,1240CT 7.7% 138.6 / 2.4 52A559V (GCC→GTC) PP_3388hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/22);  new base T (2/2);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

ACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACGGCTGTCGATATTCGCCTGAGTGAAGCTGCACG  >  NC_002947/3837984‑3838268
                                                                                                                                            |                                                                                                                                                 
aCCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGt                                                                                                                                          >  2:1020935/1‑150 (MQ=39)
 ccgccgTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTc                                                                                                                                         <  1:1020935/150‑1 (MQ=39)
  cgccgTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTcg                                                                                                                                        >  2:209068/1‑150 (MQ=38)
       tCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTg                                                                                                                                   <  2:999778/150‑1 (MQ=255)
                  cTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGt                                                                                                                        >  1:1142591/1‑150 (MQ=255)
                              ggAAGCGCTTGCCAGGCAAGCCCATGGCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGc                                                                                                            >  1:586501/1‑150 (MQ=255)
                                   cgcTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCAccc                                                                                                       >  1:373729/1‑150 (MQ=255)
                                           aGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGc                                                                                              >  1:781680/1‑150 (MQ=255)
                                             gCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCAt                                                                                            <  2:373729/150‑1 (MQ=255)
                                               aaGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATgg                                                                                          >  1:855895/1‑150 (MQ=255)
                                                  cccATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGAcc                                                                                       <  1:251346/150‑1 (MQ=255)
                                                    cATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCATACCGAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACcgc                                                                                     <  1:209068/150‑1 (MQ=255)
                                                     aTGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATACAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACcgcg                                                                                    >  1:8634/1‑150 (MQ=255)
                                                              ccAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATAGCCCTATATGCCTGGTAAGTCAAGCTGTGGCGGCGAGTGCTTCGCACCTGCACCCCCA‑ACAGCATGGACCGCGCGCTGCAcc                                                                           >  2:135232/1‑150 (MQ=14)
                                                                aaTGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAg                                                                         >  2:1067061/1‑150 (MQ=255)
                                                                    ccTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGtt                                                                     <  2:781680/150‑1 (MQ=255)
                                                                      tttGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACACCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTg                                                                   >  1:268991/1‑150 (MQ=255)
                                                                          ggCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTGGCGCTGAGGCGGCAGGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGCTTTGCGGc                                                               >  2:270920/1‑150 (MQ=255)
                                                                           gcgcATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGTCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCt                                                              >  2:920613/1‑150 (MQ=255)
                                                                           gcgcATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCt                                                              <  2:1142591/150‑1 (MQ=255)
                                                                                  aaCACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTg                                                       >  2:266582/1‑150 (MQ=255)
                                                                                  aaCACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTg                                                       >  2:862713/1‑150 (MQ=255)
                                                                                   acacCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGt                                                      >  1:651818/1‑150 (MQ=255)
                                                                                    cacCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGtt                                                     <  1:926783/150‑1 (MQ=255)
                                                                                       ccTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGAGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGc                                                  >  2:310249/1‑150 (MQ=255)
                                                                                        cTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGTACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCt                                                 <  1:342195/150‑1 (MQ=255)
                                                                                        cTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGACTGGTAAGACGCGCTGCGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCCGGGTTTGCGGCTGACTTTGTTCGCt                                                 >  2:590044/1‑150 (MQ=14)
                                                                                         tGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGAGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCtt                                                <  1:310249/150‑1 (MQ=255)
                                                                                              gAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGtt                                                     <  2:530144/140‑1 (MQ=255)
                                                                                              gAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGtt                                                     >  1:530144/1‑140 (MQ=255)
                                                                                               aaCAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGGCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGc                                          <  2:215396/150‑1 (MQ=255)
                                                                                                aCAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGc                                                               >  2:449283/1‑128 (MQ=255)
                                                                                                aCAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGc                                                               <  1:449283/128‑1 (MQ=255)
                                                                                                   gATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGa                                      >  2:712586/1‑150 (MQ=255)
                                                                                                    aTGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGaa                                     >  2:1101955/1‑150 (MQ=255)
                                                                                                        aGGCACAGCGGCGCACGGCCAGCCTCGGTCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACgg                                 <  2:251220/150‑1 (MQ=14)
                                                                                                        aGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACgg                                 >  2:574669/1‑150 (MQ=255)
                                                                                                        aGGCACAGAGCAGCACGGCCAGCATCGGCATGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACgg                                 <  1:57902/150‑1 (MQ=255)
                                                                                                            acaGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACGGCTGt                             <  1:1101955/150‑1 (MQ=255)
                                                                                                            acaGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACGGCTGt                             <  2:855895/150‑1 (MQ=255)
                                                                                                              agagCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTCCGCTTCGGTGCATGAAACGGCTGTCg                           >  1:543659/1‑150 (MQ=255)
                                                                                                                  cagcaCGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCt                                                 >  2:145483/1‑124 (MQ=255)
                                                                                                                  cagcaCGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCt                                                 <  1:145483/124‑1 (MQ=255)
                                                                                                                     caCGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCg                                                                 <  2:1113556/105‑1 (MQ=255)
                                                                                                                     caCGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCg                                                                 >  1:1113556/1‑105 (MQ=255)
                                                                                                                      aCGGCCAGCATCGGCCTGTATGTCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACGGCTGTCGATATTCGc                   <  1:920613/150‑1 (MQ=255)
                                                                                                                         gCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACGGCTGTCGATATTCGCCTg                >  1:110464/1‑150 (MQ=255)
                                                                                                                            aGCATCGTGCTCTATGCCTGGTACGTCGCGCTGGGCCGGCAAGTGCCGCGTCCCTTTTCATGCAGACGGC‑TGGACCGCCCCCTGTACCATTTTTTGCGGCTGCCTTTGTTCTCTTCGCTGCATGAAACGGCTGTCGATATTCGCCTGAGt             <  1:574669/150‑1 (MQ=2)
                                                                                                                                tCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACGGCTGTCGATATTCGCCTGAGTGAAg         <  1:1067061/150‑1 (MQ=255)
                                                                                                                                      tGTATGTCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACGGCTGTCGATATTCGCCTGAGTGAAGCTGCAc   <  2:71504/150‑1 (MQ=255)
                                                                                                                                      tGTATGTCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACGGCTGTCGATATTCGCCTGAGTGAAGCTGCAc   >  1:71504/1‑150 (MQ=255)
                                                                                                                                       gTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACTGCAGTGGCCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACGGCTGTCGATATTCGCCTGAGTGAAGCTGCACg  <  2:586501/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
ACCGCCGTCAGCGTCATCCTTCGCCCTGCCGGAAGCGCTTGCCAGGCAAGCCCATGCCGCTGCCAATGCCTTTGGGCGCATCAACACCCTGTACGAACAGATGCAGGCACAGAGCAGCACGGCCAGCATCGGCCTGTATGCCTGGTACGTCGCGCTGGGGCGGCAAGTGCTGCGCACCGCCACCCCCA‑ACAGCATGGACCGCGCGCTGCACCAGGGTTTGCGGCTGACTTTGTTCGCTTCGGTGCATGAAACGGCTGTCGATATTCGCCTGAGTGAAGCTGCACG  >  NC_002947/3837984‑3838268

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