Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,563 G→T 17.2% 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,5630GT17.2% 412.2 / 72.3 163intergenic (‑221/+70)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (63/72);  new base T (13/15);  total (76/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

Only 100 of 172 total aligned reads displayed.

TTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATC  >  NC_000913/1579464‑1579662
                                                                                                   |                                                                                                   
ttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAg                                                                                                     <  2:2953843/100‑1 (MQ=255)
 tatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAgg                                                                                                    <  1:3715620/100‑1 (MQ=255)
 tatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAgg                                                                                                    <  2:2368511/100‑1 (MQ=255)
 tatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAgg                                                                                                    <  1:164031/100‑1 (MQ=255)
   tgCTTTTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGaa                                                                                                  <  1:414122/100‑1 (MQ=255)
   tgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGaa                                                                                                  >  2:2753480/1‑100 (MQ=255)
    gCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAt                                                                                                 >  2:2569605/1‑100 (MQ=255)
     cTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAtt                                                                                                >  1:127387/1‑100 (MQ=255)
        ttGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTAt                                                                                             <  1:1687280/100‑1 (MQ=255)
          gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGt                                                                                           >  2:1789357/1‑100 (MQ=255)
          gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGt                                                                                           >  2:776861/1‑100 (MQ=255)
              aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCGTGTCggg                                                                                       <  2:2307758/100‑1 (MQ=255)
              aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCggg                                                                                       <  1:385539/100‑1 (MQ=255)
              aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                        >  1:3798320/1‑99 (MQ=255)
               tcgaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGACGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGt                                                                                      >  1:3681208/4‑100 (MQ=255)
               tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                      <  2:1376089/100‑1 (MQ=255)
               tCAAGGATATTGCGAATGCCCTGATGCATACCACGTAAGCCAGGAGTTTGGTATGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGt                                                                                      >  1:3681142/1‑100 (MQ=255)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                     <  2:661691/100‑1 (MQ=38)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                     <  2:2150737/100‑1 (MQ=39)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                     <  1:2282100/100‑1 (MQ=38)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCggg                                                                                       >  1:2756659/1‑98 (MQ=255)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTg                                                                                     <  1:3974866/100‑1 (MQ=255)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTg                                                                                     <  1:2053851/100‑1 (MQ=255)
                 aaGGATATTTCTAATGCCGTGGTGCATACCACGTAAGCCAGGATTTGCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                     <  1:1125375/99‑1 (MQ=37)
                 aaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGTTTTTCGCAAGGGAAGATGCTGCTGGCGAAACACCCGACAGGAATTTATGTCGGGTGc                                                                                    <  1:3262206/100‑1 (MQ=255)
                 aaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                     <  2:2730188/99‑1 (MQ=38)
                  aGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgtt                                                                                   <  1:77545/100‑3 (MQ=38)
                   ggATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgttt                                                                                  <  2:1450716/100‑4 (MQ=38)
                   ggATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgttt                                                                                  <  1:547243/100‑4 (MQ=38)
                     atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgtttcg                                                                                <  1:3134287/100‑6 (MQ=37)
                     atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgtttcg                                                                                <  2:4225345/100‑6 (MQ=37)
                     atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTg                                                                                <  1:3634592/100‑1 (MQ=255)
                      tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTga                                                                               >  1:3998383/1‑99 (MQ=255)
                      tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGt                                                                               <  1:1652658/100‑1 (MQ=255)
                        ttGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTa                                                                             <  1:962576/100‑1 (MQ=255)
                          gCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAg                                                                           <  2:2834837/100‑1 (MQ=255)
                           cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                          <  2:4038217/100‑1 (MQ=255)
                            tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                         >  2:3379637/1‑100 (MQ=255)
                             aaTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                        >  2:4102686/1‑100 (MQ=255)
                              aTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCa                                                                       <  1:196061/100‑1 (MQ=255)
                                gCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATa                                                                     >  2:2227940/1‑100 (MQ=255)
                                 cccTTATACATACCACGTAAGCCAGGATTTTCCCAAGGGAAGACGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATa                                                                     <  1:3889535/99‑1 (MQ=255)
                                 cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                    <  2:85822/100‑1 (MQ=255)
                                 cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                    >  1:691808/1‑100 (MQ=255)
                                    tGATGCATACCTCGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGCATTTATGTCGGGTGCCTTGTTAAGGTCATaagaa                                                                 >  1:2318914/1‑100 (MQ=255)
                                       tGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagga                                                              >  2:407834/1‑100 (MQ=255)
                                       tGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagga                                                              <  2:1642371/100‑1 (MQ=255)
                                        gCATACGACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACAAGCGACAGGAATGGATGTCGGGTGCCTTGTTAAGGTCATAAGAagga                                                              <  1:3798712/99‑1 (MQ=255)
                                         cATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggagg                                                            <  2:2461519/100‑1 (MQ=255)
                                         cATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggagg                                                            <  2:464254/100‑1 (MQ=255)
                                         cATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggagg                                                            >  2:3455386/1‑100 (MQ=255)
                                         cATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggagg                                                            >  1:2410825/1‑100 (MQ=255)
                                          aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                           >  2:196860/1‑100 (MQ=255)
                                          aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                           >  2:3054431/1‑100 (MQ=255)
                                           tACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCt                                                          >  2:2529227/1‑100 (MQ=255)
                                            aCCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGATGCCTTGTTAAGGTCATAAGAAGGAGGCTa                                                         >  1:1278117/1‑100 (MQ=255)
                                            aCCACGTAAGCCAGGAGTTTCGCAAGGGAAGATGATGATGACGCAACACCCGACAGGAAGTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTa                                                         <  1:882798/100‑1 (MQ=255)
                                             ccACGTAAGCCAGGATATTCGCAGGGTGAGATGATGATGACGACACACCCGACAGGAAATTATGTCGGATGCCTTATTAAGGTCATAAGAAGGTGGCTa                                                         >  1:150408/1‑99 (MQ=255)
                                              cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAg                                                       <  2:691809/100‑1 (MQ=255)
                                              cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAg                                                       <  1:3684189/100‑1 (MQ=255)
                                                cGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGTCTTGTTAAGGTCATAAGAAGGAGGCTAAGaa                                                     >  2:3491501/1‑100 (MQ=255)
                                                cGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGaa                                                     >  2:2887654/1‑100 (MQ=255)
                                                  tAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATg                                                   <  2:837420/100‑1 (MQ=255)
                                                    aGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTACGAATGGa                                                 <  2:2318924/100‑1 (MQ=255)
                                                      ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGt                                               >  1:1627243/1‑100 (MQ=255)
                                                        aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                             <  1:2529214/100‑1 (MQ=255)
                                                         ggATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaa                                            >  2:2803945/1‑100 (MQ=255)
                                                          gATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaa                                           <  1:474308/100‑1 (MQ=255)
                                                           aTTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTAAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaaa                                          <  2:3113079/100‑1 (MQ=255)
                                                            ttttCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAag                                         >  2:4054816/1‑100 (MQ=255)
                                                                cGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagag                                     >  1:701264/1‑100 (MQ=255)
                                                                 gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGc                                    <  2:1434344/100‑1 (MQ=255)
                                                                 gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGc                                    <  2:639868/100‑1 (MQ=255)
                                                                 gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGc                                    >  1:1752299/1‑100 (MQ=255)
                                                                  cAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCg                                   >  2:531022/1‑100 (MQ=255)
                                                                     gggAAGATGATGATGACGAAACACCCGACAGGAATTTGTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGtta                                >  1:2356284/1‑100 (MQ=255)
                                                                     gggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGtta                                <  2:3624862/100‑1 (MQ=255)
                                                                         aGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaa                            >  2:2408043/1‑100 (MQ=255)
                                                                            tgatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGTAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatta                         >  2:1186679/1‑100 (MQ=255)
                                                                            tgatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatt                          >  1:1818551/1‑99 (MQ=255)
                                                                             gatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattat                        <  1:1189222/100‑1 (MQ=255)
                                                                                cttaacaaggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtct                     >  2:4007132/11‑100 (MQ=25)
                                                                                cttaacaaggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtct                     >  2:3421365/11‑100 (MQ=25)
                                                                                  tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTcc                   >  2:3218916/1‑100 (MQ=255)
                                                                                  tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTcc                   >  1:3034876/1‑100 (MQ=255)
                                                                                    aCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAt                 <  2:719540/100‑1 (MQ=255)
                                                                                       aggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTg              >  1:972731/4‑100 (MQ=38)
                                                                                        ggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGt             >  2:1789875/3‑100 (MQ=38)
                                                                                          cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt           >  2:3738161/1‑100 (MQ=38)
                                                                                          cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt           >  1:4282195/1‑100 (MQ=38)
                                                                                          cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt           >  1:337706/1‑100 (MQ=38)
                                                                                           acCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAGGAGAGCGTTATTAATTATTCTCCATTTGttt           <  1:3396891/99‑1 (MQ=255)
                                                                                           acCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt          >  1:1203660/1‑100 (MQ=255)
                                                                                           acCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt            <  1:2951510/98‑1 (MQ=255)
                                                                                             ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGc        >  2:949836/1‑100 (MQ=255)
                                                                                             ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGc        <  2:2066676/100‑1 (MQ=255)
                                                                                              cGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCa       >  1:3384795/1‑100 (MQ=255)
                                                                                              cGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCa       <  2:2281580/100‑1 (MQ=255)
                                                                                                   ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGATTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc  >  2:516111/1‑100 (MQ=255)
                                                                                                   ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc  >  2:425553/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATC  >  NC_000913/1579464‑1579662

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