Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,180,811 G→A 7.1% W211* (TGG→TAG)  PP_3678 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,180,8110GA 7.1% 156.4 / 2.4 56W211* (TGG→TAG) PP_3678hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base A (2/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

GAATTAAAAATAAATCCTCATAATGGGCGACAGTATCTAGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTAT  >  NC_002947/4180666‑4180951
                                                                                                                                                 |                                                                                                                                            
gAATTAAAAATAAATCCTCATAATGGGCGACAGTATCTAGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGaa                                                                                                                                          <  1:1047355/150‑1 (MQ=255)
       aaataaatCCTCATAATGGGCGACAGTATCTAGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGctcc                                                                                                                                   <  1:381559/150‑1 (MQ=255)
               ccTCATAATGGGCGACAGTATCTAGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGgc                                                                                                                           >  2:1599033/1‑150 (MQ=255)
                 tCATAATGGGCGACAGTATCTAGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGgcgc                                                                                                                         <  1:1231140/150‑1 (MQ=255)
                 tCATAATGGGCGACAGTATCTAGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATACAACAAATGAACGTCGTGATGTTAAAGAGTGACGATTGACAAGTGCTCCCTCCAGgcgc                                                                                                                         >  2:1413089/1‑150 (MQ=255)
                    tAATGGGCGACAGTATCTAGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAAc                                                                                                                      <  2:1637264/150‑1 (MQ=255)
                     aaTGGGCGACAGTATCTAGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAg                                                                                                                             <  1:1042332/142‑1 (MQ=255)
                     aaTGGGCGACAGTATCTAGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAg                                                                                                                             >  2:1042332/1‑142 (MQ=255)
                          gCGACAGTATCTAGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCaa                                                                                                                <  2:491187/150‑1 (MQ=255)
                            gACAGTATCTAGACCATCCCGACCTTGTGCCTCCTGAAGATGAAGTATATGACAAACCGCCTATTGATGAGCATACCATCCAACAACTGCAAGATGACATATTTAGAAGGCACGACTGTGAAAAGCTCCCTCCAGGCGCAACCCTCAAGt                                                                                                              >  1:2686/1‑150 (MQ=255)
                              cAGTATCTAGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCATGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTAc                                                                                                            >  1:1361739/1‑150 (MQ=255)
                                      aGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACaa                                                                                                    >  2:292330/1‑150 (MQ=255)
                                             cccGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGAt                                                                                             >  2:1194887/1‑150 (MQ=255)
                                                 aCCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCtat                                                                                         <  1:1091473/150‑1 (MQ=255)
                                                           tACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCcac                                                                               <  2:103729/150‑1 (MQ=255)
                                                                      aaGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTatat                                                                                       <  1:979260/131‑1 (MQ=255)
                                                                      aaGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTatat                                                                                       >  2:979260/1‑131 (MQ=255)
                                                                      aaGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTtata                                                                    >  2:19300/1‑150 (MQ=255)
                                                                              gACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGaa                                                            >  2:1215223/1‑150 (MQ=255)
                                                                               aCAAAGCGCCTATTGATGAGCATCCCAACCAAATGCCTTAGTATGGAATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGaaa                                                           <  1:1599033/150‑1 (MQ=255)
                                                                               aCAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGaaa                                                           <  1:1194887/150‑1 (MQ=255)
                                                                                 aaaGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATg                                                              <  2:570966/145‑1 (MQ=255)
                                                                                 aaaGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATg                                                              >  1:570966/1‑145 (MQ=255)
                                                                                 aaaGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATAAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGa                                                         <  1:1413089/150‑1 (MQ=255)
                                                                                           ttGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTa                                               <  1:292330/150‑1 (MQ=255)
                                                                                              atgaGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTatca                                            >  2:1311436/1‑150 (MQ=255)
                                                                                              atgaGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTatca                                            >  2:1018622/1‑150 (MQ=255)
                                                                                              atgaGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTACCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTatca                                            >  1:186848/1‑150 (MQ=255)
                                                                                                  gCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTCTATGTCGCCACAAATCTATATAGTAATGAAATACGCATGGTGATCATCAg                                        >  2:1250211/1‑150 (MQ=255)
                                                                                                     tACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAACGACGCATGTTTATCATCAGtat                                     >  2:325333/1‑150 (MQ=255)
                                                                                                     tACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGgat                                     >  1:1100225/1‑150 (MQ=255)
                                                                                                      aCCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGCAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATTTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTat                                              <  1:268773/140‑1 (MQ=255)
                                                                                                      aCCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTat                                              >  2:268773/1‑140 (MQ=255)
                                                                                                            cAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGaaa                              <  1:1311436/150‑1 (MQ=255)
                                                                                                            cAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGaaa                              >  2:386786/1‑150 (MQ=255)
                                                                                                                  cTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATg                        <  1:1018622/150‑1 (MQ=255)
                                                                                                                   tGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATAttt                                                                 >  2:64867/1‑108 (MQ=255)
                                                                                                                   tGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATgg                       <  1:386786/150‑1 (MQ=255)
                                                                                                                   tGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTACCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATgg                       <  2:186848/150‑1 (MQ=255)
                                                                                                                   tGCAAGATGAGATATTTAGAAGGCACGACCGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATAttt                                                                 <  1:64867/108‑1 (MQ=255)
                                                                                                                     cAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGc                     >  2:713063/1‑150 (MQ=255)
                                                                                                                        gATGAGATATTTAGAAGGCACGACTAGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATGTGTCGCCACAAATCttt                                                                       <  1:349050/97‑1 (MQ=255)
                                                                                                                        gATGAGATATTTAGAAGGCACGACTAGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATGTGTCGCCACAAATCttt                                                                       >  2:349050/1‑97 (MQ=255)
                                                                                                                        gATGAGATATTTAGAAGGCACGACTAGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCttt                                                                       <  1:349047/97‑1 (MQ=255)
                                                                                                                        gATGAGATATTTAGAAGGCACGACTAGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCttt                                                                       >  2:349047/1‑97 (MQ=255)
                                                                                                                         aTGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGc                 <  1:1215223/150‑1 (MQ=255)
                                                                                                                          tGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCa                    >  2:302939/1‑146 (MQ=255)
                                                                                                                          tGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCa                    <  1:302939/146‑1 (MQ=255)
                                                                                                                              atatTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCa            <  1:713063/150‑1 (MQ=255)
                                                                                                                              accactAGAAAGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCa            >  2:1426776/6‑150 (MQ=255)
                                                                                                                                atTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAACTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCaga          <  2:1100225/150‑1 (MQ=255)
                                                                                                                                     gAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAAc     <  1:19300/150‑1 (MQ=255)
                                                                                                                                     gAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCATCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAAc     <  1:1250211/150‑1 (MQ=255)
                                                                                                                                        ggCACGACTGTGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACtat  <  1:1295482/150‑1 (MQ=255)
                                                                                                                                             gACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGc                 >  2:1486251/1‑130 (MQ=255)
                                                                                                                                             gACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGc                 <  1:1486251/130‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
GAATTAAAAATAAATCCTCATAATGGGCGACAGTATCTAGACCATCCCGACCTTGTGCCTACTGAAGATGAAGTATATGACAAAGCGCCTATTGATGAGCATACCATCCAAAAACTGCAAGATGAGATATTTAGAAGGCACGACTGGGAAGAGCTCCCTCCAGGCGCAACCCTCAAGTACAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTAT  >  NC_002947/4180666‑4180951

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