Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,563 G→T 15.5% 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,5630GT15.5% 309.7 / 48.6 125intergenic (‑221/+70)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (51/53);  new base T (5/14);  total (58/67)
Fisher's exact test for biased strand distribution p-value = 8.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-02

Only 100 of 132 total aligned reads displayed.

TTATGCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATC  >  NC_000913/1579464‑1579662
                                                                                                   |                                                                                                   
ttatgCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAg                                                                                                     <  2:2257098/100‑1 (MQ=255)
    gCTATTGATAATCAAGGATATTGCTAATGCCCTGCTGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGaa                                                                                                  >  2:337544/1‑99 (MQ=255)
    gCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAt                                                                                                 >  1:2074533/1‑100 (MQ=255)
     cTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAtt                                                                                                >  2:181128/1‑100 (MQ=255)
      tATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAttt                                                                                               <  1:2142715/100‑1 (MQ=255)
      tATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAAtt                                                                                                <  2:2271533/99‑1 (MQ=255)
       aTTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTa                                                                                              <  2:320158/100‑1 (MQ=255)
        ttGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTAt                                                                                             <  2:755064/100‑1 (MQ=255)
         tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                            >  2:3291830/1‑100 (MQ=255)
             aaTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                        <  2:1778827/100‑1 (MQ=255)
              aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCggg                                                                                       >  1:2601729/1‑100 (MQ=255)
               tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGTTGTCGAAACACCCGACATAAATTCCTGTCGGGt                                                                                      <  2:2388897/100‑1 (MQ=255)
               tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                      <  1:325205/100‑1 (MQ=255)
               tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                      <  1:3236333/100‑1 (MQ=255)
               tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                      <  1:240659/100‑1 (MQ=255)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                     <  2:1150217/100‑1 (MQ=38)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                     <  1:1550950/100‑1 (MQ=39)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                     <  1:3062146/100‑1 (MQ=38)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                     <  2:3306543/100‑1 (MQ=38)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                     <  1:854063/100‑1 (MQ=38)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTg                                                                                     <  1:1184602/100‑1 (MQ=255)
                cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTg                                                                                     <  1:3306518/100‑1 (MQ=255)
                  aGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGcc                                                                                   <  1:647456/100‑1 (MQ=255)
                     atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgtttcg                                                                                <  1:1193356/100‑6 (MQ=37)
                     atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTg                                                                                <  1:1994755/100‑1 (MQ=255)
                        ttGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgtttcgtca                                                                             <  1:2815858/100‑9 (MQ=25)
                         tGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGGAACACGCGACAGGACTTTATGTCGGGTGCCTTGTTaa                                                                            >  2:2460990/1‑100 (MQ=255)
                         tGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTaa                                                                            >  2:2001382/1‑100 (MQ=255)
                           cTAATGCCCTTATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGTCAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                          <  2:1519499/100‑1 (MQ=255)
                           cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                          <  2:3252236/100‑1 (MQ=255)
                            tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                         <  2:1048950/100‑1 (MQ=255)
                            tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                         <  1:592632/100‑1 (MQ=255)
                             aaTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                        <  2:3308516/100‑1 (MQ=255)
                               tGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCAt                                                                      >  2:1867896/1‑100 (MQ=255)
                                gCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATa                                                                     <  1:1604219/100‑1 (MQ=255)
                                 cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                    >  2:2753622/1‑100 (MQ=255)
                                   cTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaaga                                                                  >  2:3182584/1‑100 (MQ=255)
                                       tGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagga                                                              <  1:3490809/100‑1 (MQ=255)
                                       tGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAagga                                                              >  1:1538957/1‑100 (MQ=255)
                                       ngCATATCACGTAAGCCAGGCTTTTCGCACGGGAAGATGATGATGACTAAACACCCGACAGGGATTTATGAAGGGTGCCTTGTTCAAGtaataa                                                                    >  2:3103432/2‑94 (MQ=255)
                                        gCATACCACGTAAGCCAGGATTTTCGCAAGTGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggag                                                             <  2:442775/100‑1 (MQ=255)
                                        gCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggag                                                             >  2:306479/1‑100 (MQ=255)
                                        ccgatctaCGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggag                                                             <  1:1849922/93‑1 (MQ=255)
                                          aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                           <  1:2895233/100‑1 (MQ=255)
                                          aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                           <  1:219357/100‑1 (MQ=255)
                                           tACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGACTTTTTAAGGTCATAAGAAGAAGGCt                                                          >  2:2527411/1‑100 (MQ=255)
                                            aCCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTa                                                         <  2:3598148/100‑1 (MQ=255)
                                             ccACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTaa                                                        <  2:2089632/100‑1 (MQ=255)
                                             ccACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTaa                                                        >  1:3130544/1‑100 (MQ=255)
                                              cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAg                                                       <  1:688486/100‑1 (MQ=255)
                                               aCGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGa                                                      >  1:2399531/1‑100 (MQ=255)
                                                    aGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGa                                                 <  2:738584/100‑1 (MQ=255)
                                                    aGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGACGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGa                                                 >  2:1564949/1‑100 (MQ=255)
                                                     gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                >  2:360074/1‑100 (MQ=255)
                                                       cAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGtt                                              <  1:897315/100‑1 (MQ=255)
                                                       cAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTGATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                >  2:1368185/1‑98 (MQ=255)
                                                         ggATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaa                                            >  1:1887480/1‑100 (MQ=255)
                                                         ggATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTCTTAAGGTCATAAGAAGGAGGATAAGAATGACGTTaa                                            >  2:2954225/1‑100 (MQ=255)
                                                            ttttCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAag                                         <  2:1444543/100‑1 (MQ=255)
                                                            ttttCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAag                                         <  2:1887471/100‑1 (MQ=255)
                                                            ttttCGCAAGGAAAGATGATGATGACGAAACACCAGACAGGAATTTATGTCGGGTGCATTGTTAAGGTCATAAGAAGGAGGTTAAGAATAGAGTTAAAag                                         >  2:1471208/1‑100 (MQ=255)
                                                              ttCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagag                                       >  1:3350469/1‑100 (MQ=255)
                                                                cGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagagag                                     <  1:1680913/100‑1 (MQ=255)
                                                                    aGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGtt                                 >  2:3482849/1‑100 (MQ=255)
                                                                    aGGGAAGATGATGAGGACGAAACACCCGACAGGAAGTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGt                                  <  2:2981475/98‑1 (MQ=255)
                                                                      ggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttat                               >  2:181166/1‑100 (MQ=255)
                                                                      ggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttat                               <  1:2727740/100‑1 (MQ=255)
                                                                      ggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttat                               <  2:2055711/100‑1 (MQ=255)
                                                                        aaGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatta                             >  1:1029917/1‑100 (MQ=255)
                                                                            tgatgattaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatta                         >  1:1906022/1‑100 (MQ=255)
                                                                            tgatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatta                         >  1:2603473/1‑100 (MQ=255)
                                                                             gatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattat                        <  2:826933/100‑1 (MQ=255)
                                                                              atgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGCTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattat                        <  2:344261/99‑1 (MQ=255)
                                                                               tgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtc                      >  1:1435171/1‑100 (MQ=255)
                                                                               tgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtc                      >  1:3559987/1‑100 (MQ=255)
                                                                                  tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTcc                   <  1:1115746/100‑1 (MQ=255)
                                                                                  tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTcc                   >  2:422152/1‑100 (MQ=255)
                                                                                   gaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                  >  2:1861276/1‑100 (MQ=255)
                                                                                   gaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                  <  1:1424069/100‑1 (MQ=255)
                                                                                   gaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCa                  <  1:781215/100‑1 (MQ=255)
                                                                                    aCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAt                 <  1:801608/100‑1 (MQ=255)
                                                                                      gAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAttt               <  1:306481/100‑1 (MQ=255)
                                                                                      aaggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAttt               >  2:2749768/5‑100 (MQ=37)
                                                                                       aaaCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTg              >  1:986162/1‑100 (MQ=255)
                                                                                       aaaCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTg              <  1:2975577/100‑1 (MQ=255)
                                                                                          cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt           >  1:2703015/1‑100 (MQ=39)
                                                                                          cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt           >  1:17683/1‑100 (MQ=38)
                                                                                           acCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt          >  2:2313261/1‑100 (MQ=255)
                                                                                            cccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTg         <  1:1861285/100‑1 (MQ=255)
                                                                                            cccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTg         >  2:433687/1‑100 (MQ=255)
                                                                                            cccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTg         >  1:884853/1‑100 (MQ=255)
                                                                                              cGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCa       <  2:3504707/100‑1 (MQ=255)
                                                                                              cGACAAGAACTTATGTCGGGTGCACTGTTAAGGTCATAAGACGGACGCTAAAACTAGAGCTAAGAGCGCGCGATATTAATCACTCTCCCTTTGtttctca       >  2:1617315/1‑96 (MQ=255)
                                                                                               gACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCaa      >  2:1257482/1‑100 (MQ=255)
                                                                                               gACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCaa      >  2:1057157/1‑100 (MQ=255)
                                                                                                 cAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACa    <  1:1993435/100‑1 (MQ=255)
                                                                                                   ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACAt   >  2:2815528/1‑99 (MQ=255)
                                                                                                   ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc  <  1:2403992/100‑1 (MQ=255)
                                                                                                   ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc  >  1:1735003/1‑100 (MQ=255)
                                                                                                   ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc  <  1:1221482/100‑1 (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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: