Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,770,203 T→A 15.7% intergenic (+17/‑412) ydiK → / → ydiL UPF0118 family inner membrane protein/putative HTH domain DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,770,2030TA15.7% 248.5 / 38.0 102intergenic (+17/‑412)ydiK/ydiLUPF0118 family inner membrane protein/putative HTH domain DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (39/47);  new base A (8/8);  total (47/55)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.93e-01

Only 100 of 104 total aligned reads displayed.

GGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGC  >  NC_000913/1770105‑1770297
                                                                                                  |                                                                                              
ggtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGctc                                                                                               >  2:1180946/1‑100 (MQ=255)
 gtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGctct                                                                                              <  2:655290/100‑1 (MQ=255)
  tggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTg                                                                                             >  1:1687612/1‑100 (MQ=255)
   ggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgc                                                                                            >  2:374957/1‑100 (MQ=255)
   ggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgc                                                                                            <  2:646567/100‑1 (MQ=255)
    gaagaaGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCAGCTCTgcc                                                                                           >  1:2006252/1‑100 (MQ=255)
     aagaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccg                                                                                          <  2:2118831/100‑1 (MQ=255)
     aagaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccg                                                                                          <  2:2118782/100‑1 (MQ=255)
        aagTCCCGCCGCCGACGGACGAACCGGAAGAAATTCTCGAAGAACTTGTCGAAATCGAGAAACCGAATAAGTAATTTCTCCTCAGGCGGCTCTgccgcc                                                                                        <  1:284312/99‑1 (MQ=255)
        aagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCt                                                                                       >  2:2363983/1‑100 (MQ=255)
         agTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCt                                                                                       <  1:911767/99‑1 (MQ=255)
         agTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTg                                                                                      <  2:395818/100‑1 (MQ=255)
          gTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGTAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGa                                                                                     <  1:1749876/100‑1 (MQ=255)
          gTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGa                                                                                     <  2:2437273/100‑1 (MQ=255)
           tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                    <  2:2404842/100‑1 (MQ=255)
           tCCCGCCGCCGACGGACCAACCGGAAGAAAGTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                    <  1:384944/100‑1 (MQ=255)
            cccGCCGCCGACGGACCATCCGGAATAAATTCTCGAAGAACTTGGCGAGATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtg                                                                                   <  1:462470/100‑2 (MQ=255)
            cccGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtg                                                                                   <  2:1395201/100‑2 (MQ=255)
            cccGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtg                                                                                   <  1:2374400/100‑2 (MQ=255)
             ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAt                                                                                    <  1:1608745/98‑1 (MQ=255)
             ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTg                                                                                  >  2:2209211/1‑100 (MQ=255)
             ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTg                                                                                  >  1:2332929/1‑100 (MQ=255)
             ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTg                                                                                  <  2:931169/100‑1 (MQ=255)
             ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTAATTg                                                                                  >  1:173876/1‑100 (MQ=255)
             ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATAAGGCGGCTCTGAAGCCTGATTg                                                                                  >  2:1120002/1‑100 (MQ=255)
             ccgccgccgACGGACCAACCAGAAGATATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAATAATTTCTCATCAGGCGGCAGAGCCGCCTGAtg                                                                                   <  1:1556599/99‑2 (MQ=255)
                ccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTa                                                                               <  2:1637920/100‑1 (MQ=255)
                 cgccgACGGACCAACCGGAAGAAATTCTCCAAGCACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTa                                                                               <  2:690783/99‑1 (MQ=255)
                  gccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAac                                                                             >  1:120351/1‑100 (MQ=255)
                   cgatctgACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAacc                                                                            <  2:2539360/94‑1 (MQ=255)
                   ccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAacc                                                                            <  1:242774/100‑1 (MQ=255)
                    cgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAacca                                                                           >  1:67315/1‑100 (MQ=255)
                     gacttgccAGCCGTAAGAAATTGTCGAAGAACTTGGCGAAATCGAGAAACCGAATACGTAATTTCTCATCAGGCGGCTCTGCCGCATGATTGTTAaccac                                                                          >  1:1125916/7‑100 (MQ=255)
                      aCGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAaccacc                                                                         <  2:1874072/100‑1 (MQ=255)
                      aCGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAaccacc                                                                         >  2:727347/1‑100 (MQ=255)
                        ggACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGc                                                                       <  2:2209681/100‑1 (MQ=255)
                        ggACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGc                                                                       <  2:1444339/100‑1 (MQ=255)
                         gACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAaccagatcg                                                                      >  1:2539341/1‑95 (MQ=255)
                         gACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCt                                                                      >  2:2147246/1‑100 (MQ=255)
                           ccAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCtt                                                                     <  2:1030547/99‑2 (MQ=255)
                           ccAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTaa                                                                    <  2:1687621/100‑1 (MQ=255)
                             aaCCGGAAGAAATTCTCGAAGCACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAtt                                                                  >  2:833827/1‑100 (MQ=255)
                             aaCAGGAAGACATTCTCGAAGAACTTGGCGAACTAGCGACACCGAATAAGTAATCTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAta                                                                  <  1:519736/100‑2 (MQ=255)
                              aCCGGAAGAAATTCTCGCAGAACTTGGCGCACTCGAGAAACCGAATAAGTCATTTCTCATCAGGAGGCTCTGCCGCCTGATTGTTCACACCCACTACTTg                                                                 >  1:2371034/1‑100 (MQ=255)
                              aCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTg                                                                 >  1:1722617/1‑100 (MQ=255)
                                 ggAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGata                                                              <  2:909486/100‑1 (MQ=255)
                                 ggAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGata                                                              >  1:792053/1‑100 (MQ=255)
                                     gaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttg                                                          <  1:80126/100‑1 (MQ=255)
                                       aaTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgct                                                        >  2:2564122/1‑100 (MQ=255)
                                        aTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgctt                                                       <  2:2057025/100‑1 (MQ=255)
                                          tctcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgcttgc                                                     <  2:2494635/100‑1 (MQ=255)
                                           ctcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTtgcttgct                                                    >  1:1166146/1‑100 (MQ=255)
                                             cGAAGAACTTGGCGACATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                  <  2:2550830/100‑1 (MQ=255)
                                             cGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                  >  2:843130/1‑100 (MQ=255)
                                                  aaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTcccc                                             >  2:1601389/1‑100 (MQ=255)
                                                   aCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCa                                            >  2:1201648/1‑100 (MQ=255)
                                                      tGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCatat                                         <  2:1935309/100‑1 (MQ=255)
                                                         cGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAtttt                                      >  2:1930695/1‑100 (MQ=255)
                                                           aaaTCTAGAAACCGAATAAGGAGTGTCTGATCAGACGGCTCTGCCGCCTTATTGTTACCCGCAGCTCATTGATATTTGCTTGCTCGTCCCCATGTTGTcc                                    <  1:437545/100‑1 (MQ=255)
                                                           aaaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTcc                                    <  2:1318272/100‑1 (MQ=255)
                                                           aaaTCGAGAAACCAAATCAGTAATTTCTCATCAGGCGGCTCTACCGCCGAATTGTTAACCACCGCTCATTGATATTTGCTTGCTCTTCCCCata                                          >  1:2002562/1‑94 (MQ=255)
                                                            aaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCt                                   <  1:2552233/100‑1 (MQ=255)
                                                            aaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCt                                   <  1:51591/100‑1 (MQ=255)
                                                            aaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCt                                   <  1:1914238/100‑1 (MQ=255)
                                                             aTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTg                                  >  2:176450/1‑100 (MQ=255)
                                                                gagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCtt                               >  2:2314828/1‑100 (MQ=255)
                                                                    aaCCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCg                           <  2:218963/100‑1 (MQ=255)
                                                                     aCCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCg                           <  1:885488/99‑1 (MQ=255)
                                                                     aCCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGc                          >  2:68516/1‑100 (MQ=255)
                                                                     aCCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGc                          <  1:4257/100‑1 (MQ=255)
                                                                     aCCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGc                          <  2:558136/100‑1 (MQ=255)
                                                                       cGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCa                        <  2:635073/100‑1 (MQ=255)
                                                                        gAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAt                       <  2:1903001/100‑1 (MQ=255)
                                                                        gAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAt                       >  2:1018284/1‑100 (MQ=255)
                                                                         aaTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAtt                      <  2:1376062/100‑1 (MQ=255)
                                                                           tAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCa                    >  2:2138785/1‑100 (MQ=255)
                                                                            aaGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAg                   <  2:1303917/100‑1 (MQ=255)
                                                                             aGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTCCCGCCATTCAGc                  >  2:979504/1‑100 (MQ=255)
                                                                             aGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGc                  >  1:407739/1‑100 (MQ=255)
                                                                                aaTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTgg               <  1:2240148/100‑1 (MQ=255)
                                                                                aaTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTgg               >  2:315259/1‑100 (MQ=255)
                                                                                 aTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgt              <  2:164361/100‑1 (MQ=255)
                                                                                 aTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgt              <  2:1587117/100‑1 (MQ=255)
                                                                                  ttaacaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgta             >  1:2573073/7‑100 (MQ=38)
                                                                                    tctcATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagt           >  2:20968/1‑100 (MQ=255)
                                                                                      tcatcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAc         <  2:1256266/100‑1 (MQ=255)
                                                                                      tcatcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAc         <  1:2146175/100‑1 (MQ=255)
                                                                                      caatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAc         >  2:1345654/3‑100 (MQ=255)
                                                                                      caatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAc         >  2:1965320/3‑100 (MQ=255)
                                                                                       catcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc        >  2:855543/1‑100 (MQ=255)
                                                                                       catcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc        <  1:1665136/100‑1 (MQ=255)
                                                                                       aatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc        >  2:1975564/2‑100 (MQ=255)
                                                                                        atcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt       >  1:805656/1‑100 (MQ=255)
                                                                                        atcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt       >  2:2175945/1‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt       >  2:1201882/1‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt       >  1:35487/1‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt       >  1:1138162/1‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCACCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt       >  1:623057/1‑100 (MQ=255)
                                                                                           aGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTc    <  1:1075476/100‑1 (MQ=255)
                                                                                             gcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGc  <  2:2273568/100‑1 (MQ=255)
                                                                                                  |                                                                                              
GGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGC  >  NC_000913/1770105‑1770297

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: