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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,770,2030TA18.2% 327.2 / 66.9 143intergenic (+17/‑412)ydiK/ydiLUPF0118 family inner membrane protein/putative HTH domain DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (50/66);  new base A (17/9);  total (67/76)
Fisher's exact test for biased strand distribution p-value = 5.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

Only 100 of 148 total aligned reads displayed.

GGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAAT  >  NC_000913/1770104‑1770301
                                                                                                   |                                                                                                  
 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)
      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:2458375/100‑1 (MQ=255)
          agTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTg                                                                                          >  1:2144984/1‑100 (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)
                    ccgACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAacc                                                                                >  1:3319656/1‑100 (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                                                                    >  1:1248453/1‑100 (MQ=255)
                                ccGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                    <  2:859520/100‑1 (MQ=255)
                                ccGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                    <  2:1082944/100‑1 (MQ=255)
                                ccGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGa                                                                    <  1:2688728/100‑1 (MQ=255)
                                  ggAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGata                                                                  <  2:3449703/100‑1 (MQ=255)
                                  ggAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGata                                                                  <  1:3560499/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:960514/1‑100 (MQ=255)
                                      gaaATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttg                                                              <  1:187359/100‑1 (MQ=255)
                                        aaTTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGCTATttgct                                                            <  2:2498932/100‑1 (MQ=255)
                                          ttCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgcttg                                                          <  1:2576280/100‑1 (MQ=255)
                                          ttCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATttgcttg                                                          >  1:2078326/1‑100 (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)
                                                      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)
                                                            aaaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTcc                                        <  1:2962793/100‑1 (MQ=255)
                                                            aaaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCAATATTTTcc                                        <  1:366153/100‑1 (MQ=255)
                                                            aaaTCGAGAAACCGAATAAGTAATTTCGCATCAGGCGGGTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTCTcc                                        <  2:3119175/100‑1 (MQ=255)
                                                             aaTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCt                                       <  1:3166353/100‑1 (MQ=255)
                                                              aTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTg                                      <  1:2893361/100‑1 (MQ=255)
                                                              aTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTg                                      >  1:82209/1‑100 (MQ=255)
                                                              aTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTATTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTg                                      <  1:1290591/100‑1 (MQ=255)
                                                               tCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGc                                     >  2:812376/1‑100 (MQ=255)
                                                                cGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTTCTCTTCCCCATATTTTCCTGCt                                    <  2:2432978/100‑1 (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)
                                                                                 aaTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTgg                   <  2:3299139/100‑1 (MQ=255)
                                                                                  aTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgt                  <  1:2877721/100‑1 (MQ=255)
                                                                                  aTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgt                  <  1:3828347/100‑1 (MQ=255)
                                                                                  aTTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgt                  <  2:3060551/100‑1 (MQ=255)
                                                                                     tctcATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagt               >  1:3391036/1‑100 (MQ=255)
                                                                                      ctcATCAGGCGGCTCTGCCGCCTGAATGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgt                  <  2:1558537/96‑1 (MQ=255)
                                                                                      acaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagta              >  2:1194848/4‑100 (MQ=38)
                                                                                      acaatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagta              >  1:3788773/4‑100 (MQ=38)
                                                                                        catcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc            <  2:1916705/100‑1 (MQ=255)
                                                                                        catcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc            <  1:3011333/100‑1 (MQ=255)
                                                                                        catcaGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc            <  1:683523/100‑1 (MQ=255)
                                                                                        aatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc            >  2:3759171/2‑100 (MQ=255)
                                                                                        aatcaGGCGGCAGAGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTAcc            >  1:3000780/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)
                                                                                               cggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCa     >  2:2614065/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GGGTGGAAGAAGTCCCGCCGCCGACGGACCAACCGGAAGAAATTCTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAAT  >  NC_000913/1770104‑1770301

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: