Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 565,993 T→G 5.2% intergenic (+30/‑134) katA → / → bfr‑I catalase/bacterioferritin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947565,9930TG 5.2% 854.5 / 21.7 254intergenic (+30/‑134)katA/bfr‑Icatalase/bacterioferritin
Reads supporting (aligned to +/- strand):  ref base T (124/116);  new base G (7/7);  total (131/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

GCCGGTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCACCC  >  NC_002947/565844‑566140
                                                                                                                                                     |                                                                                                                                                   
  cGGTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTg                                                                                                                                                   <  1:5152379/150‑1 (MQ=255)
    gTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCGCATTTggg                                                                                                                                                 >  2:5386277/1‑150 (MQ=255)
        cATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGt                                                                                                                                             >  2:2852539/1‑150 (MQ=255)
          tGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCAAATGAGGGCGGttc                                                                                                                                           <  2:4311354/150‑2 (MQ=255)
             cgggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTc                                                                                                                                        <  1:529672/150‑1 (MQ=255)
             cgggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCAAATGAGGGCGGttcttc                                                                                                                                        <  2:4669083/150‑1 (MQ=255)
              gggcgTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCa                                                                                                                                       <  1:174562/150‑1 (MQ=255)
               ggcgTCAGCGAAGTTGTCATTCTGCGTCAGCTACTGTACTTCTTCTAGGCGGATCCGGCTTTCGGCGAAGGGTTTGCAAAGGCTTTGGGCATCCATCTCGCCTAAGTCGAAGTGTTAAGAAGATCCGGCCTGCTTTGGGCGGTTTTTCaa                                                                                                                                      <  2:3748465/150‑1 (MQ=255)
                    cAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAat                                                                                                                                 <  2:4828035/150‑1 (MQ=255)
                         aaGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCac                                                                                                                            <  1:3324067/150‑1 (MQ=255)
                          aGATGTCATTCATCGTCAGCTACAGTACTTCTTCTAGTCGGATCCGGCTTACGTCGAAGGGATTGCATAGGCTTTGGGCTTCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCact                                                                                                                           <  2:2711607/150‑1 (MQ=255)
                                 aTTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGttt                                                                                                                    >  1:5066247/1‑150 (MQ=255)
                                       cGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGa                                                                                                              <  1:2271455/150‑1 (MQ=255)
                                       cGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGa                                                                                                              >  2:698060/1‑150 (MQ=255)
                                          cAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGAttt                                                                                                           <  2:2116527/150‑1 (MQ=255)
                                            gCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGAttttt                                                                                                         <  1:886070/150‑1 (MQ=255)
                                             gttgaGTACTGCTTCTAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGGCTAAGTCGAAGTGATAAGAAGAAGGGCCCTCCTTTGGGCGGTTTTTGAATGAATATCACTACTGTTTAACCGAtttttt                                                                                                        <  1:3789152/146‑1 (MQ=255)
                                                cAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCg                                                                                                     >  1:3796071/1‑150 (MQ=255)
                                                  gTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTc                                                                                                   >  1:4699145/1‑150 (MQ=255)
                                                  gTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTc                                                                                                   >  1:4732193/1‑150 (MQ=255)
                                                     cttcttCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGa                                                                                                <  1:4709180/150‑1 (MQ=255)
                                                     cttcttCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGa                                                                                                <  2:4483013/150‑1 (MQ=255)
                                                     cttcttCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGa                                                                                                <  2:2827852/150‑1 (MQ=255)
                                                     cttcttCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGa                                                                                                <  2:4015003/150‑1 (MQ=255)
                                                      ttcttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGAt                                                                                               >  2:1655714/1‑150 (MQ=255)
                                                       actacAACGCGGATCCGGCTTACCGCGAAGGGATTGCAAAGGCATTGGGCATCAATCACGCCAAAGTCGAACTGATAAGAAGAACCGCCCTCATTTGGGCGGTATTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGAtt                                                                                              <  2:609437/146‑1 (MQ=255)
                                                        cttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCATAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTg                                                                                             <  2:41349/150‑1 (MQ=255)
                                                         ttcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTgc                                                                                            >  1:308933/1‑150 (MQ=255)
                                                          tcAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTgcg                                                                                           >  1:2633811/1‑150 (MQ=255)
                                                                 ggATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGttttt                                                                                    <  1:74202/150‑1 (MQ=255)
                                                                   aTCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCa                                                                                  <  2:3755792/150‑1 (MQ=255)
                                                                   aTCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCa                                                                                  <  2:1717553/150‑1 (MQ=255)
                                                                   aTCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCa                                                                                  <  1:209741/150‑1 (MQ=255)
                                                                    tCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAg                                                                                 >  1:3948451/1‑150 (MQ=255)
                                                                     ccGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGt                                                                                <  1:454431/150‑1 (MQ=255)
                                                                     ccGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGt                                                                                <  1:4495216/150‑1 (MQ=255)
                                                                       ggCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGa                                                                              >  2:4885866/1‑150 (MQ=255)
                                                                       ggCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGa                                                                              >  1:4428350/1‑150 (MQ=255)
                                                                           tACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCAt                                                                          <  1:4928625/150‑1 (MQ=255)
                                                                            aCGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATa                                                                         <  1:1284601/150‑1 (MQ=255)
                                                                            aCGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATa                                                                         <  1:4951893/150‑1 (MQ=255)
                                                                              ggCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGa                                                                       <  1:5494408/150‑1 (MQ=255)
                                                                                 gAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACaa                                                                    >  2:3138076/1‑150 (MQ=255)
                                                                                  aaGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACaaa                                                                   <  2:3371518/150‑1 (MQ=255)
                                                                                     ggATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGc                                                                >  2:3280992/1‑150 (MQ=255)
                                                                                     ggATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGc                                                                >  2:2646034/1‑150 (MQ=255)
                                                                                       aTTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTg                                                              <  1:4154091/150‑1 (MQ=255)
                                                                                       aTTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTg                                                              <  1:2035883/150‑1 (MQ=255)
                                                                                       aTTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTg                                                              <  1:1655714/150‑1 (MQ=255)
                                                                                        ttGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGTCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTgg                                                             <  2:4938911/150‑1 (MQ=255)
                                                                                        ttGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTgg                                                             <  1:775337/150‑1 (MQ=255)
                                                                                        ttGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTgg                                                             <  2:5504534/150‑1 (MQ=255)
                                                                                        ttGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTgg                                                             <  1:2741825/150‑1 (MQ=255)
                                                                                          gCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGtt                                                           >  2:1116995/1‑150 (MQ=255)
                                                                                          gCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGtt                                                           <  2:2449088/150‑1 (MQ=255)
                                                                                            aaaGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTCCTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTca                                                         >  2:5547878/1‑150 (MQ=255)
                                                                                            aaaGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTca                                                         >  2:4801821/1‑150 (MQ=255)
                                                                                            aaaGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTca                                                         >  2:4387690/1‑150 (MQ=255)
                                                                                             aaGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTcac                                                        <  1:1969583/150‑1 (MQ=255)
                                                                                             aaGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTcac                                                        >  1:4391202/1‑150 (MQ=255)
                                                                                               ggCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCacac                                                      <  1:1863610/150‑1 (MQ=255)
                                                                                                  tttGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTAt                                                   >  1:1741315/1‑150 (MQ=255)
                                                                                                  tttGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTAt                                                   >  2:1943258/1‑150 (MQ=255)
                                                                                                     gggCATCAATCTCGCCTAAGGCGAAGTGATAAGAAGAACCGCCCTCATGTGGGCGGTTTTTCAATGAATATCACTACTGATCAACCGATTTTGTGCGTCTGATTGCGCGTTTTTCAGTGACCAGAGACAAAAGCTGGTTCACACTAtgtg                                                <  1:631001/150‑1 (MQ=255)
                                                                                                      ggCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTgg                                               <  2:91718/150‑1 (MQ=255)
                                                                                                        cATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCa                                             <  1:2058047/150‑1 (MQ=255)
                                                                                                        cATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCa                                             <  1:349860/150‑1 (MQ=255)
                                                                                                        cATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCa                                             <  1:4020886/150‑1 (MQ=255)
                                                                                                         aTCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGa                                                                       >  1:1228273/1‑123 (MQ=255)
                                                                                                         aTCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGa                                                                       <  2:1228273/123‑1 (MQ=255)
                                                                                                         aTCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCAt                                            <  2:5241822/150‑1 (MQ=255)
                                                                                                            aaTCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGa                                         <  1:1517244/150‑1 (MQ=255)
                                                                                                             aTCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGAc                                        <  1:1115866/150‑1 (MQ=255)
                                                                                                             aTCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGAc                                        <  1:959575/150‑1 (MQ=255)
                                                                                                             aTCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGAc                                        <  1:4143828/150‑1 (MQ=255)
                                                                                                                tcGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGtt                                     >  2:355850/1‑150 (MQ=255)
                                                                                                                  gCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCTTTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGtttt                                   <  1:3754254/150‑1 (MQ=255)
                                                                                                                  gCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGtttt                                   >  1:2506510/1‑150 (MQ=255)
                                                                                                                    cTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTca                                 <  1:1888951/150‑1 (MQ=255)
                                                                                                                      aaGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTcaca                               <  1:2852539/150‑1 (MQ=255)
                                                                                                                          cGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGga                           <  2:3508818/150‑1 (MQ=255)
                                                                                                                           gAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGgag                          >  1:1459569/1‑150 (MQ=255)
                                                                                                                           gAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGgag                          >  2:2974685/1‑150 (MQ=255)
                                                                                                                              attcattgaaaaaCCGCCCAAATGAGGGCGGTTCTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCGTAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAAt                       >  1:2216328/8‑150 (MQ=38)
                                                                                                                               tGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCGCACTATGTGGCATAGACGTTTTCACATGGAGAATc                      >  2:416676/1‑150 (MQ=255)
                                                                                                                                 aTAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAg                    <  1:3384029/150‑1 (MQ=255)
                                                                                                                                   ttgaaaaaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAggg                  >  1:1927055/3‑150 (MQ=255)
                                                                                                                                   ttgaaaaaCCGCCCAAATGAGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTca                                                         >  2:3791095/3‑111 (MQ=38)
                                                                                                                                     gaagaaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTCTGTGGCATAGACGTTTTCACATGGAGAATCAGGGCg                >  1:3238382/1‑150 (MQ=255)
                                                                                                                                     gaagaaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCg                >  1:3666076/1‑150 (MQ=255)
                                                                                                                                     gaagaaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCg                >  2:2347719/1‑150 (MQ=255)
                                                                                                                                      aagaaCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGa               >  1:1037417/1‑150 (MQ=255)
                                                                                                                                       aaaaCCGCCCAAATGAGGGCGGTTCTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGAt              >  2:2272551/3‑150 (MQ=255)
                                                                                                                                          aCCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCa           <  2:3796071/150‑1 (MQ=255)
                                                                                                                                            cGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAg         <  2:2997961/150‑1 (MQ=255)
                                                                                                                                             gCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGACTCAGGGCGATGCAAgg        >  1:3856559/1‑150 (MQ=255)
                                                                                                                                              cccTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGc       >  2:1420979/1‑150 (MQ=255)
                                                                                                                                                cTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCa     >  1:3724743/1‑150 (MQ=255)
                                                                                                                                                   aTTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCAccc  >  1:2974309/1‑150 (MQ=255)
                                                                                                                                                    tttGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAAt                       <  2:1331979/128‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCGGTACCATGGCGGGCGTCAGCGAAGATGTCATTCAGCGTCAGCTACAGTACTTCTTCAAGGCGGATCCGGCTTACGGCGAAGGGATTGCAAAGGCTTTGGGCATCAATCTCGCCTAAGTCGAAGTGATAAGAAGAACCGCCCTCATTTGGGCGGTTTTTCAATGAATATCACTACTGTTTAACCGATTTTTTGCGTCTGATTGCGCGTTTTTCAGTGACCATAGACAAAAGCTGGTTCACACTATGTGGCATAGACGTTTTCACATGGAGAATCAGGGCGATGCAAGGCCACCC  >  NC_002947/565844‑566140

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