Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,770,204 C→G 18.2% 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.2% 349.3 / 74.6 143intergenic (+18/‑411)ydiK/ydiLUPF0118 family inner membrane protein/putative HTH domain DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (50/67);  new base G (17/9);  total (67/76)
Fisher's exact test for biased strand distribution p-value = 4.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

Only 100 of 148 total aligned reads displayed.

GGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTC  >  NC_000913/1770105‑1770303
                                                                                                   |                                                                                                   
ggtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGCATAAGTAATTTCTCATCCGGCGGctc                                                                                                     >  2:1232054/1‑100 (MQ=255)
ggtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGctc                                                                                                     <  1:3888761/100‑1 (MQ=255)
 gtggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAGATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGctct                                                                                                    >  1:1709894/1‑100 (MQ=255)
  tggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTg                                                                                                   >  1:201683/1‑100 (MQ=255)
   ggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGATAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgc                                                                                                  <  2:2083023/100‑1 (MQ=255)
   ggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgc                                                                                                  >  2:2893352/1‑100 (MQ=255)
   ggAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgc                                                                                                  >  1:2875369/1‑100 (MQ=255)
     aagaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccg                                                                                                <  2:71260/100‑1 (MQ=255)
       gaagTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTgccgcc                                                                                              <  2:111464/100‑1 (MQ=255)
        aagTCCCGCGGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCt                                                                                             <  2:2765543/100‑1 (MQ=255)
         agTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCt                                                                                             <  2:3231127/99‑1 (MQ=255)
         agTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTg                                                                                            >  1:2144984/1‑100 (MQ=255)
         agTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTg                                                                                            <  1:2458375/100‑1 (MQ=255)
          gTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTgn                                                                                           <  1:1465387/100‑2 (MQ=255)
          gTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGa                                                                                           <  1:1194133/100‑1 (MQ=255)
           tCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAt                                                                                          <  1:668424/100‑1 (MQ=255)
            cccGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGAtt                                                                                         >  1:2468336/1‑100 (MQ=255)
             ccgccgccgACGGACCAGCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtga                                                                                        <  1:1815687/100‑3 (MQ=255)
             ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGCAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTg                                                                                        >  2:2325693/1‑100 (MQ=255)
             ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGACTg                                                                                        >  2:3909328/1‑100 (MQ=255)
             ccgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtga                                                                                        <  1:2583880/100‑3 (MQ=255)
              cgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGt                                                                                       >  2:3011320/1‑100 (MQ=255)
              cgccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCAGAGCCGCCTGAtgag                                                                                       <  2:2412677/100‑4 (MQ=38)
               gccgccgACGCACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTACTTTCTCATCAGGCGGCAGAGCCGCCTGAtgaga                                                                                      <  2:2752683/100‑5 (MQ=37)
                ccgccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTa                                                                                     <  1:873497/100‑1 (MQ=255)
                  gccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTa                                                                                     <  2:2228199/98‑1 (MQ=255)
                  gccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAac                                                                                   >  2:3828322/1‑100 (MQ=255)
                  gccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAac                                                                                   <  1:4155604/100‑1 (MQ=255)
                    cgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAacca                                                                                 >  2:2157684/1‑100 (MQ=255)
                             aaCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAtt                                                                        >  2:2994030/1‑100 (MQ=255)
                             aaCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAAtt                                                                        <  2:556109/100‑1 (MQ=255)
                              aCCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTg                                                                       <  2:1421243/100‑1 (MQ=255)
                               ccGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                      <  2:859520/100‑1 (MQ=255)
                               ccGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                      <  2:1082944/100‑1 (MQ=255)
                               ccGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                      >  1:1248453/1‑100 (MQ=255)
                                 ggAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGata                                                                    <  1:3560499/100‑1 (MQ=255)
                                 ggAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGata                                                                    <  2:3449703/100‑1 (MQ=255)
                                 ggAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAAAAAGTAATTTCTCATCAGTCGGCTCAGCCGCCTGATTGTTAACCCCCGCTAATTGata                                                                    >  1:38962/1‑100 (MQ=255)
                                   aagaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATAtt                                                                  >  1:2773821/1‑100 (MQ=255)
                                     gaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttg                                                                >  1:2970604/1‑100 (MQ=255)
                                     gaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttg                                                                <  1:187359/100‑1 (MQ=255)
                                     gaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttg                                                                >  1:960514/1‑100 (MQ=255)
                                       aaTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGCTATttgct                                                              <  2:2498932/100‑1 (MQ=255)
                                         ttCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgcttg                                                            >  1:2078326/1‑100 (MQ=255)
                                         ttCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgcttg                                                            <  1:2576280/100‑1 (MQ=255)
                                            tcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctc                                                         >  1:302457/1‑100 (MQ=255)
                                             cGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgcttgc                                                           >  2:4023640/1‑97 (MQ=255)
                                             cGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                        >  1:1517716/1‑100 (MQ=255)
                                              gaagaaCTGGGCGGAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGACTGTTAACCACCGCTAATTGATATTTGCTTGctct                                                        <  2:1071315/99‑1 (MQ=255)
                                                 gaaCTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTccc                                                    <  2:2859132/100‑1 (MQ=255)
                                                    cTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCat                                                 >  1:3756460/1‑100 (MQ=255)
                                                     ttagCGAAATCTAGAAACCAAATAACTAATTCGTCATAAGGCGGCTCTGCCGTTTTACTGTTAACCACCGCTAATTGATATGTGCTTGCTCTTCGCCata                                                <  2:2820851/97‑1 (MQ=255)
                                                     ttGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCata                                                <  1:306962/100‑1 (MQ=255)
                                                     ttGACGAAATCGAGAAACCGTATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCa                                                  <  1:771914/98‑1 (MQ=255)
                                                        gCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAttt                                             >  2:4099960/1‑100 (MQ=255)
                                                        gCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAttt                                             >  1:3226609/1‑100 (MQ=255)
                                                           aaaTCGAGAAACCGAATAAGTAATTTCTCATCGGGCGGCGCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTcc                                          <  2:2941980/100‑1 (MQ=255)
                                                           aaaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAAGAACCGCTAATTGATATTTGCTTGCTCTTTCCCATATTTTcc                                          >  1:2353485/1‑100 (MQ=255)
                                                           aaaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTcc                                          <  2:3076816/100‑1 (MQ=255)
                                                           aaaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCAATATTTTcc                                          <  1:366153/100‑1 (MQ=255)
                                                           aaaTCGAGAAACCGAATAAGTAATTTCGCATCAGGCGGGTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTCTcc                                          <  2:3119175/100‑1 (MQ=255)
                                                             aTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTg                                        >  1:82209/1‑100 (MQ=255)
                                                             aTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTg                                        <  1:2893361/100‑1 (MQ=255)
                                                             aTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTg                                        <  1:2995057/100‑1 (MQ=255)
                                                             aTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTATTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTg                                        <  1:1290591/100‑1 (MQ=255)
                                                              tCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGc                                       >  2:812376/1‑100 (MQ=255)
                                                               cGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTTCTCTTCCCCATATTTTCCTGCt                                      <  2:2432978/100‑1 (MQ=255)
                                                               cGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCt                                      >  1:2956917/1‑100 (MQ=255)
                                                                     aCCGAATAAGTAATTTCTCATCAGGCGGTTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGc                                >  2:591523/1‑100 (MQ=255)
                                                                     aCCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGc                                <  2:426399/100‑1 (MQ=255)
                                                                       cGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATCTTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCa                              >  2:4092039/1‑100 (MQ=255)
                                                                       cGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCa                              <  1:1611167/100‑1 (MQ=255)
                                                                       cGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCa                              <  1:1610411/100‑1 (MQ=255)
                                                                       cGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTCCTTCCTCTTCCCCATATTTTCCTGCTCACCGCCa                              >  2:3562447/1‑100 (MQ=255)
                                                                        gAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACTGCCAt                             >  2:3475871/1‑100 (MQ=255)
                                                                         aaTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCAtt                            <  2:1587666/100‑1 (MQ=255)
                                                                          aTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTc                           <  1:2331547/100‑1 (MQ=255)
                                                                              gTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCt                       >  2:2598009/1‑100 (MQ=255)
                                                                               tAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTg                      <  1:3230600/100‑1 (MQ=255)
                                                                                aaTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTgg                     <  2:3299139/100‑1 (MQ=255)
                                                                                 aTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgt                    <  2:3060551/100‑1 (MQ=255)
                                                                                 aTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgt                    <  1:3828347/100‑1 (MQ=255)
                                                                                    tctcATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagt                 >  1:3391036/1‑100 (MQ=255)
                                                                                     ctcATCAGGCGGCTCTGCCGCCTGAATGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgt                    <  2:1558537/96‑1 (MQ=255)
                                                                                     acaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagta                >  1:3788773/4‑100 (MQ=38)
                                                                                     acaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagta                >  2:1194848/4‑100 (MQ=38)
                                                                                       catcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc              <  2:1916705/100‑1 (MQ=255)
                                                                                       catcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc              <  1:683523/100‑1 (MQ=255)
                                                                                       catcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc              <  1:3011333/100‑1 (MQ=255)
                                                                                       aatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc              >  2:3759171/2‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAGTTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt             >  2:4182900/1‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTCCCt             >  2:708241/1‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt             >  1:963632/1‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt             >  1:726621/1‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt             >  1:3443889/1‑100 (MQ=255)
                                                                                        atcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCt             >  1:559254/1‑100 (MQ=255)
                                                                                         tcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTg            <  2:3251393/100‑1 (MQ=255)
                                                                                           aGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTc          <  2:204779/100‑1 (MQ=255)
                                                                                             gcggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGc        >  1:3413113/1‑100 (MQ=255)
                                                                                              cggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCa       >  2:2614065/1‑100 (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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: