Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,563 G→T 6.4% intergenic (‑221/+70) yddA ← / ← ydeM putative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,579,5630GT 6.4% 332.0 / 9.8 117intergenic (‑221/+70)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (47/62);  new base T (5/3);  total (52/65)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 118 total aligned reads displayed.

TAACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT  >  NC_000913/1579418‑1579712
                                                                                                                                                 |                                                                                                                                                     
tAACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCAAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAt                                                                                                                                                   >  1:408469/1‑150 (MQ=255)
  aCAATATGAAATATTAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATAACACTTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAATCACCCGACTGGAAttt                                                                                                                                                 >  2:1589042/1‑150 (MQ=255)
  aCAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAttt                                                                                                                                                 >  1:691330/1‑150 (MQ=255)
  aCAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGAGGACGAAACACCCGACAGGAAttt                                                                                                                                                 >  2:1595960/1‑150 (MQ=255)
   cAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTa                                                                                                                                                <  1:465237/150‑1 (MQ=255)
     atatGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                                                                              <  2:1803001/150‑1 (MQ=255)
        tGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCg                                                                                                                                           >  1:1470860/1‑150 (MQ=255)
         taaaTATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                                                                          <  1:200312/149‑1 (MQ=255)
         gAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                                                                          >  2:665285/1‑150 (MQ=255)
           aaTATAAAGTAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGt                                                                                                                                        <  1:1360957/150‑1 (MQ=255)
            atatAAAGGATTATATCATTTTCTGTGGTCATGGTTATGATATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCTCAAGGTAAGATGATGATGACGAAACACCCGACATAAATTCCTGTcttgtg                                                                                                                                       <  1:87556/150‑6 (MQ=255)
             tataAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGc                                                                                                                                      <  1:1506100/150‑1 (MQ=255)
             tataAAGGAATATATCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGCTGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGc                                                                                                                                      <  1:761205/150‑1 (MQ=255)
              ataAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGcc                                                                                                                                     <  1:780193/150‑1 (MQ=255)
                aaaGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCtt                                                                                                                                   >  2:1529470/1‑150 (MQ=255)
                 aaGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCGAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTg                                                                                                                                  >  2:520510/1‑150 (MQ=255)
                 aaGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGAGGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTg                                                                                                                                  <  1:1668530/150‑1 (MQ=255)
                   ggAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGtt                                                                                                                                <  2:122784/150‑1 (MQ=255)
                    gAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTa                                                                                                                               <  1:486921/150‑1 (MQ=255)
                    gAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTa                                                                                                                               >  2:683601/1‑150 (MQ=255)
                      atatTTCATTTTCAGTGGTTATCGTTATGCTATTGATAATCATTGATATTGCTAATGCCCTGATTCTTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATTATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAg                                                                                                                             <  1:205268/150‑1 (MQ=255)
                       tatTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                                                                            <  1:1232211/150‑1 (MQ=255)
                        atTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                                                                           <  2:1380209/150‑1 (MQ=255)
                        atTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                                                                           <  2:385509/150‑1 (MQ=255)
                         tttCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                                                                          >  1:982010/1‑150 (MQ=255)
                         tttCATTTTCAGTGGTTATGGTTATGCAATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATATTCGCAAGGTAAGATGATGAAGACGAAACACCCGACAGGTATTTATGTCGCGTGCCTTTATAATGTc                                                                                                                          >  2:1137688/1‑150 (MQ=255)
                            cATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATa                                                                                                                       >  1:1711817/1‑150 (MQ=255)
                            cATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGAGCCTTGTTAAGTTCATa                                                                                                                       >  2:130886/1‑150 (MQ=255)
                            cATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCAAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATa                                                                                                                       <  1:578369/150‑1 (MQ=255)
                              ttttCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTCAAGGTCATaag                                                                                                                     <  1:566117/150‑1 (MQ=255)
                               tttCAGTGGTTATGGTTATGCTATTGATAATCAATGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTTAGGTCATaaga                                                                                                                    <  2:1389829/150‑1 (MQ=255)
                                    gTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggag                                                                                                               <  1:1396106/150‑1 (MQ=255)
                                     tggttatggttaTGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggagg                                                                                                              >  1:61213/1‑150 (MQ=255)
                                      ggttatggttatGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                                                                             <  1:1263747/150‑1 (MQ=255)
                                      ggttatggttatGCTATTGATAATCAAGGATATTGCTAATGACCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                                                                             <  2:66290/150‑1 (MQ=255)
                                       gttatggttatgTTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCATGGGAAGATGTTGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCt                                                                                                            <  1:355064/150‑1 (MQ=255)
                                             gttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAAt                                                                                                      >  2:873060/1‑150 (MQ=255)
                                              ttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATg                                                                                                     <  1:1157779/150‑1 (MQ=255)
                                                 tgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                                                                  <  2:1679489/150‑1 (MQ=255)
                                                 tgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                                                                  <  2:1431866/150‑1 (MQ=255)
                                                   cTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGtt                                                                                                >  2:1732203/1‑150 (MQ=255)
                                                    tATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                                                                               <  2:1741265/150‑1 (MQ=255)
                                                    tATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                                                                               <  2:369973/150‑1 (MQ=255)
                                                       tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaaa                                                                                            <  2:286214/150‑1 (MQ=255)
                                                        gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAag                                                                                           <  1:1170150/150‑1 (MQ=255)
                                                            aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagag                                                                                       >  2:459975/1‑150 (MQ=255)
                                                              cAAGTATATTGCTAATTCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGAGGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCg                                                                                     <  2:436468/150‑1 (MQ=255)
                                                                 ggATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGGGCCTTGTTAAGGTCATAAGACGGAGGCTAAGAAGGTAGTTAAAATAGAGCGttt                                                                                  >  2:723000/1‑149 (MQ=255)
                                                                  gATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGAATTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAAATATGTCGGGTGCCATGTTAAGGTCATAAGAAGGAGGCAAAGAATGCAGTTAAAAGAGAGCGttat                                                                                 >  2:613726/1‑150 (MQ=255)
                                                                    tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatta                                                                               >  1:991471/1‑150 (MQ=255)
                                                                      ttGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaat                                                                             <  1:1688561/150‑1 (MQ=255)
                                                                      ttGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaat                                                                             <  1:771704/150‑1 (MQ=255)
                                                                      ttGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaat                                                                             <  1:501760/150‑1 (MQ=255)
                                                                       tGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatt                                                                            >  2:1249940/1‑150 (MQ=255)
                                                                       tGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatt                                                                            >  1:543590/1‑150 (MQ=255)
                                                                           aaTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattatta                                                                        >  2:1492477/1‑149 (MQ=255)
                                                                           aaTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtc                                                                        >  2:1044238/1‑150 (MQ=255)
                                                                            aTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtct                                                                       >  1:1525248/1‑150 (MQ=255)
                                                                                 cTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTAATTAATTATTCTCCAtt                                                                  <  2:237875/150‑1 (MQ=255)
                                                                                        aTACCACGTAATCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTg                                                           <  2:1370322/150‑1 (MQ=255)
                                                                                        aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTg                                                           <  2:1665655/150‑1 (MQ=255)
                                                                                        aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGAAAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAATGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTg                                                           <  2:1601187/150‑1 (MQ=255)
                                                                                              cGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACAt                                                     <  1:1529470/150‑1 (MQ=255)
                                                                                                tAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATcc                                                   <  2:1080482/150‑1 (MQ=255)
                                                                                                tAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAACGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATcc                                                   <  1:1587681/150‑1 (MQ=255)
                                                                                                       ggATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAAtgt                                            >  1:585988/1‑150 (MQ=255)
                                                                                                          ttttCGCAAGTGAAGATGATGATGACGAAACACCCGATAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGt                                         <  2:1063929/150‑1 (MQ=255)
                                                                                                            ttCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGATCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTaa                                       <  2:1545000/150‑1 (MQ=255)
                                                                                                                cAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTAATGGCGGGTGCCTTGTTAAGGTCAAAAGAAGGAGGCTAAGAATGGAGATAAAAGACAGCGTAATTAATTATTATCCATTTTTTTTGAAACATCCATAATTTGGAAAAGCg                                   >  1:158436/1‑150 (MQ=255)
                                                                                                                 aaGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCgg                                  <  2:1726793/150‑1 (MQ=255)
                                                                                                                 aaGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCgg                                  <  1:665285/150‑1 (MQ=255)
                                                                                                                 aaGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCgg                                  <  2:926493/150‑1 (MQ=255)
                                                                                                                  aGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGt                                 <  2:519569/150‑1 (MQ=255)
                                                                                                                   gggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTa                                <  2:902572/150‑1 (MQ=255)
                                                                                                                    ggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCCTTATTAAATATTATACATTTGTATTGAAAAATCCATAATGTTGTAAAGCGGTaa                               >  2:271174/1‑149 (MQ=255)
                                                                                                                      aaGATTATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCATTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTActc                             <  1:1595960/150‑1 (MQ=255)
                                                                                                                      aaGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTActc                             <  2:691330/150‑1 (MQ=255)
                                                                                                                      aaGATGATGATGAAGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCAAAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATAATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTActc                             <  1:1589042/150‑1 (MQ=255)
                                                                                                                        tatgatgatgACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTg                           <  2:408469/149‑1 (MQ=255)
                                                                                                                         atgatgatgACGAAGCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGt                          >  2:1362099/1‑150 (MQ=255)
                                                                                                                         atgatgatgACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGt                          <  1:520510/150‑1 (MQ=255)
                                                                                                                           tatgattaCTAAACTCCCTACATGAATTTATGTCTGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTc                        <  1:312704/149‑1 (MQ=255)
                                                                                                                            ttgttgaCGAAACACCCGACAGGAATTTATGTCGGGTGACTTGTTAAGGTCATAAGAAGGAGGCTAAGATTGGAGTTAAAAGAGAGCGTTTTTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCt                       <  2:176615/146‑1 (MQ=255)
                                                                                                                            atgttgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCt                       >  1:1723683/1‑150 (MQ=255)
                                                                                                                             tgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCtt                      <  2:982010/150‑1 (MQ=255)
                                                                                                                               atgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATCTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAg                    >  1:1769460/1‑150 (MQ=255)
                                                                                                                                tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGc                   >  2:1241747/1‑150 (MQ=255)
                                                                                                                                 gaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCt                  >  2:453191/1‑150 (MQ=255)
                                                                                                                                   caaggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTaa                >  1:1669123/6‑150 (MQ=255)
                                                                                                                                   cGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTaa                >  2:878743/1‑150 (MQ=255)
                                                                                                                                    gAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAAc               >  1:477438/1‑150 (MQ=255)
                                                                                                                                       gcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTAAGCTAActct            >  1:888015/2‑150 (MQ=255)
                                                                                                                                        cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTa           >  2:1454071/1‑150 (MQ=255)
                                                                                                                                        cacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTa           >  1:331257/1‑150 (MQ=255)
                                                                                                                                        cacCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAAATCTa           >  2:1682482/1‑150 (MQ=255)
                                                                                                                                          cccGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTATTTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTAcc         >  1:590559/1‑150 (MQ=255)
                                                                                                                                           ccGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCa        >  2:1322315/1‑150 (MQ=255)
                                                                                                                                            cGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCAt       >  1:1575926/1‑150 (MQ=255)
                                                                                                                                             gACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCAtt      >  1:1139544/1‑149 (MQ=255)
                                                                                                                                                 ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGt  <  1:1003532/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TAACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT  >  NC_000913/1579418‑1579712

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