Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,984,227 G→A 6.4% intergenic (+34/+125) aslB → / ← aslA putative anaerobic sulfatase maturation enzyme AslB/arylsulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,984,2270GA 6.4% 493.1 / 14.5 155intergenic (+34/+125)aslB/aslAputative anaerobic sulfatase maturation enzyme AslB/arylsulfatase
Reads supporting (aligned to +/- strand):  ref base G (80/65);  new base A (6/4);  total (86/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 162 total aligned reads displayed.

ATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGC  >  NC_000913/3984078‑3984372
                                                                                                                                                     |                                                                                                                                                 
aTCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATg                                                                                                                                                   <  2:709392/150‑1 (MQ=255)
 tCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgc                                                                                                                                                  >  1:58530/1‑150 (MQ=255)
  cAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgcg                                                                                                                                                 >  2:1880256/1‑150 (MQ=255)
   aGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgcgg                                                                                                                                                <  2:2043748/150‑1 (MQ=255)
   aGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCGGCGGTTGATTGCCGGATgcgg                                                                                                                                                <  1:1872591/150‑1 (MQ=255)
      gTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGt                                                                                                                                             <  1:938620/150‑1 (MQ=255)
      gTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGt                                                                                                                                             <  1:862324/150‑1 (MQ=255)
          tttCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAAc                                                                                                                                         <  2:1784683/150‑1 (MQ=255)
          tttCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAAc                                                                                                                                         >  1:2433827/1‑150 (MQ=255)
           ttCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACg                                                                                                                                        >  1:369740/1‑150 (MQ=255)
             ccGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGcc                                                                                                                                      >  2:2267543/1‑150 (MQ=255)
              cGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCt                                                                                                                                     <  2:917206/150‑1 (MQ=255)
              cGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCt                                                                                                                                     >  1:2350333/1‑150 (MQ=255)
              cGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATAAGGCGTTTACGccg                                                                                                                                     <  1:1122177/150‑2 (MQ=255)
                      aCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCTGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                             >  2:2252621/1‑150 (MQ=255)
                      aCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                             <  2:797079/150‑1 (MQ=255)
                      aCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                             <  2:1157985/150‑1 (MQ=255)
                       ccGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTa                                                                                                                                   >  1:114040/1‑143 (MQ=255)
                        cGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                           >  1:1522721/1‑150 (MQ=255)
                          ccATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                         <  2:2099337/150‑1 (MQ=255)
                           cATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                                                                        <  2:676428/150‑1 (MQ=255)
                            atatCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAt                                                                                                                       <  1:42481/150‑1 (MQ=255)
                            atatCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCGGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAt                                                                                                                       <  2:338066/150‑1 (MQ=255)
                                 ttAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCg                                                                                                                  <  1:345533/150‑1 (MQ=255)
                                 ttAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCg                                                                                                                  <  2:1387606/150‑1 (MQ=255)
                                   aaaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTg                                                                                                                >  1:208780/1‑150 (MQ=255)
                                       cAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCaaa                                                                                                            <  1:108028/150‑1 (MQ=255)
                                        aaTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAAt                                                                                                           >  2:2239603/1‑150 (MQ=255)
                                         aTGTCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAAtt                                                                                                          <  1:1793733/150‑1 (MQ=255)
                                              tGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATa                                                                                                     >  1:1388810/1‑150 (MQ=255)
                                               gATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATaa                                                                                                    >  2:654889/1‑150 (MQ=255)
                                                  ttGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAAtt                                                                                                 >  1:546259/1‑150 (MQ=255)
                                                   tGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTg                                                                                                >  1:2195351/1‑150 (MQ=255)
                                                    gCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGc                                                                                               <  2:144696/150‑1 (MQ=255)
                                                      tGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAg                                                                                             <  2:646793/150‑1 (MQ=255)
                                                       ggCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGc                                                                                            >  2:833093/1‑150 (MQ=255)
                                                       ggCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGc                                                                                            >  2:76464/1‑150 (MQ=255)
                                                       ggCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGc                                                                                            >  2:108843/1‑150 (MQ=255)
                                                       ggCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGc                                                                                            <  2:608949/150‑1 (MQ=255)
                                                                 cGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGTGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAgg                                                                                  >  1:595209/1‑150 (MQ=255)
                                                                  gCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGc                                                                                 >  2:1912235/1‑150 (MQ=255)
                                                                    ccGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTg                                                                               <  1:700642/150‑1 (MQ=255)
                                                                     cGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAACGTTCTGTAGGCTgg                                                                              <  2:1142373/150‑1 (MQ=255)
                                                                            cGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCTGCCGCCTGCGGTTTATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGAt                                                                       <  2:1589733/150‑1 (MQ=255)
                                                                                  tATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCa                                                                 <  2:208780/150‑1 (MQ=255)
                                                                                   atgcatgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCCGgatgcgatgctgacg                                                                <  2:48592/150‑6 (MQ=255)
                                                                                    tgcatgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgc                                                               >  1:2482279/1‑150 (MQ=255)
                                                                                    tgcatgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgc                                                               <  2:1835764/150‑1 (MQ=255)
                                                                                     gcatgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgca                                                              <  2:1619942/150‑1 (MQ=255)
                                                                                     gcatgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCGGCGGGTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgca                                                              <  1:1367325/150‑1 (MQ=255)
                                                                                       atgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatc                                                            <  2:506319/150‑1 (MQ=255)
                                                                                        tgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcg                                                           >  1:1375243/1‑150 (MQ=255)
                                                                                        tgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcg                                                           >  2:1213066/1‑150 (MQ=255)
                                                                                         gcgcATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcgc                                                          <  1:662480/150‑1 (MQ=255)
                                                                                            cATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcgcatc                                                       >  1:2225056/1‑150 (MQ=255)
                                                                                             aTTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATcc                                                      >  1:1891474/1‑150 (MQ=255)
                                                                                               ttGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCgg                                                    <  1:743857/150‑1 (MQ=255)
                                                                                                    ggtggtGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAg                                               >  1:1212189/1‑150 (MQ=255)
                                                                                                    ggtggtGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAg                                               <  2:1522721/150‑1 (MQ=255)
                                                                                                      tggtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGc                                             >  1:2398745/1‑150 (MQ=255)
                                                                                                       ggtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCa                                            <  1:2019774/150‑1 (MQ=255)
                                                                                                        gtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAg                                           >  2:354395/1‑150 (MQ=255)
                                                                                                        gtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAg                                           <  1:443898/150‑1 (MQ=255)
                                                                                                         tgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGa                                          >  2:200778/1‑150 (MQ=255)
                                                                                                         tgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGa                                          >  2:64080/1‑150 (MQ=255)
                                                                                                         tgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGa                                          >  1:2326083/1‑150 (MQ=255)
                                                                                                            gtaagtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGAtct                                       >  2:1174952/1‑150 (MQ=255)
                                                                                                             taagtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGAtctc                                      >  2:419637/1‑150 (MQ=255)
                                                                                                               agtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGAGGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAg                                    >  2:231534/1‑150 (MQ=255)
                                                                                                                gtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGc                                   <  1:2030225/150‑1 (MQ=255)
                                                                                                                gtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGGTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGc                                   <  1:554581/150‑1 (MQ=255)
                                                                                                                 taGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCg                                  <  2:1356209/150‑1 (MQ=255)
                                                                                                                   gAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGAt                                <  1:711866/150‑1 (MQ=255)
                                                                                                                    aaaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATa                               >  2:2056809/1‑150 (MQ=255)
                                                                                                                     aaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAg                              <  2:595209/150‑1 (MQ=255)
                                                                                                                     aaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAg                              <  1:2375352/150‑1 (MQ=255)
                                                                                                                     aaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAg                              <  1:2280024/150‑1 (MQ=255)
                                                                                                                     aaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAg                              >  1:1789721/1‑150 (MQ=255)
                                                                                                                      aTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAgc                             >  1:1975180/1‑150 (MQ=255)
                                                                                                                         ggcggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGcc                          >  2:2093321/1‑150 (MQ=255)
                                                                                                                         ggcggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGcc                          >  2:1738710/1‑150 (MQ=255)
                                                                                                                          gcggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCg                         <  1:2006190/150‑1 (MQ=255)
                                                                                                                           cggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCgg                        >  2:1782812/1‑150 (MQ=255)
                                                                                                                            ggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGc                       <  2:704539/150‑1 (MQ=255)
                                                                                                                            ggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTGCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGc                       <  2:712919/150‑1 (MQ=255)
                                                                                                                             gccgccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCt                      >  1:599698/1‑150 (MQ=255)
                                                                                                                               cgccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCGGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTa                    >  2:1396129/1‑150 (MQ=255)
                                                                                                                                gccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAg                   >  1:1160432/1‑150 (MQ=255)
                                                                                                                                    gCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGGCAg               >  2:2474087/1‑150 (MQ=255)
                                                                                                                                        ttGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGAATCGCATCCGGCAAAGGAAGATATCAGCGATAGCGCCGGCTTAGTCAGAttt           >  2:2101965/1‑150 (MQ=255)
                                                                                                                                          gATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGGAGCGTTCTGTAGGCTGGATAAGATGAGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTaa         >  1:508584/1‑150 (MQ=255)
                                                                                                                                          gATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTaa         <  1:1942889/150‑1 (MQ=255)
                                                                                                                                            ttGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATc       >  1:1977570/1‑150 (MQ=255)
                                                                                                                                              gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTg     <  2:1109058/150‑1 (MQ=255)
                                                                                                                                              gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTg     >  1:1661633/1‑150 (MQ=255)
                                                                                                                                               ccGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgc    >  2:1208362/1‑150 (MQ=255)
                                                                                                                                               ccGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgc    >  2:1701802/1‑150 (MQ=255)
                                                                                                                                               ccGGATAAGGCGTTTACGCCGCATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgc    >  1:2366302/1‑150 (MQ=39)
                                                                                                                                               ccGGATAAGGCGTTTACGCCGCATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgc    >  2:2327910/1‑150 (MQ=39)
                                                                                                                                               ccGGATAAGGCGTTTACGCCGCATCCGGCCTACATGATCGTGCAAATTAAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgc    >  2:726981/1‑150 (MQ=38)
                                                                                                                                                     |                                                                                                                                                 
ATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGC  >  NC_000913/3984078‑3984372

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