Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,770,205 T→A 18.7% intergenic (+19/‑410) 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,2050TA18.7% 334.7 / 69.7 144intergenic (+19/‑410)ydiK/ydiLUPF0118 family inner membrane protein/putative HTH domain DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (52/65);  new base A (15/12);  total (67/77)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

Only 100 of 150 total aligned reads displayed.

GTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCT  >  NC_000913/1770106‑1770304
                                                                                                   |                                                                                                   
gtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGctct                                                                                                     >  1:3307503/1‑100 (MQ=255)
    aagaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccg                                                                                                 <  1:3532242/100‑1 (MQ=255)
    aagaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccg                                                                                                 <  2:3083583/100‑1 (MQ=255)
      gaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgcc                                                                                               <  2:3055873/100‑1 (MQ=255)
        agTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTACCGCCTg                                                                                             >  2:1265033/1‑100 (MQ=255)
         gTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGa                                                                                            >  2:3390794/1‑100 (MQ=255)
          tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAGATCGAGAAACCGAATACGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                           <  2:2518998/100‑1 (MQ=255)
          tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                           <  2:673661/100‑1 (MQ=255)
          tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                           <  1:2789011/100‑1 (MQ=255)
          tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                           <  1:1779899/100‑1 (MQ=255)
           cccGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAtt                                                                                          >  2:3165314/1‑100 (MQ=255)
           cccGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAtt                                                                                          >  1:226716/1‑100 (MQ=255)
            ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTg                                                                                         <  2:3406514/100‑1 (MQ=255)
            ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTg                                                                                         >  2:3513823/1‑100 (MQ=255)
            ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTg                                                                                         <  2:1667424/100‑1 (MQ=255)
            ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtga                                                                                         <  1:2275528/100‑3 (MQ=255)
             cgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGt                                                                                        <  1:3390820/100‑1 (MQ=255)
               ccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAt                                                                                           <  2:2709386/95‑1 (MQ=255)
               ccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTa                                                                                      <  2:1157015/100‑1 (MQ=255)
                cgccgAGGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAGTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAt                                                                                           <  2:2657048/94‑1 (MQ=255)
                     aCGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAaccacc                                                                                <  1:499739/100‑1 (MQ=255)
                        gACCAGCCGGAAGAAAGTCTCGAAGAACGTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGc                                                                              <  2:3308049/99‑1 (MQ=255)
                           cAACCGGAAGAAATTCTCGAAGAACTTGGCGCAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTa                                                                            >  2:3429688/1‑98 (MQ=255)
                              ccGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                       >  1:1836783/1‑100 (MQ=255)
                              ccGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                       <  1:424039/100‑1 (MQ=255)
                              ccGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                       <  1:1446832/100‑1 (MQ=255)
                                  aagaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATAtt                                                                   >  1:163169/1‑100 (MQ=255)
                                   agaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCATGATTGTTAACCACCGCTAATTGATAttt                                                                  >  1:3118251/1‑100 (MQ=255)
                                    gaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCATGATTGTTAACCACCGCTAATTGATAttt                                                                  <  2:563621/99‑1 (MQ=255)
                                      aaTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgct                                                               >  1:2212139/1‑100 (MQ=255)
                                      aaTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgct                                                               >  1:863002/1‑100 (MQ=255)
                                       aTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgctt                                                              <  1:1448272/100‑1 (MQ=255)
                                       aTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCAGCCTGATTGTTAACCTCCGCTAATTGATATttgctt                                                              >  2:3089712/1‑100 (MQ=255)
                                         tctcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgcttgc                                                            <  2:285756/100‑1 (MQ=255)
                                          ctcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTtgcttgct                                                           <  1:1391511/100‑1 (MQ=255)
                                           tcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctc                                                          <  2:110918/100‑1 (MQ=255)
                                            cGAAGAACTTGGTGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                         >  1:151488/1‑100 (MQ=255)
                                            cGAAGAACTTGGCGAGATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                         <  2:2205703/100‑1 (MQ=255)
                                            cGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                         <  2:1388041/100‑1 (MQ=255)
                                            cGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                         >  1:3440023/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)
                                                    ttGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCata                                                 >  2:1519687/1‑100 (MQ=255)
                                                     tGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCatat                                                <  1:1936103/100‑1 (MQ=255)
                                                      ggCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAtt                                               <  2:1323147/100‑1 (MQ=255)
                                                      ggCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAtt                                               <  2:862993/100‑1 (MQ=255)
                                                       gCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAttt                                              >  2:1715610/1‑100 (MQ=255)
                                                        cGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAtttt                                             <  2:2826202/100‑1 (MQ=255)
                                                         gAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTc                                            >  2:1792659/1‑100 (MQ=255)
                                                           aaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTCATTGATATTTGCTTGCTCTTCCCCATATTTTCCt                                          >  2:1215951/1‑100 (MQ=255)
                                                            aTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTTCTCTTCCCCATATTTTCCTg                                         <  2:2338726/100‑1 (MQ=255)
                                                             tCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGc                                        >  1:1196447/1‑100 (MQ=255)
                                                               gagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCtt                                      >  2:3223069/1‑100 (MQ=255)
                                                               gagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCtt                                      >  1:2125574/1‑100 (MQ=255)
                                                               gagaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCtt                                      >  2:2913190/1‑100 (MQ=255)
                                                                 gaAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTa                                     >  2:922531/1‑99 (MQ=255)
                                                                  aaaCCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTAcc                                   <  1:1381728/100‑1 (MQ=255)
                                                                        aaTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAtt                             <  1:1349469/100‑1 (MQ=255)
                                                                        aaTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAtt                             <  1:2097300/100‑1 (MQ=255)
                                                                        aaTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAtt                             <  2:3778875/100‑1 (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                        <  1:2491562/100‑1 (MQ=255)
                                                                             gTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCt                        <  2:49938/100‑1 (MQ=255)
                                                                             gTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCt                        >  2:2200512/1‑100 (MQ=255)
                                                                              tAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTg                       <  1:1075762/100‑1 (MQ=255)
                                                                              tAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTg                       <  1:2136572/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)
                                                                                  ttCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGTTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtag                   <  1:2708996/100‑1 (MQ=255)
                                                                                  taacaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtag                   >  2:2143985/6‑100 (MQ=38)
                                                                                   tctcATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAGTCAGCTGgtagt                  >  1:3432190/1‑100 (MQ=255)
                                                                                    acaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagta                 >  1:561514/4‑100 (MQ=38)
                                                                                     caatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAc                >  2:1573949/3‑100 (MQ=255)
                                                                                     caatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCCTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAc                >  2:359427/3‑100 (MQ=255)
                                                                                      aatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc               >  1:126232/2‑100 (MQ=255)
                                                                                       atcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt              >  2:114145/1‑100 (MQ=255)
                                                                                       atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt              >  1:1535665/1‑100 (MQ=255)
                                                                                       atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt              >  2:1096329/1‑100 (MQ=255)
                                                                                       atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt              >  1:275915/1‑100 (MQ=255)
                                                                                          aGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTc           <  2:2151156/100‑1 (MQ=255)
                                                                                           ggcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTTCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCg          >  2:2713296/1‑100 (MQ=255)
                                                                                           ggcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCg          >  1:164247/1‑100 (MQ=255)
                                                                                           ggcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCg          >  1:3822191/1‑100 (MQ=255)
                                                                                            gcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGc         >  1:137079/1‑100 (MQ=255)
                                                                                               gcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaaa      >  1:1830441/1‑100 (MQ=255)
                                                                                               gcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATAGTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCaaa      <  1:1761949/100‑1 (MQ=255)
                                                                                                ctctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAAt     >  1:1954557/1‑100 (MQ=255)
                                                                                                 tctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCCTTCAGCTGGTAGTACCTGTCGCAAAtt    >  2:1699520/1‑100 (MQ=255)
                                                                                                  ctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCACATTc   >  2:3220653/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCT  >  NC_000913/1770106‑1770304

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: