Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,787,672 C→G 5.7% intergenic (+76/‑49) 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,6720CG 5.7% 219.0 / 2.5 69intergenic (+76/‑49)ahpC/ahpFperoxiredoxin/alkylhydroperoxide reductase small subunit/alkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base C (29/36);  new base G (2/2);  total (31/38)
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.84e-01

AAGTGTGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCG  >  NC_002947/2787523‑2787813
                                                                                                                                                     |                                                                                                                                             
aaGTGTGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGAc                                                                                                                                               <  2:1424831/150‑1 (MQ=255)
      gCCCGGCCACGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGCCCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCCACCGCCAACGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGcc                                                                                                                                         >  1:757751/1‑150 (MQ=255)
          ggCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCggg                                                                                                                                     <  2:1088311/150‑1 (MQ=255)
             cAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGt                                                                                                                                  <  2:1547341/150‑1 (MQ=255)
               aGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttt                                                                                                                                <  2:13941/150‑1 (MQ=255)
               aGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttt                                                                                                                                <  2:2139171/150‑1 (MQ=255)
                gTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGtttt                                                                                                                               <  2:1281897/150‑1 (MQ=255)
                 tGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttttt                                                                                                                              <  2:711557/150‑1 (MQ=255)
                     aagaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTAt                                                                                                                          >  2:1603669/1‑150 (MQ=255)
                      agaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATc                                                                                                                         <  1:1990799/150‑1 (MQ=255)
                       gaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATcc                                                                                                                        <  2:1853595/150‑1 (MQ=255)
                         agGCGAAGCGGCTCTGGCGCCGTGGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCGGCCAAAGCTGAAGTACCTACCGCCCCGTCAGACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGa                                                                                                                      <  2:757751/150‑1 (MQ=255)
                         agGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGa                                                                                                                      <  2:370753/150‑1 (MQ=255)
                         agGCGAAGCCACTCTGGCGCCGTCGCTGGACCCTGCAGACCAGATCTAAGGTCAAGTTCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGTCTTCGCCTGGGCGTTTTTTTATCCGa                                                                                                                      <  1:967206/150‑1 (MQ=255)
                           gCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttttt                                                                                                                              >  1:1537385/1‑140 (MQ=255)
                           gCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttttt                                                                                                                              <  2:1537385/140‑1 (MQ=255)
                           gCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGt                                                                                                                    >  1:2338145/1‑150 (MQ=255)
                             gAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGttt                                                                                                                  <  1:1176017/150‑1 (MQ=255)
                              aaGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGtttt                                                                                                                 <  1:1505577/150‑1 (MQ=255)
                                  cACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaa                                                                                                             >  2:424054/1‑150 (MQ=255)
                                   accctGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACTCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCGGGCGt                                                                                                                                  <  2:210162/125‑1 (MQ=255)
                                   aCTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCGGGCGt                                                                                                                                  >  1:210162/1‑128 (MQ=255)
                                    ctctGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                           <  1:31749/150‑1 (MQ=255)
                                     tctGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                           <  1:1077141/149‑1 (MQ=255)
                                     tctGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                           >  2:1077141/1‑149 (MQ=255)
                                     tctGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAg                                                                                                          >  2:976470/1‑150 (MQ=255)
                                              gTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGccc                                                                                                 <  1:178429/150‑1 (MQ=255)
                                              gTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGccc                                                                                                 <  1:976470/150‑1 (MQ=255)
                                              gTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGccc                                                                                                 <  2:2338145/150‑1 (MQ=255)
                                               tCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGcc                                                                                                  <  2:210187/148‑1 (MQ=255)
                                               tCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGcc                                                                                                  >  1:210187/1‑148 (MQ=255)
                                               tCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCg                                                                                                >  1:1540552/1‑150 (MQ=255)
                                                  cTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTAt                                                                                             >  2:71008/1‑150 (MQ=255)
                                                   tGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                           <  2:2023844/135‑1 (MQ=255)
                                                   tGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                           >  1:2023844/1‑135 (MQ=255)
                                                        cTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGa                                                                                       >  1:387989/1‑150 (MQ=255)
                                                              ggCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCAc                                                                                 <  2:1962255/150‑1 (MQ=255)
                                                                cAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGc                                                                               >  1:1937305/1‑150 (MQ=255)
                                                                 aaGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGccc                                                                                                 >  2:628878/1‑131 (MQ=255)
                                                                 aaGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGccc                                                                                                 <  1:628878/131‑1 (MQ=255)
                                                                 aaGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCt                                                                              >  2:942911/1‑150 (MQ=255)
                                                                       tAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATc                                                                        >  1:725135/1‑150 (MQ=255)
                                                                           gTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCaa                                                                    >  1:136215/1‑150 (MQ=255)
                                                                              aTGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGa                                                                                       <  1:1037202/128‑1 (MQ=255)
                                                                              aTGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGa                                                                                       >  2:1037202/1‑128 (MQ=255)
                                                                                                gTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCgg                                               <  1:71008/150‑1 (MQ=255)
                                                                                                      agcCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGtt                                                                                          >  1:2295580/1‑101 (MQ=255)
                                                                                                      agcCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGtt                                                                                          <  2:2295580/101‑1 (MQ=255)
                                                                                                      agcCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCg                                                                       <  2:825576/120‑1 (MQ=255)
                                                                                                      agcCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCg                                                                       >  1:825576/1‑120 (MQ=255)
                                                                                                        ccgccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCa                                       >  1:771332/1‑150 (MQ=255)
                                                                                                         cgccAAAGCTGAAGTACCTACCGCCCCGTCAATCGCCCGGGCGACCTCGCCTGGGCGGTTTTGTATCCGAGTTTTGACCAAGGAATCGCCCATATGTTTGACTCCACTCGTAAATCGCAGCTGAAAACCTACTTGGAGCGGGTCACCCAg                                      <  1:636374/150‑1 (MQ=255)
                                                                                                           ccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCa                                           <  1:273172/143‑1 (MQ=255)
                                                                                                           ccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCa                                           >  2:273172/1‑143 (MQ=255)
                                                                                                           ccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGcc                                    <  2:136215/150‑1 (MQ=255)
                                                                                                             aaaGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTGACCCAGCCGa                                  >  1:1853438/1‑150 (MQ=255)
                                                                                                               aGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATc                                <  1:424054/150‑1 (MQ=255)
                                                                                                               aGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATc                                <  1:608030/150‑1 (MQ=255)
                                                                                                               aGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATc                                <  1:107103/150‑1 (MQ=255)
                                                                                                                gCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCg                               >  1:1515035/1‑150 (MQ=255)
                                                                                                                        acctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGtt                       <  1:1402085/150‑1 (MQ=255)
                                                                                                                            accGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCtt                   >  1:366433/1‑150 (MQ=255)
                                                                                                                                   ggTCCAACACCCGGGCAACCTCGCCTGGGCGTTTTTGTATCCGACTTTTGAAAAAGGAATCGCCCGTATGTTGGACCCCACGCTTAAATCGCAACGGAATACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcga            <  2:1007572/149‑1 (MQ=255)
                                                                                                                                    gTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcgac           <  2:387989/150‑1 (MQ=255)
                                                                                                                                         aCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACgg      >  2:106811/1‑150 (MQ=255)
                                                                                                                                           gCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGgcg    >  2:2184037/1‑150 (MQ=255)
                                                                                                                                             ccGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGcgcg  >  1:1392389/1‑150 (MQ=255)
                                                                                                                                             ccGGGCGACCTCGCCTGGGCATTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGAAAAAAATCATTGCTTCCCTCGACGACGGcgcg  >  1:1170679/1‑150 (MQ=255)
                                                                                                                                                 gCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACt                                                                  <  2:210684/82‑1 (MQ=255)
                                                                                                                                                 gCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACt                                                                  >  1:210684/1‑82 (MQ=255)
                                                                                                                                                     ccTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTg                      <  1:1157259/122‑1 (MQ=255)
                                                                                                                                                     ccTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTg                      >  2:1157259/1‑122 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
AAGTGTGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCG  >  NC_002947/2787523‑2787813

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: