Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,770,204 C→G 18.8% 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,2040CG18.8% 351.3 / 76.9 144intergenic (+18/‑411)ydiK/ydiLUPF0118 family inner membrane protein/putative HTH domain DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (51/65);  new base G (15/12);  total (67/77)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

Only 100 of 148 total aligned reads displayed.

GTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTC  >  NC_000913/1770106‑1770303
                                                                                                  |                                                                                                   
gtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGctct                                                                                                    >  1:3307503/1‑100 (MQ=255)
    aagaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccg                                                                                                <  1:3532242/100‑1 (MQ=255)
    aagaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccg                                                                                                <  2:3083583/100‑1 (MQ=255)
      gaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgcc                                                                                              <  2:12536/100‑1 (MQ=255)
      gaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgcc                                                                                              <  2:3055873/100‑1 (MQ=255)
         gTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGa                                                                                           >  2:3390794/1‑100 (MQ=255)
          tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAGATCGAGAAACCGAATACGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                          <  2:2518998/100‑1 (MQ=255)
          tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                          <  2:2172925/100‑1 (MQ=255)
          tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                          <  1:2789011/100‑1 (MQ=255)
          tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                          <  2:673661/100‑1 (MQ=255)
          tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                          <  1:1779899/100‑1 (MQ=255)
          tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                          <  1:1804682/100‑1 (MQ=255)
          tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                          <  1:2348412/100‑1 (MQ=255)
           cccGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAtt                                                                                         >  2:3165314/1‑100 (MQ=255)
           cccGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtg                                                                                         <  1:1379982/100‑2 (MQ=255)
            ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTg                                                                                        >  2:3513823/1‑100 (MQ=255)
            ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTg                                                                                        <  2:3406514/100‑1 (MQ=255)
            ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtga                                                                                        <  1:2275528/100‑3 (MQ=255)
            ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtga                                                                                        <  1:1721099/100‑3 (MQ=255)
             cgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGt                                                                                       <  1:3390820/100‑1 (MQ=255)
               ccgccgACGGACCAACCGGAAGTAATACTCGAAGACCCTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtg                                                                                         <  2:3153442/96‑2 (MQ=255)
               ccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAt                                                                                          <  2:2709386/95‑1 (MQ=255)
                  ccgACGGAGCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAacc                                                                                  <  2:1147188/100‑1 (MQ=255)
                  ccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAacc                                                                                  <  2:1903983/100‑1 (MQ=255)
                   cgACGGTCCACCCGGAAGAAATTCTCGAAGAACTTGGCGAGATAGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAt                                                                                          <  2:2657163/91‑1 (MQ=255)
                     aCGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAaccacc                                                                               <  1:499739/100‑1 (MQ=255)
                        gACCAGCCGGAAGAAAGTCTCGAAGAACGTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGc                                                                             <  2:3308049/99‑1 (MQ=255)
                        gACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCt                                                                            <  2:151810/100‑1 (MQ=255)
                         aCCAACCGGAAGAAATTCTCGAATAACTTGGCTAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGTCACTTGTATGTCAACCACCg                                                                              <  2:1702300/97‑1 (MQ=255)
                           cAACCGGAAGAAATTCTCGAAGAACTTGGCGCAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTa                                                                           >  2:3429688/1‑98 (MQ=255)
                              ccGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                      >  1:1836783/1‑100 (MQ=255)
                              ccGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                      <  1:424039/100‑1 (MQ=255)
                                  aagaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATAtt                                                                  <  2:2819133/100‑1 (MQ=255)
                                   agaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCATGATTGTTAACCACCGCTAATTGATAttt                                                                 >  1:3118251/1‑100 (MQ=255)
                                    gaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCATGATTGTTAACCACCGCTAATTGATAttt                                                                 <  2:563621/99‑1 (MQ=255)
                                     aaaTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCCCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgc                                                               >  1:1617764/1‑100 (MQ=255)
                                      aaTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgct                                                              >  1:2212139/1‑100 (MQ=255)
                                      aaTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgct                                                              >  1:863002/1‑100 (MQ=255)
                                       aTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgctt                                                             <  1:1448272/100‑1 (MQ=255)
                                         tctcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgcttgc                                                           <  2:285756/100‑1 (MQ=255)
                                          ctcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTtgcttgct                                                          <  1:1391511/100‑1 (MQ=255)
                                           tcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctc                                                         <  2:110918/100‑1 (MQ=255)
                                            cGAAGAACTTGGCGAGATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                        <  2:2205703/100‑1 (MQ=255)
                                            cGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                        >  1:3440023/1‑100 (MQ=255)
                                            cGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                        >  1:1504453/1‑100 (MQ=255)
                                             naagaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGCTTATTAACCACCGCTAATTTTTATTTGCTTGTTCtt                                                       >  2:1762361/2‑100 (MQ=255)
                                                gaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTccc                                                    <  1:497920/100‑1 (MQ=255)
                                                gaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTccc                                                    <  1:3220822/100‑1 (MQ=255)
                                                      ggCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAtt                                              <  2:1323147/100‑1 (MQ=255)
                                                      ggCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAtt                                              <  2:862993/100‑1 (MQ=255)
                                                        cGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAtttt                                            <  2:2826202/100‑1 (MQ=255)
                                                           aaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTCATTGATATTTGCTTGCTCTTCCCCATATTTTCCt                                         >  2:1215951/1‑100 (MQ=255)
                                                             tCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGc                                       >  1:1196447/1‑100 (MQ=255)
                                                             tCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGc                                       <  1:2125651/100‑1 (MQ=255)
                                                               gagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCtt                                     >  1:139410/1‑100 (MQ=255)
                                                               gagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCtt                                     >  2:2913190/1‑100 (MQ=255)
                                                               gagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCtt                                     >  2:3223069/1‑100 (MQ=255)
                                                               gagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCtt                                     >  1:2125574/1‑100 (MQ=255)
                                                                 gaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTa                                    >  2:922531/1‑99 (MQ=255)
                                                                  aaaCCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTAcc                                  >  2:1381719/1‑100 (MQ=255)
                                                                       gtaTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAt                             <  2:2321811/98‑1 (MQ=255)
                                                                        aaTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAtt                            <  2:3778875/100‑1 (MQ=255)
                                                                         aTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTc                           <  2:163169/100‑1 (MQ=255)
                                                                          tAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCa                          >  1:1347588/1‑100 (MQ=255)
                                                                          tAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCa                          <  1:3566038/100‑1 (MQ=255)
                                                                           aaGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCa                          >  1:3704947/1‑99 (MQ=255)
                                                                           aaGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAg                         <  2:897503/100‑1 (MQ=255)
                                                                           aaGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAg                         <  1:363132/100‑1 (MQ=255)
                                                                            aGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGc                        >  2:1455511/1‑100 (MQ=255)
                                                                            aGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTCCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGc                        >  2:2362728/1‑100 (MQ=255)
                                                                             gTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCt                       <  2:49938/100‑1 (MQ=255)
                                                                             gTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCt                       <  1:2491562/100‑1 (MQ=255)
                                                                              tAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTg                      <  1:1075762/100‑1 (MQ=255)
                                                                              tAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTg                      >  1:1337847/1‑100 (MQ=255)
                                                                              tAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTg                      <  2:1137304/100‑1 (MQ=255)
                                                                               aaTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTgg                     <  1:590568/100‑1 (MQ=255)
                                                                               aaTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTgg                     <  1:826565/100‑1 (MQ=255)
                                                                                 tttCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgta                   <  2:2963949/100‑1 (MQ=255)
                                                                                 tttCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCACATATTTTCCTGCTTACCGCCATCCAGCTGgt                    >  2:3565799/1‑99 (MQ=255)
                                                                                 ttaacaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgta                   >  2:188805/7‑100 (MQ=38)
                                                                                  taacaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtag                  >  1:1879185/6‑100 (MQ=38)
                                                                                   tctcATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAGTCAGCTGgtagt                 >  1:3432190/1‑100 (MQ=255)
                                                                                   aacaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagt                 >  1:1207575/5‑100 (MQ=38)
                                                                                    acaatgaGGCGGCAGAGCCGCCTGATTGTTAAACACCGCTAATTGGTATTTGCTTGCTCTTCCCCATATATTCCTGCTTGCCGGCATTCAGCTGgtagta                >  2:1787840/7‑100 (MQ=255)
                                                                                    acaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagta                >  1:561514/4‑100 (MQ=38)
                                                                                    acaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagta                >  1:2255135/4‑100 (MQ=38)
                                                                                     caatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCGTGCTTTGCGCCATTCATCTTGTAGTAc               >  1:2079805/3‑100 (MQ=255)
                                                                                     caatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAc               >  2:1573949/3‑100 (MQ=255)
                                                                                     caatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCCTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAc               >  2:359427/3‑100 (MQ=255)
                                                                                       atcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt             >  2:114145/1‑100 (MQ=255)
                                                                                       atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt             >  1:1535665/1‑100 (MQ=255)
                                                                                       atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt             >  1:275915/1‑100 (MQ=255)
                                                                                       atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt             >  2:1096329/1‑100 (MQ=255)
                                                                                          aGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTc          <  2:2151156/100‑1 (MQ=255)
                                                                                           ggcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCg         >  1:164247/1‑100 (MQ=255)
                                                                                           ggcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCg         >  1:3822191/1‑100 (MQ=255)
                                                                                            gcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGc        >  1:137079/1‑100 (MQ=255)
                                                                                                ctctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAAt    >  1:1954557/1‑100 (MQ=255)
                                                                                                 tctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCCTTCAGCTGGTAGTACCTGTCGCAAAtt   >  2:1699520/1‑100 (MQ=255)
                                                                                                  ctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCACATTc  >  2:3220653/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTC  >  NC_000913/1770106‑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: