Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,770,204 C→G 19.0% intergenic (+18/‑411) 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,2040CG19.0% 260.6 / 58.8 111intergenic (+18/‑411)ydiK/ydiLUPF0118 family inner membrane protein/putative HTH domain DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (50/40);  new base G (4/17);  total (54/57)
Fisher's exact test for biased strand distribution p-value = 3.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

Only 100 of 119 total aligned reads displayed.

GGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTC  >  NC_000913/1770105‑1770303
                                                                                                   |                                                                                                   
ggtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGctc                                                                                                     >  2:3430435/1‑100 (MQ=255)
ggtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGctc                                                                                                     <  2:1876036/100‑1 (MQ=255)
 gtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGctct                                                                                                    >  1:1014526/1‑100 (MQ=255)
    gaagaaGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgcc                                                                                                 >  1:3222431/1‑100 (MQ=255)
    gaagaaGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgcc                                                                                                 >  1:3221444/1‑100 (MQ=255)
    gaagaaGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgcc                                                                                                 <  1:2248182/100‑1 (MQ=255)
    aaagaagTCCCGCTGCCGACGGACCAACCGGCAGAAATTCTCGAAGAACTCGGCGAAATCGAGAAACCGACTACGTACTTTCCCATTAGGCGGCTCTgc                                                                                                  >  1:3894950/2‑99 (MQ=255)
      agaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgc                                                                                               >  2:451590/1‑100 (MQ=255)
       gaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgcc                                                                                              <  2:2567249/100‑1 (MQ=255)
       gaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgcc                                                                                              >  2:544168/1‑100 (MQ=255)
        aagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCt                                                                                             <  1:4212729/100‑1 (MQ=255)
        aagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCt                                                                                             <  1:1613640/100‑1 (MQ=255)
          gTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGa                                                                                           <  1:1292311/100‑1 (MQ=255)
           tctcgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTTGCGCAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAt                                                                                          <  1:584938/97‑1 (MQ=255)
           tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGa                                                                                           <  1:210857/99‑1 (MQ=255)
           tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAt                                                                                          <  2:2334569/100‑1 (MQ=255)
           tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                          <  1:2665575/100‑1 (MQ=255)
           tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                          <  1:3227959/100‑1 (MQ=255)
           tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                          <  1:2809441/100‑1 (MQ=255)
            cccGCCGCCGATGTACCAACCGGAAGAAATTCTCGAAGAACTCGGCCAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGa                                                                                           <  1:696551/98‑1 (MQ=255)
            cccGCCGCCGACGGACCAACTGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtg                                                                                         <  2:483220/100‑2 (MQ=255)
            cccGCCGCCGACGGACCAACCGGAAGCAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGGAATTTCTCATCAGGCGGCAGAGCCGCCTGAtg                                                                                         <  1:2465948/100‑2 (MQ=255)
            cccGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGGATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtg                                                                                         <  2:3255682/100‑2 (MQ=255)
            cccGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtg                                                                                         <  1:2992251/100‑2 (MQ=255)
             ccgccgccgACGGACCAAGCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtg                                                                                         <  2:4301351/99‑2 (MQ=255)
             ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTg                                                                                        >  2:1078065/1‑100 (MQ=255)
             ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTg                                                                                        <  2:2030981/100‑1 (MQ=255)
             ccgccgacgACGGACCAACCGGAAGAAATTCTCGAAGTACTTGGAGACATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTTCCGCCTGAtag                                                                                        >  1:2058543/1‑98 (MQ=255)
              cgccgccgACGGACCCACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtga                                                                                        <  1:68185/99‑3 (MQ=38)
              cgccgccgACGGACCAACCGGAAGAAATTCTCGAATAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtn                                                                                         <  1:1300557/98‑2 (MQ=255)
              cgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGt                                                                                       <  2:697268/100‑1 (MQ=255)
              cgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGt                                                                                       <  2:2056123/100‑1 (MQ=255)
              cgccgccgACGGACCAAACGGAAGAATTTCTCGAAGATCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtgag                                                                                       <  1:4126577/100‑4 (MQ=37)
                 cgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTaa                                                                                    >  2:4006084/1‑100 (MQ=255)
                 cgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTaa                                                                                    <  2:1893053/100‑1 (MQ=255)
                 cgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTaa                                                                                    >  1:3510855/1‑100 (MQ=255)
                 cgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtga                                                                                        <  1:2734171/96‑3 (MQ=39)
                    cgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAacca                                                                                 >  1:3301330/1‑100 (MQ=255)
                     gACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAacc                                                                                  <  1:322400/98‑1 (MQ=255)
                      aCGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAaccacc                                                                               <  1:3516680/100‑1 (MQ=255)
                        ggACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGc                                                                             <  1:4006057/100‑1 (MQ=255)
                        ggACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGc                                                                             <  2:2045045/100‑1 (MQ=255)
                          aCCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTa                                                                           >  2:1379107/1‑100 (MQ=255)
                             aaCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAtt                                                                        >  2:1917374/1‑100 (MQ=255)
                             aaCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAtt                                                                        >  2:2161557/1‑100 (MQ=255)
                             aaCCGGAAGAAATACTCGAAGAACTTGGCGAGATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAtt                                                                        >  2:2964509/1‑100 (MQ=255)
                              aCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTg                                                                       <  1:3537453/100‑1 (MQ=255)
                                 ggAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGata                                                                    >  2:1216519/1‑100 (MQ=255)
                                 ggAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGata                                                                    <  1:3823358/100‑1 (MQ=255)
                                        aTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgctt                                                             <  2:3224430/100‑1 (MQ=255)
                                        aTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgctt                                                             <  2:2944950/100‑1 (MQ=255)
                                        aTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgctt                                                             <  1:4024402/100‑1 (MQ=255)
                                          tctcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgcttgc                                                           <  1:3107565/100‑1 (MQ=255)
                                           ctcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTtgcttgct                                                          >  2:67226/1‑100 (MQ=255)
                                             cGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                        >  2:374542/1‑100 (MQ=255)
                                              gaagaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCtt                                                       >  2:217289/1‑100 (MQ=255)
                                                 gaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTccc                                                    <  2:3646919/100‑1 (MQ=255)
                                                 gaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTccc                                                    >  1:196024/1‑100 (MQ=255)
                                                 gaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTccc                                                    <  2:1442172/100‑1 (MQ=255)
                                                 gaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTcc                                                     <  1:404503/99‑1 (MQ=255)
                                                     ttGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACGACCGCTAATTGATATTTGCTTGCTCATACCCata                                                >  1:1630875/1‑100 (MQ=255)
                                                        gCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTCATTGATATTTGCTTGCTCTTCCCCATAtt                                              >  1:3555740/1‑99 (MQ=255)
                                                         cGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAtttt                                            >  1:517589/1‑100 (MQ=255)
                                                         cGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAtttt                                            >  2:79610/1‑100 (MQ=255)
                                                          gAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTc                                           >  1:4144738/1‑100 (MQ=255)
                                                             aTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTg                                        >  1:3789555/1‑100 (MQ=255)
                                                              tCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGc                                       <  1:2141630/100‑1 (MQ=255)
                                                              tCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGc                                       <  2:2860405/100‑1 (MQ=255)
                                                               cGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCt                                      >  2:3668927/1‑100 (MQ=255)
                                                               cGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCt                                      >  2:2935746/1‑100 (MQ=255)
                                                               agagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGACATTTGCTTGCTCCTCCCCATATTTTCCTGCt                                      >  1:1611457/2‑100 (MQ=255)
                                                                gagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCtt                                     >  2:488341/1‑100 (MQ=255)
                                                                gagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCtt                                     >  2:1093340/1‑100 (MQ=255)
                                                                 agaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTa                                    <  1:226733/100‑1 (MQ=255)
                                                                 agaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTa                                    >  2:3179925/1‑100 (MQ=255)
                                                                      ccGAATAAGTAATTTCTCATCAGGTGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGcc                               >  1:2684252/1‑100 (MQ=255)
                                                                      ccGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGcc                               <  1:2958061/100‑1 (MQ=255)
                                                                        gAATAAGTAATTTCTCATCAGGCGGCTCTGCCGTCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAt                             <  1:1154303/100‑1 (MQ=255)
                                                                          aTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTc                           <  1:2924499/100‑1 (MQ=255)
                                                                              gTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCt                       <  1:4316272/100‑1 (MQ=255)
                                                                              gTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCt                       >  2:2281170/1‑100 (MQ=255)
                                                                                  tttCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgta                   >  2:2612721/1‑100 (MQ=255)
                                                                                  tttCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgta                   <  1:909195/100‑1 (MQ=255)
                                                                                    tctcATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagt                 >  1:1880085/1‑100 (MQ=255)
                                                                                     acaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagta                >  1:439145/4‑100 (MQ=38)
                                                                                      tcatcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCACTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAc               >  1:2811444/1‑100 (MQ=255)
                                                                                      tcatcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCATGCTCTTCCACATATATTCCTGCTCACCGCCATTCAGCTGGTAGTAc               >  1:2976928/1‑100 (MQ=255)
                                                                                       catcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc              >  1:3459018/1‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt             >  1:1150158/1‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt             >  2:3877555/1‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt             >  2:2785352/1‑100 (MQ=255)
                                                                                           aGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTc          <  2:112881/100‑1 (MQ=255)
                                                                                            ggcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCg         >  2:3750900/1‑100 (MQ=255)
                                                                                            ggcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCg         >  2:2495007/1‑100 (MQ=255)
                                                                                               ggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaa      >  1:618403/1‑100 (MQ=255)
                                                                                                gcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaaa     <  2:4075715/100‑1 (MQ=255)
                                                                                                 ctctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAAt    <  2:336419/100‑1 (MQ=255)
                                                                                                 ctctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAAt    <  2:2259590/100‑1 (MQ=255)
                                                                                                  tctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAAtt   <  2:3510873/100‑1 (MQ=255)
                                                                                                   ctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTc  <  1:1216514/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTC  >  NC_000913/1770105‑1770303

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: