Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,816,200 A→C 6.0% intergenic (+72/+186) katE → / ← chbG catalase HPII, heme d‑containing/chito‑oligosaccharide deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,816,2000AC 6.0% 780.2 / 32.6 234intergenic (+72/+186)katE/chbGcatalase HPII, heme d‑containing/chito‑oligosaccharide deacetylase
Reads supporting (aligned to +/- strand):  ref base A (110/108);  new base C (4/10);  total (116/118)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

Only 100 of 235 total aligned reads displayed.

GGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGA  >  NC_000913/1816053‑1816349
                                                                                                                                                   |                                                                                                                                                     
  aTGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATGGCGCCGCCGCAAATGATTTACAACGTGCGCTTGGTTTATGCCGGATGCGCCGGGAACGc                                                                                                                                                   >  1:3593364/1‑150 (MQ=255)
        tGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                             >  1:3116403/1‑150 (MQ=255)
        tGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                             >  1:1318853/1‑150 (MQ=255)
        tGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                             >  1:1957385/1‑150 (MQ=255)
           tAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                          >  1:1076987/1‑150 (MQ=255)
           tAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                          <  2:1938112/150‑1 (MQ=255)
           tAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                          <  2:2439264/150‑1 (MQ=255)
           gaaCGCTGAGGGCAGCACACCGCGTGTGGTCACGCATTCCTACGATTGACAAACTTCCGGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                          <  2:1649398/149‑1 (MQ=255)
            aaCGCTGATGGCAGCACCCCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                         <  2:1700798/150‑1 (MQ=255)
            aaCGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACCAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                         <  1:2172013/150‑1 (MQ=255)
             aCGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                                        <  2:675132/150‑2 (MQ=255)
             aCGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                                        <  2:1629622/150‑2 (MQ=255)
             aCGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                                        <  2:3667982/150‑2 (MQ=255)
             aCGCTGAGGGCAGCACCCCGCGTGTGGTCACGCATTCCTCAGCTTGACACAATTCCTGCCGGATGGGAGGCGGCGCCATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                                        <  2:233747/150‑2 (MQ=255)
              cGCTGATGGCAGCACACCGCGAGTGGTCACGCATTCCTAAGATTGACAAAATTCCTCCCTGAGGGGAGGCGGCGCAAATGCGCCGCCTCAAATGATTTACAAAGTGCGCTGTGTTTATGCCGGATGCGGCGTGAACGCCATATCCGGcca                                                                                                                                       >  2:3636905/1‑149 (MQ=38)
                cTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                     >  2:3098374/1‑150 (MQ=255)
                cTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                     >  2:532660/1‑150 (MQ=255)
                   aTGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaa                                                                                                                                  <  2:1381869/150‑1 (MQ=255)
                      gcagcaCACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACt                                                                                                                               <  1:1799108/150‑1 (MQ=255)
                       cagcaCCCCGCGTGTGGTCACGCATTCCTCAGAGTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAAGTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTGATCCGGCCTACAAAACtg                                                                                                                              <  2:3428697/150‑1 (MQ=38)
                        agcaCACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACtgt                                                                                                                             >  1:3530258/1‑150 (MQ=255)
                        agcaCACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACtgt                                                                                                                             <  2:1604950/150‑1 (MQ=255)
                        agcaCACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACtgt                                                                                                                             <  2:2969082/150‑1 (MQ=255)
                           acacCGCGTGTGGTCACGCATTCCTAAGATTGATAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCGCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCGTATCCGGCCTACAAAACTGTGcc                                                                                                                          >  1:125231/1‑149 (MQ=255)
                           acacCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCa                                                                                                                          >  2:2588067/1‑150 (MQ=255)
                                  gtgtGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCaa                                                                                                                   <  2:1129581/150‑1 (MQ=255)
                                    gtgGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAata                                                                                                                 <  1:477207/150‑1 (MQ=255)
                                        tCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCCAATTATTTACATAGTGCGCTTTGTTTATTCCGGCTTCGGCGTGTACGCCGTATCCGGCCTACAAAACGGTGCAAATTCAAGAAAtt                                                                                                             >  1:1042865/1‑150 (MQ=255)
                                        acACGCATTCCTAAGAATGACAAAATACCTGACTGATGGGAGGCCTCGCAATAGCGCCGCCTCAAATTATTTACATAGTGCGCTATGTTTATGCCGCATGCGGCGTTAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATAtt                                                                                                             >  2:3860125/2‑150 (MQ=25)
                                         cACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTg                                                                                                            <  2:3304441/150‑1 (MQ=255)
                                          aCGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGc                                                                                                           >  2:1209441/1‑150 (MQ=255)
                                          aCGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACAAAACTGTGCACATTCAATATATTGc                                                                                                           >  1:3806017/1‑150 (MQ=255)
                                           cGCAGTCCTAAGCTTGACAAAATTCCTGCCTGCTGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATGTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCa                                                                                                          <  2:3010129/150‑1 (MQ=255)
                                              aTTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCCCCGCAAATGATTTACATAGTGCGCTTTGTTTATTACGGCTGCGGCGTGAACGCCATATACGGCCTAAAAAACAGTGCAAATTAAAAAAATTGCAGGa                                                                                                       >  1:1475534/1‑150 (MQ=255)
                                              aTTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAAATGCGCCGCCTCAAATGATTTACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAAAGCCATATCCGGCCTACAAAACTGTGCAAATTAAAAATATTGCAGGa                                                                                                       >  1:2933221/1‑150 (MQ=37)
                                                   tAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAacac                                                                                                  <  1:2333034/150‑1 (MQ=255)
                                                   tAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAacac                                                                                                  <  1:492273/150‑1 (MQ=255)
                                                     aGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTTCGCCGCCGCAAAAGATTTAAATAGTGCGCTTTGTTTTTTGCGGATGCGGCGTGAACGCCGTATCCGGCCTACAAAAAAGTTCAAATTAAAAATATTGCAGGAAAAACGt                                                                                                >  1:1629857/1‑150 (MQ=255)
                                                      gATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTa                                                                                               >  2:452448/1‑150 (MQ=255)
                                                        ttGACAAAATTCCGGCCTGATGGGTGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATTTATTGCAGGAAACACGTAgg                                                                                             <  2:2989990/150‑1 (MQ=255)
                                                          gACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGcc                                                                                           >  2:1303970/1‑150 (MQ=255)
                                                          gACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGcc                                                                                           >  2:2074980/1‑150 (MQ=255)
                                                           aCAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCt                                                                                          <  2:843270/150‑1 (MQ=255)
                                                            cAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTg                                                                                         >  2:3435473/1‑150 (MQ=255)
                                                             acaaTTCCTGCCTGATGGGAGGCGGCGCAAGTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGa                                                                                        <  2:2690491/148‑1 (MQ=255)
                                                                aTTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCAACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATaa                                                                                     <  1:3417769/150‑1 (MQ=255)
                                                                            tGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATc                                                                         <  1:542357/150‑1 (MQ=255)
                                                                            tGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATc                                                                         <  2:1076987/150‑1 (MQ=255)
                                                                             gggAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCa                                                                        <  1:3881005/150‑1 (MQ=255)
                                                                             gggAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCa                                                                        <  2:697016/150‑1 (MQ=255)
                                                                               gAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGGTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAgg                                                                      >  2:3164194/1‑150 (MQ=255)
                                                                                  gcggcgCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAg                                                                   >  2:1779373/1‑150 (MQ=255)
                                                                                     gcgcAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGttt                                                                >  1:1415687/1‑150 (MQ=255)
                                                                                       gcAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTg                                                              >  2:280012/1‑150 (MQ=255)
                                                                                       gcAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTg                                                              <  1:3261684/150‑1 (MQ=255)
                                                                                        cAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGc                                                             <  1:31941/150‑1 (MQ=255)
                                                                                          aTTGCGCCGCCTCATATGATTTACATAGTGCGCTGTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGt                                                           <  1:2841896/150‑1 (MQ=255)
                                                                                          aTTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGt                                                           <  2:3384991/150‑1 (MQ=255)
                                                                                          aTTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGt                                                           <  2:2079016/150‑1 (MQ=255)
                                                                                             gcgcCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCCGTTTTGCGTTTg                                                        >  2:3814766/1‑150 (MQ=255)
                                                                                               gccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTc                                                      <  2:135829/150‑1 (MQ=255)
                                                                                                ccgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTca                                                     <  1:183944/150‑1 (MQ=255)
                                                                                                 cgccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCAACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcag                                                    >  1:1384504/1‑150 (MQ=255)
                                                                                                   ccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagca                                                  >  1:3409534/1‑150 (MQ=255)
                                                                                                   ccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagca                                                  >  1:882542/1‑150 (MQ=255)
                                                                                                   ccTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTcagca                                                  >  2:2484458/1‑150 (MQ=255)
                                                                                                        aaTGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGtctc                                             >  1:2667746/1‑150 (MQ=255)
                                                                                                        aaTGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGtctc                                             >  1:2924765/1‑150 (MQ=255)
                                                                                                        aaTGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGtctc                                             >  1:3608145/1‑150 (MQ=255)
                                                                                                        aaTGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGtctc                                             <  2:3350623/150‑1 (MQ=255)
                                                                                                         aTGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGAAGTCTCa                                            >  1:303847/1‑150 (MQ=255)
                                                                                                             tttACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAgcg                                        <  1:2465595/150‑1 (MQ=255)
                                                                                                              ttACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAgcgg                                       <  1:2105643/150‑1 (MQ=255)
                                                                                                                    aGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAGACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGt                                 >  2:3937664/1‑150 (MQ=255)
                                                                                                                    aGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGt                                 >  2:2531755/1‑150 (MQ=255)
                                                                                                                     gTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGtt                                <  1:1091241/150‑1 (MQ=255)
                                                                                                                        cgcTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTAcg                             <  2:3211189/150‑1 (MQ=255)
                                                                                                                          cTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATCTTGCAGGAAACACGTAGGCCTGATACGCGAAGCGCATCAGGCAGTATTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTAcgcc                           <  1:656266/150‑1 (MQ=255)
                                                                                                                          cTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTAcgcc                           >  1:797234/1‑150 (MQ=255)
                                                                                                                          cTTTGTTTAGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTAcgcc                           <  1:2698341/150‑1 (MQ=255)
                                                                                                                           tttGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTAcgccg                          >  2:2055755/1‑150 (MQ=255)
                                                                                                                                ttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAGGCGGCGTAACTGcgccgcttga                     <  1:1375745/150‑5 (MQ=18)
                                                                                                                                 tATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTg                    >  1:3748429/1‑150 (MQ=255)
                                                                                                                                   tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTa                  <  1:2125151/150‑1 (MQ=255)
                                                                                                                                   tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTa                  <  1:604174/150‑1 (MQ=255)
                                                                                                                                   tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTa                  <  2:3489238/150‑1 (MQ=255)
                                                                                                                                   tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATAGTGCCGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTa                  <  2:826689/150‑1 (MQ=255)
                                                                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAg                 >  2:1454498/1‑150 (MQ=21)
                                                                                                                                     ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAgg                >  2:3290308/1‑150 (MQ=255)
                                                                                                                                     ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAgg                >  2:3542080/1‑150 (MQ=21)
                                                                                                                                     ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTATTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAgg                >  2:1192422/1‑150 (MQ=25)
                                                                                                                                      cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGa               >  1:3141634/1‑150 (MQ=255)
                                                                                                                                      cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGa               >  2:2840559/1‑150 (MQ=255)
                                                                                                                                       ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGaa              <  1:2423381/150‑1 (MQ=255)
                                                                                                                                             ggcgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATc        <  1:3443975/150‑1 (MQ=255)
                                                                                                                                               cgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGACTTAATCGc      <  2:182957/150‑1 (MQ=255)
                                                                                                                                               cgTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGc      >  2:3573133/1‑150 (MQ=255)
                                                                                                                                                gTGAACGCCTTATCCGGCCTACCAAACTGTGCAAATTCAATATAGTGCAGGAAACACGGAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGcc     <  2:1372989/150‑1 (MQ=255)
                                                                                                                                                  gAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCgg   >  1:3689287/1‑150 (MQ=255)
                                                                                                                                                  gAACGCCTCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCgg   <  1:2268147/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGATGAACTGCTAACGCTGATGGCAGCACACCGCGTGTGGTCACGCATTCCTAAGATTGACAAAATTCCTGCCTGATGGGAGGCGGCGCAATTGCGCCGCCTCAAATGATTTACATAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCAAGCGGCGCAGTTACGCCGCCTTTGTAGGAATTAATCGCCGGA  >  NC_000913/1816053‑1816349

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