Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,564 G→A 15.9% intergenic (‑222/+69) 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,5640GA15.9% 300.5 / 42.5 125intergenic (‑222/+69)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (52/53);  new base A (6/14);  total (58/67)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

Only 100 of 133 total aligned reads displayed.

GCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCC  >  NC_000913/1579468‑1579663
                                                                                                |                                                                                                   
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)
   aTTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTa                                                                                               <  2:320158/100‑1 (MQ=255)
    ttGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTAt                                                                                              <  2:755064/100‑1 (MQ=255)
          aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCggg                                                                                        >  1:2601729/1‑100 (MQ=255)
           tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGTTGTCGAAACACCCGACATAAATTCCTGTCGGGt                                                                                       <  2:2388897/100‑1 (MQ=255)
           tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                       <  1:240659/100‑1 (MQ=255)
           tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                       <  1:3236333/100‑1 (MQ=255)
           tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                       <  1:325205/100‑1 (MQ=255)
            cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                      <  1:1550950/100‑1 (MQ=39)
            cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                      <  1:3062146/100‑1 (MQ=38)
            cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                      <  2:3306543/100‑1 (MQ=38)
            cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTg                                                                                      <  1:1184602/100‑1 (MQ=255)
              aGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGcc                                                                                    >  1:3058332/1‑100 (MQ=255)
              aGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGcc                                                                                    <  1:647456/100‑1 (MQ=255)
               ggATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCt                                                                                   >  1:3317370/1‑100 (MQ=255)
                gATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCtt                                                                                  <  1:755981/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)
                       cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                           <  2:3252236/100‑1 (MQ=255)
                        tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                          <  1:592632/100‑1 (MQ=255)
                        tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                          <  2:1048950/100‑1 (MQ=255)
                        tAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGt                                                                          >  2:1357285/1‑100 (MQ=255)
                         aaTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                         <  2:3308516/100‑1 (MQ=255)
                           tGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCAt                                                                       >  2:1867896/1‑100 (MQ=255)
                            gCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATa                                                                      <  1:1604219/100‑1 (MQ=255)
                             cccTGATGCATACCACGTAAGCCAGTGTTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCAt                                                                       <  2:2380360/98‑1 (MQ=255)
                             cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                     >  2:1500114/1‑100 (MQ=255)
                              ccTGTTGCATACCACGTAAGCCAGAATTTTCGCAAGGGAAGATGATGATGACTCAACACCCGACAAATTTTTATAACGGATGCCttttaaattt                                                                          >  2:3103918/1‑86 (MQ=255)
                               cTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaaga                                                                   >  2:2702845/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:219357/100‑1 (MQ=255)
                                      aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                            <  1:2895233/100‑1 (MQ=255)
                                       tACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                            >  1:881220/1‑99 (MQ=255)
                                       tACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGACTTTTTAAGGTCATAAGAAGAAGGCt                                                           >  2:2527411/1‑100 (MQ=255)
                                        aCCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTa                                                          <  2:3598148/100‑1 (MQ=255)
                                         ccACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTaa                                                         <  2:2089632/100‑1 (MQ=255)
                                           aCGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAaggagagatcgg                                                       >  2:1849913/1‑93 (MQ=255)
                                           aCGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGa                                                       >  1:2399531/1‑100 (MQ=255)
                                           aCGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGa                                                       >  2:2979790/1‑100 (MQ=255)
                                                aGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGa                                                  <  1:2753640/100‑1 (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)
                                                    aGGATTTTCGCAAGTGAAGATGAAGATGACGAAACACCCAACAGGAATGTATGTCGAGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAgg                                                <  2:3267479/98‑2 (MQ=255)
                                                     ggATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaa                                             >  1:1887480/1‑100 (MQ=255)
                                                     ggATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTCTTAAGGTCATAAGAAGGAGGATAAGAATGACGTTaa                                             >  2:2954225/1‑100 (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)
                                                              cAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCg                                    <  1:3174455/100‑1 (MQ=255)
                                                                aGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGtt                                  >  2:3482849/1‑100 (MQ=255)
                                                                  ggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttat                                >  2:181166/1‑100 (MQ=255)
                                                                  ggAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttat                                <  1:360077/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)
                                                                              tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTcc                    <  1:1115746/100‑1 (MQ=255)
                                                                              tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTcc                    >  2:422152/1‑100 (MQ=255)
                                                                              tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTcc                    <  2:1019712/100‑1 (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)
                                                                                      cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt            >  1:2703015/1‑100 (MQ=39)
                                                                                      cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt            >  1:17683/1‑100 (MQ=38)
                                                                                       acCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtttt           >  2:2313261/1‑100 (MQ=255)
                                                                                        cccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTg          >  1:884853/1‑100 (MQ=255)
                                                                                        cccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTg          >  2:433687/1‑100 (MQ=255)
                                                                                        cccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTg          <  1:1861285/100‑1 (MQ=255)
                                                                                          cGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCa        <  2:3504707/100‑1 (MQ=255)
                                                                                          cGACAAGAACTTATGTCGGGTGCACTGTTAAGGTCATAAGACGGACGCTAAAACTAGAGCTAAGAGCGCGCGATATTAATCACTCTCCCTTTGtttctca        >  2:1617315/1‑96 (MQ=255)
                                                                                           gACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCaa       >  2:1057157/1‑100 (MQ=255)
                                                                                           gACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCaa       >  2:1257482/1‑100 (MQ=255)
                                                                                             cAGGAATTTATGTCGGGTGCCTTGTTAAGGTCGTAAGGAGGAGGCTAAGAATGGAGTTAAGAGAGAGCGTTATTAATTATTCTCCATTTGATATGCAACa     >  2:1958230/1‑100 (MQ=255)
                                                                                             cAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAAc      <  2:3317345/99‑1 (MQ=255)
                                                                                             cAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACa     <  1:1993435/100‑1 (MQ=255)
                                                                                              aGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACAt    >  2:1733978/1‑100 (MQ=255)
                                                                                               ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACAt    >  2:2815528/1‑99 (MQ=255)
                                                                                               ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc   >  1:1735003/1‑100 (MQ=255)
                                                                                               ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc   <  1:2403992/100‑1 (MQ=255)
                                                                                               ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc   <  1:1221482/100‑1 (MQ=255)
                                                                                               ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc   >  2:224740/1‑100 (MQ=255)
                                                                                                gAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATcc  >  2:1445563/1‑100 (MQ=255)
                                                                                                |                                                                                                   
GCTATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCC  >  NC_000913/1579468‑1579663

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: