Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,787,673 C→G 5.2% intergenic (+77/‑48) ahpC → / → ahpF peroxiredoxin/alkylhydroperoxide reductase small subunit/alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,787,6730CG 5.2% 268.4 / 2.4 86intergenic (+77/‑48)ahpC/ahpFperoxiredoxin/alkylhydroperoxide reductase small subunit/alkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base C (44/37);  new base G (3/2);  total (47/39)
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.99e-01

CGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGAAGTCCCGC  >  NC_002947/2787532‑2787822
                                                                                                                                             |                                                                                                                                                     
cGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTgg                                                                                                                                               >  1:191636/1‑150 (MQ=255)
cGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTgg                                                                                                                                               <  1:581055/150‑1 (MQ=255)
      aGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttt                                                                                                                                         <  2:1809462/150‑1 (MQ=255)
        tGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttttt                                                                                                                                       <  2:1302458/150‑1 (MQ=255)
        tGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttttt                                                                                                                                       >  2:1039533/1‑150 (MQ=255)
        tGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCCGGCGACCTCGCCTGGGCGttttt                                                                                                                                       >  2:611010/1‑150 (MQ=255)
             agaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCTGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATc                                                                                                                                  <  1:764740/150‑1 (MQ=255)
             agaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATc                                                                                                                                  <  1:1339743/150‑1 (MQ=255)
              gaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATcc                                                                                                                                 <  1:1493081/150‑1 (MQ=255)
                    gAAGCCACTCTGGCGCCGTCGCAGGACCTCGTAGGCAAGAAATAAGGTCATGCCCATAACGCGGTCGGTAATCAACCGCCAAAGCTGAAGTACCTACCGCCCCCTCAAACTCCCGCGCGACCTCGCCTAGTCGCTTTTGTATCCTAGctt                                                                                                                           >  2:2253197/1‑150 (MQ=255)
                     aaGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGtttt                                                                                                                          <  2:191636/150‑1 (MQ=255)
                       gCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGGCGACATCGGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGa                                                                                                                        <  1:1621621/150‑1 (MQ=255)
                        ccACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaa                                                                                                                       >  1:491726/1‑150 (MQ=255)
                          aCTCTGGCGCCGTCGCTCGGGCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaa                                                                                                                     >  2:2394402/1‑150 (MQ=255)
                            tctGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGa                                                                                                                        >  2:354996/1‑145 (MQ=255)
                            tctGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGa                                                                                                                        <  1:354996/145‑1 (MQ=255)
                            tctGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAg                                                                                                                   >  1:1461119/1‑150 (MQ=255)
                            cctGGCGCCGTCGACAGAACTCGTAGGCAAGATCTAAGGTCATGGCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCTCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTCTATCCGAGTTTTGaaaacg                                                                                                                   <  2:140723/149‑3 (MQ=255)
                                      tCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCg                                                                                                         >  2:1365215/1‑150 (MQ=255)
                                        gCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTa                                                                                                       >  2:544613/1‑150 (MQ=255)
                                         cTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttt                                                                                                                                         <  1:1323724/115‑1 (MQ=255)
                                         cTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttt                                                                                                                                         >  2:1323724/1‑115 (MQ=255)
                                         cTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTAt                                                                                                      <  2:1461119/150‑1 (MQ=255)
                                          tGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATg                                                                                                     <  1:1039533/150‑1 (MQ=255)
                                                           aTCTAAGGTCATGCGCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTgg                                                                                                                                               >  2:1299288/1‑91 (MQ=255)
                                                           aTCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTgg                                                                                                                                               <  1:1299288/91‑1 (MQ=255)
                                                           aTCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTaa                                                                                    <  1:1365215/150‑1 (MQ=255)
                                                           aTCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTAAAACGCCCGGGCGCCGGCGCCTGGCCGCTTTGTTAAGGGAGTTTTGAAAAAGGAATCGCGCGTATGTTGGACGTCACGcctta                                                                                    >  2:1477786/1‑146 (MQ=255)
                                                             cTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAAt                                                                                  >  1:2215471/1‑150 (MQ=255)
                                                              tAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCGGGc                                                                                                                                             >  1:1860176/1‑90 (MQ=38)
                                                              tAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCGGGc                                                                                                                                             <  2:1860176/90‑1 (MQ=38)
                                                                 ggTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCa                                                                              >  1:1410129/1‑150 (MQ=255)
                                                                  gTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGtttt                                                                                                                          <  1:2177270/105‑1 (MQ=255)
                                                                  gTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGtttt                                                                                                                          >  2:2177270/1‑105 (MQ=255)
                                                                        cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttttt                                                                                                                                       <  2:1611716/86‑1 (MQ=255)
                                                                        cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttttt                                                                                                                                       >  1:1611716/1‑86 (MQ=255)
                                                                        cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCaa                                                                             <  2:1312764/144‑1 (MQ=255)
                                                                        cccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCaa                                                                             >  1:1312764/1‑144 (MQ=255)
                                                                         ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                                    >  2:1250597/1‑104 (MQ=255)
                                                                         ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGaaaa                                                                      >  1:1632054/1‑150 (MQ=255)
                                                                         ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGaaaa                                                                      >  2:1088707/1‑150 (MQ=255)
                                                                         ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGaaaa                                                                      >  2:1693398/1‑150 (MQ=255)
                                                                         ccATAACGCGGGCGGTAAGCAGACGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                                    <  1:1250597/104‑1 (MQ=255)
                                                                           aTAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAAcc                                                                    >  1:688683/1‑150 (MQ=255)
                                                                           aTAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAAcc                                                                    >  1:1059764/1‑150 (MQ=255)
                                                                                gcgggcggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACtt                                                               <  2:1410129/150‑1 (MQ=255)
                                                                                 cgggcggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTg                                                              <  2:2215471/150‑1 (MQ=255)
                                                                                    gcggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTa                                                                                                       <  2:992789/106‑1 (MQ=255)
                                                                                    gcggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTa                                                                                                       >  1:992789/1‑106 (MQ=255)
                                                                                     cggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGc                                                          <  1:184279/150‑1 (MQ=255)
                                                                                       gTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCgg                                                        <  2:2367285/150‑1 (MQ=255)
                                                                                           gcagcCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCa                                                    >  1:516371/1‑150 (MQ=255)
                                                                                              gccgccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGCCGCACCACGAACAACGCAACGCAATACCTACGGCGAGCGGGTCAgcc                                                 >  2:1832240/1‑150 (MQ=255)
                                                                                              gccgccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGGTTTTGTATCCGAGTTTTGAAAAACGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCAccc                                                 >  2:1711724/1‑150 (MQ=255)
                                                                                                    aaaGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAAc                                                                            >  2:685743/1‑117 (MQ=255)
                                                                                                    aaaGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAAc                                                                            <  1:685743/117‑1 (MQ=255)
                                                                                                     aaGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGAt                                          <  2:1263581/150‑1 (MQ=255)
                                                                                                      aGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATg                                                                                                     <  1:147884/90‑1 (MQ=255)
                                                                                                      aGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATg                                                                                                     >  2:147884/1‑90 (MQ=255)
                                                                                                       gCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGtttt                                                                                                                          <  1:1754342/68‑1 (MQ=255)
                                                                                                       gCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGtttt                                                                                                                          >  2:1754342/1‑68 (MQ=255)
                                                                                                              tacctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACtt                                                               <  2:917355/120‑1 (MQ=255)
                                                                                                              tacctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACtt                                                               >  1:917355/1‑120 (MQ=255)
                                                                                                               acctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGtt                                >  2:1333764/1‑150 (MQ=255)
                                                                                                               acctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAAGAAGCGCTTGGATTGTGGTCGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGtt                                <  2:196350/150‑1 (MQ=255)
                                                                                                                cctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCgaga                                     <  1:1281590/144‑1 (MQ=255)
                                                                                                                cctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCgaga                                     >  2:1281590/1‑144 (MQ=255)
                                                                                                                 ctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGtt                                                                                                   <  2:1572824/81‑1 (MQ=255)
                                                                                                                 ctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGtt                                                                                                   >  1:1572824/1‑81 (MQ=255)
                                                                                                                 ctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGc                              >  2:2236588/1‑150 (MQ=255)
                                                                                                                  taccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCt                             >  1:170843/1‑150 (MQ=255)
                                                                                                                    ccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGt                                                                                                    >  2:1380628/1‑77 (MQ=255)
                                                                                                                    ccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGt                                                                                                    <  1:1380628/77‑1 (MQ=255)
                                                                                                                       ccccGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGc                                                          >  2:585112/1‑116 (MQ=255)
                                                                                                                       ccccGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGc                                                          <  1:585112/116‑1 (MQ=255)
                                                                                                                       ccccGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCACGTATGTTGGACTCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCt                        >  2:1891352/1‑150 (MQ=255)
                                                                                                                           atCAGACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGCGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcgac                    >  1:1667374/2‑150 (MQ=255)
                                                                                                                            tCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTGACATCCAAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcgaca                   >  1:1998287/1‑149 (MQ=255)
                                                                                                                              aaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATc                                                                                 >  2:2097954/1‑86 (MQ=37)
                                                                                                                              aaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATc                                                                                 <  1:2097954/86‑1 (MQ=37)
                                                                                                                               aaCGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGacgacg                >  1:1689403/1‑150 (MQ=255)
                                                                                                                                   cccGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGc                                              <  1:2398677/116‑1 (MQ=255)
                                                                                                                                   cccGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGc                                              >  2:2398677/1‑116 (MQ=255)
                                                                                                                                      gggCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAACACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGaa         <  2:491726/150‑1 (MQ=255)
                                                                                                                                      gggCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGaa         >  1:1554543/1‑150 (MQ=255)
                                                                                                                                           aCCTCGCCTGGGCGTTTTTGTATCCGAGTTATCAAAAATGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGAAGTccc    <  2:1998287/150‑1 (MQ=255)
                                                                                                                                            ccTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGAAGTCCcg   <  2:688683/150‑1 (MQ=255)
                                                                                                                                             cTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGAGTCACCCAGCCGATCTAGAACTTTGCTTCCCTCGACGACGGCGCGAAGTCCcga  >  2:1330709/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
CGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGAAGTCCCGC  >  NC_002947/2787532‑2787822

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: