Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,579,569 T→C 16.3% intergenic (‑227/+64) 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,5690TC16.3% 389.3 / 67.5 172intergenic (‑227/+64)yddA/ydeMputative multidrug ABC transporter permease/ATPase/putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base T (70/73);  new base C (13/15);  total (83/89)
Fisher's exact test for biased strand distribution p-value = 8.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-02

Only 100 of 177 total aligned reads displayed.

ATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATG  >  NC_000913/1579471‑1579669
                                                                                                  |                                                                                                    
aTTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTa                                                                                                     <  1:532673/100‑1 (MQ=255)
 ttGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTAt                                                                                                    <  2:578259/100‑1 (MQ=255)
  tGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATg                                                                                                   <  1:751804/100‑1 (MQ=255)
   gATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGt                                                                                                  >  2:1789357/1‑100 (MQ=255)
       aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCGTGTCggg                                                                                              <  2:2307758/100‑1 (MQ=255)
       aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCggg                                                                                              <  1:384714/100‑1 (MQ=255)
       aTCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCgg                                                                                               >  1:3798320/1‑99 (MQ=255)
        tcgaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGACGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGt                                                                                             >  1:3681208/4‑100 (MQ=255)
        tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                             <  1:3220606/100‑1 (MQ=255)
        tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGt                                                                                             <  2:1376089/100‑1 (MQ=255)
        tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGt                                                                                             >  2:1652665/1‑100 (MQ=255)
        tCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGt                                                                                             <  2:4391788/100‑1 (MQ=255)
         cAAGGATATTGCTAATGCCCTGCTGCATACCCCGTAATCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                            <  1:2970419/100‑1 (MQ=37)
         cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                            <  2:2150737/100‑1 (MQ=39)
         cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                            <  2:661691/100‑1 (MQ=38)
         cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCggg                                                                                              >  1:2756659/1‑98 (MQ=255)
         cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTg                                                                                            <  1:2053851/100‑1 (MQ=255)
         cAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTg                                                                                            <  1:3974866/100‑1 (MQ=255)
         cAAGGATATTGCTAATGCCCTGATGCATAACACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGt                                                                                             <  1:4046534/99‑1 (MQ=255)
          aaGGATATTTCTAATGCCGTGGTGCATACCACGTAAGCCAGGATTTGCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTg                                                                                            <  1:1125375/99‑1 (MQ=37)
          aaGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGTTTTTCGCAAGGGAAGATGCTGCTGGCGAAACACCCGACAGGAATTTATGTCGGGTGc                                                                                           <  1:3262206/100‑1 (MQ=255)
              atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACATAAATTCCTGTCGGGTgtttcg                                                                                       <  1:3134287/100‑6 (MQ=37)
              atatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTg                                                                                       <  1:3634592/100‑1 (MQ=255)
               tatTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGt                                                                                      <  1:1652658/100‑1 (MQ=255)
                 ttGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTa                                                                                    <  1:962576/100‑1 (MQ=255)
                 ttGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTa                                                                                    >  2:2548146/1‑100 (MQ=255)
                   gCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACCCCCGACAGGAATTTATGTCGGGTACCTTGTTCAg                                                                                  >  2:1465794/1‑100 (MQ=255)
                   gCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAg                                                                                  <  2:2834837/100‑1 (MQ=255)
                    cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                                 <  2:4038217/100‑1 (MQ=255)
                    cTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                                 <  1:776860/100‑1 (MQ=255)
                    cTAATGCCATGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAgg                                                                                 <  1:1230412/100‑1 (MQ=255)
                      aaTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTc                                                                               >  2:4102686/1‑100 (MQ=255)
                       aTGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCa                                                                              <  1:196061/100‑1 (MQ=255)
                          cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                           <  2:85822/100‑1 (MQ=255)
                          cccTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaa                                                                           >  1:691808/1‑100 (MQ=255)
                            cTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGGTGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaaga                                                                         >  1:2370703/1‑100 (MQ=255)
                             tGATGCATACCTCGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGCATTTATGTCGGGTGCCTTGTTAAGGTCATaagaa                                                                        >  1:2318914/1‑100 (MQ=255)
                              gATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATaagaag                                                                       <  1:3977834/100‑1 (MQ=255)
                                   aTACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGc                                                                  >  2:3054431/1‑100 (MQ=255)
                                    tACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCt                                                                 >  2:2529227/1‑100 (MQ=255)
                                     aCCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTa                                                                >  2:3379926/1‑100 (MQ=255)
                                     aCCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTa                                                                <  1:666028/100‑1 (MQ=255)
                                      ccACGTAAGCCAGGATATTCGCAGGGTGAGATGATGATGACGACACACCCGACAGGAAATTATGTCGGATGCCTTATTAAGGTCATAAGAAGGTGGCTa                                                                >  1:150408/1‑99 (MQ=255)
                                       cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAg                                                              <  2:2369655/100‑1 (MQ=255)
                                       cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAg                                                              <  1:3684189/100‑1 (MQ=255)
                                       cACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAg                                                              <  2:691809/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)
                                              gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                       <  2:281032/100‑1 (MQ=255)
                                              gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                       <  2:2410836/100‑1 (MQ=255)
                                              gCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAg                                                       >  2:3851948/1‑100 (MQ=255)
                                               ccAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGt                                                      >  1:1627243/1‑100 (MQ=255)
                                                 aGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTa                                                    <  1:2529214/100‑1 (MQ=255)
                                                    aTTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTAAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTaaaa                                                 <  2:3113079/100‑1 (MQ=255)
                                                       ttCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAagag                                              <  2:3427146/100‑1 (MQ=255)
                                                          gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGc                                           <  2:1434344/100‑1 (MQ=255)
                                                          gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGc                                           >  1:1752299/1‑100 (MQ=255)
                                                          gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGc                                           >  2:1342366/1‑100 (MQ=255)
                                                           cAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCg                                          >  2:531022/1‑100 (MQ=255)
                                                                gAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatt                                     <  1:1469408/100‑1 (MQ=255)
                                                                  aGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGttatta                                    <  1:3859912/99‑1 (MQ=255)
                                                                   gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaat                                  >  2:3665859/1‑100 (MQ=255)
                                                                   gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaat                                  <  1:3665834/100‑1 (MQ=255)
                                                                     tgatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGTAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatta                                >  2:1186679/1‑100 (MQ=255)
                                                                     tgatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatta                                <  2:4063266/100‑1 (MQ=255)
                                                                     tgatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTattaatt                                 >  1:1818551/1‑99 (MQ=255)
                                                                      gatgatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAattat                               <  2:4234679/100‑1 (MQ=255)
                                                                         gatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtct                            <  2:269160/100‑1 (MQ=255)
                                                                         gatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtct                            >  2:315431/1‑100 (MQ=255)
                                                                         gatgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtct                            <  1:2803928/100‑1 (MQ=255)
                                                                          atgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATtctc                           >  2:2149378/1‑100 (MQ=255)
                                                                           tgaCGAAACACCCGACAGGAATTTATGTTGGGTGCCTTGTTAAGGTCATAAGAAAGAGGCTCAGATTGTCCTTAAAAGAGAGCGTTATTAattat                               >  1:3724439/1‑95 (MQ=255)
                                                                           tgaCGAAACACCCGACAGGAATTTATGTCGTGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGCATGGAGTTAAAAGAGCGCGGTATTAATTATGCTcc                          >  1:1198868/1‑100 (MQ=255)
                                                                           tgaCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTcc                          >  2:747548/1‑100 (MQ=255)
                                                                             aCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCAt                        <  2:719540/100‑1 (MQ=255)
                                                                                aaaCACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTg                     <  2:2060481/100‑1 (MQ=255)
                                                                                 ggcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGt                    >  2:1789875/3‑100 (MQ=38)
                                                                                  gcacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGtt                   >  2:2755134/2‑100 (MQ=38)
                                                                                   cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt                  >  1:4282195/1‑100 (MQ=38)
                                                                                   cacCCGACATAAATTCCTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttt                  >  1:337706/1‑100 (MQ=38)
                                                                                     cccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTg                <  1:2329301/100‑1 (MQ=255)
                                                                                      ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGttttn                <  1:891519/99‑2 (MQ=255)
                                                                                      ccGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGc               >  2:949836/1‑100 (MQ=255)
                                                                                       cGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCa              <  2:2281580/100‑1 (MQ=255)
                                                                                       cGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCa              >  1:3384795/1‑100 (MQ=255)
                                                                                        gACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAATGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCaa             <  2:4054195/100‑1 (MQ=255)
                                                                                          cAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACa           >  2:359034/1‑100 (MQ=255)
                                                                                            ggAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGATTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATc         >  2:516111/1‑100 (MQ=255)
                                                                                             gAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATcc        >  2:2551530/1‑100 (MQ=255)
                                                                                              aaTTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGTAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCa       <  1:2140636/100‑1 (MQ=255)
                                                                                              aaTTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCa       >  2:190727/1‑100 (MQ=255)
                                                                                              aaTTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCa       >  1:992244/1‑100 (MQ=255)
                                                                                              aaTTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCa       >  1:3570399/1‑100 (MQ=255)
                                                                                               aTTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCAt      >  2:4217575/1‑100 (MQ=255)
                                                                                               aTTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCAt      >  2:3360875/1‑100 (MQ=255)
                                                                                                tttATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATa     >  2:4209474/1‑100 (MQ=255)
                                                                                                  tATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTCAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAAt   >  1:1050294/1‑100 (MQ=255)
                                                                                                  tATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAAt   <  2:2257743/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
ATTGATAATCAAGGATATTGCTAATGCCCTGATGCATACCACGTAAGCCAGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTTGCAACATCCATAATG  >  NC_000913/1579471‑1579669

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