Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,218,221 T→C 7.8% intergenic (‑47/+193) PP_1958 ← / ← PP_5495 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,218,2210TC 7.8% 142.7 / 2.0 51intergenic (‑47/+193)PP_1958/PP_5495hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/22);  new base C (2/2);  total (27/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

GTGAAGGGATCAGAGGAGCTATGGCATTCCCTCTAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGCCCATGG  >  NC_002947/2218081‑2218364
                                                                                                                                             |                                                                                                                                               
gtgAAGGGATCAGAGGAGCTATGGCATTCCCTCTAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGctc                                                                                                                                         >  6:294985/1‑149 (MQ=255)
  gAAGGGATCAGAGGAGCTATGGCATTCCCTCTAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTg                                                                                                                                       >  2:42547/1‑149 (MQ=255)
   aaGGGATCAGAGGAGCTATGGCATTCCCTCTAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTg                                                                                                                                       <  6:189441/148‑1 (MQ=255)
           agagGAGCTATGGCATTCCCTCTAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCggg                                                                                                                               >  4:276280/1‑148 (MQ=255)
              ggagCTATGGCATTCCCTCTAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATg                                                                                                                           <  5:294985/149‑1 (MQ=255)
              ggagCTATGGCATTCCCTCTAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATg                                                                                                                           <  5:178139/149‑1 (MQ=255)
               gagCTATGGCATTCCCTCTAAGACCATCCCCCGGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATg                                                                                                                           >  1:47545/1‑149 (MQ=255)
                         aTTCCCTCTAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATCGa                                                                                                                                               <  4:158175/118‑1 (MQ=255)
                         aTTCCCTCTAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATCGa                                                                                                                                               >  3:158175/1‑118 (MQ=255)
                              ctctAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGgcgc                                                                                                                    >  7:64439/1‑140 (MQ=255)
                              ctctAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGt                                                                                                           >  8:5294/1‑149 (MQ=255)
                                ctAAGACCATCCCCCGTGTACGCGGCCTCAACCAGCCACACACTCTGCCTCGCTGGCCTCATTCCGCTCTGCGGCTATTTCGATGTGTGTGTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGCAGTTc                                                                                                         >  5:76074/1‑149 (MQ=255)
                                actaGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTc                                                                                                         <  5:362699/146‑1 (MQ=255)
                                 tAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATCGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCg                                                                                                        >  2:171594/1‑149 (MQ=255)
                                   aGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTc                                                                                                         <  6:163071/146‑1 (MQ=255)
                                         tCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAgg                                                                                                >  6:278767/1‑149 (MQ=255)
                                            cccGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAagag                                                                                                 <  6:241371/145‑1 (MQ=255)
                                            cccGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAagag                                                                                                 >  5:241371/1‑145 (MQ=255)
                                                 tgtACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGt                                                                                         <  4:137400/149‑1 (MQ=255)
                                                      gcgTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGt                                                                                     >  2:125821/1‑148 (MQ=255)
                                                            tCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTc                                                                                                         >  2:326668/1‑122 (MQ=255)
                                                            tCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTc                                                                                                         <  1:326668/122‑1 (MQ=255)
                                                               cccAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAAtt                                                                           <  7:5294/149‑1 (MQ=255)
                                                                ccAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCTTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTAc                                                                                       >  6:273181/1‑136 (MQ=255)
                                                                ccAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTAc                                                                                       <  5:273181/136‑1 (MQ=255)
                                                                ccAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGAc                                                                                  <  1:258928/141‑1 (MQ=255)
                                                                ccAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGAc                                                                                  >  2:258928/1‑141 (MQ=255)
                                                                   gCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCTGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTa                                                                       <  1:271529/149‑1 (MQ=255)
                                                                          ccTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCg                                                                >  3:6833/1‑149 (MQ=255)
                                                                           ctctGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTAtt                                                                                                                                                 >  8:71039/1‑67 (MQ=255)
                                                                           ctctGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGTTACCGGTACGTGACAGTAATTTCTACGACTCgg                                                               <  1:42547/149‑1 (MQ=255)
                                                                                ccTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATCGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGtata                                                          <  1:171594/149‑1 (MQ=255)
                                                                                 cTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAg                                                         <  3:276280/149‑1 (MQ=255)
                                                                                  tCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAgg                                                        >  6:181330/1‑149 (MQ=255)
                                                                                        gCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAaga                                                                                                  <  4:252870/101‑1 (MQ=255)
                                                                                        gCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAaga                                                                                                  >  3:252870/1‑101 (MQ=255)
                                                                                          cTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTAc                                                  >  4:128483/1‑147 (MQ=255)
                                                                                          cTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTAc                                                  >  3:47087/1‑147 (MQ=255)
                                                                                           tCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAg                                               >  6:81092/1‑149 (MQ=255)
                                                                                           tCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACAGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAg                                               >  7:211522/1‑149 (MQ=255)
                                                                                             aTGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAg                                               <  5:181330/147‑1 (MQ=255)
                                                                                                 cGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAg                                         >  6:135808/1‑149 (MQ=255)
                                                                                                           cTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCt                               >  3:319059/1‑149 (MQ=255)
                                                                                                            tATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTg                              >  3:147859/1‑149 (MQ=255)
                                                                                                                tCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTAc                                                  >  1:200940/1‑125 (MQ=255)
                                                                                                                tCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTAc                                                  <  2:200940/125‑1 (MQ=255)
                                                                                                                          ttttACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGc                 <  5:81092/148‑1 (MQ=255)
                                                                                                                            ttACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTg               <  4:6833/148‑1 (MQ=255)
                                                                                                                                gAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAg            >  6:327601/1‑147 (MQ=255)
                                                                                                                                  aTCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGc        <  2:47545/149‑1 (MQ=255)
                                                                                                                                         gTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGCCCATgg  <  4:147859/148‑1 (MQ=255)
                                                                                                                                         gTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGCCCATgg  <  3:128483/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GTGAAGGGATCAGAGGAGCTATGGCATTCCCTCTAAGACCATCCCCCGTGTACGCGTCCTCACCCAGCCACACCCTCTGCCTCGCTGGCCTCATGCCGCTCTGCGGCTATTTCGATGTGTGTTTTACGAATCGTTAGTATTGAGCGCTCTGGTAGCGGGGATGGAGGCGCGGAAGAAGTTCGCCAAGAGGTACCGGTACGTGACAGTAATTTCTACGACTCGGGTATAGGTCTTACGAGCATCAGGTTGCTTTCTGCCCCTCGCAATGCTGGAGTAGCCCATGG  >  NC_002947/2218081‑2218364

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