Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913857,5760TC20.8% 248.6 / 63.0 120intergenic (+21/+220)ybiT/ybiUABC‑F family putative regulatory ATPase/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base T (43/52);  new base C (14/11);  total (57/63)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-01

Only 100 of 120 total aligned reads displayed.

CTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCA  >  NC_000913/857478‑857672
                                                                                                  |                                                                                                
    aaaTCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGt                                                                                             >  2:1784649/1‑100 (MQ=255)
        cACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCggg                                                                                         >  2:2644699/1‑100 (MQ=255)
        cACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCggg                                                                                         <  2:2164219/100‑1 (MQ=255)
         aCCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCgggg                                                                                        <  2:1637919/100‑1 (MQ=255)
          ccccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                       <  2:2630205/100‑1 (MQ=255)
          ccccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                       >  1:1059325/1‑100 (MQ=255)
          ccccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGt                                                                                       >  1:757614/1‑100 (MQ=255)
           cccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTg                                                                                      >  2:308669/1‑100 (MQ=255)
           cccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTg                                                                                      <  1:1571642/100‑1 (MQ=255)
           cccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTg                                                                                      <  1:1065018/100‑1 (MQ=255)
            ccGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGa                                                                                     >  2:1126528/1‑100 (MQ=255)
             cGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAc                                                                                    <  2:875275/100‑1 (MQ=39)
             cGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAc                                                                                    <  1:1170673/100‑1 (MQ=39)
             cGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTAGCCGTCACCCCGACAGTATCGGGGTGAc                                                                                    <  1:1953186/100‑1 (MQ=38)
              ggAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGAcg                                                                                   >  1:326385/1‑100 (MQ=255)
              ggAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAcg                                                                                   <  1:1972609/100‑1 (MQ=39)
              ggAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAcg                                                                                   <  2:1551132/100‑1 (MQ=39)
              ggAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAcg                                                                                   <  1:383705/100‑1 (MQ=39)
              ggAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAc                                                                                    <  2:535286/99‑1 (MQ=255)
                   gcgTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACggcggc                                                                              <  2:2523355/100‑1 (MQ=255)
                   gcgTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACggcggc                                                                              <  2:192265/100‑1 (MQ=255)
                     gTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGACAGTATCGGGGTGAcggataatc                                                                            <  2:362608/100‑7 (MQ=37)
                      tGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtga                                                                           <  2:1450460/100‑1 (MQ=255)
                      tGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtga                                                                           >  2:2656303/1‑100 (MQ=255)
                       gATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCtgat                                                                          >  1:1344696/1‑100 (MQ=255)
                          cgccTTTAGCGGTAATTACGACGATTCCCCGCGTCGTACCGGGATCGAGTCAGTGATTATCCGTCACCCCGTTAATGTCGGGGCGTCGGCGCCTGATTAc                                                                       <  2:2246279/97‑1 (MQ=255)
                          cGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGAc                                                                       <  2:423160/100‑1 (MQ=255)
                           gACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACa                                                                      <  2:1595262/100‑1 (MQ=255)
                             cTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACaaa                                                                    <  2:1799749/100‑1 (MQ=255)
                             cTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACa                                                                      <  2:1880073/98‑1 (MQ=255)
                               ttAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATg                                                                  <  1:1684401/100‑1 (MQ=255)
                               ttAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATg                                                                  <  2:844496/100‑1 (MQ=255)
                                  gCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGcga                                                               >  1:392431/1‑98 (MQ=255)
                                  gCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaa                                                               >  1:2018720/1‑100 (MQ=255)
                                    ggTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCaaaa                                                             >  2:2484869/1‑100 (MQ=255)
                                        aTTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGc                                                         <  2:864833/100‑1 (MQ=255)
                                        aTTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGc                                                         <  1:2619397/100‑1 (MQ=255)
                                         ttACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGcc                                                        <  1:2091523/100‑1 (MQ=255)
                                          tACGAAGATTACCTGCGTAGTAAAGGGATCGTGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATATcc                                                        >  2:1281376/1‑99 (MQ=255)
                                          tACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                       >  1:2214300/1‑100 (MQ=255)
                                          tACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCt                                                       >  1:389497/1‑100 (MQ=255)
                                           aCGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTg                                                      >  2:635515/1‑100 (MQ=255)
                                            cGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGa                                                     >  1:972134/1‑100 (MQ=255)
                                                  ttACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCt                                               >  1:956612/1‑100 (MQ=255)
                                                    aCCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTAc                                             <  2:1860967/100‑1 (MQ=255)
                                                    aCCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTAc                                             >  1:971157/1‑100 (MQ=255)
                                                     ccTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACg                                            <  1:953409/100‑1 (MQ=255)
                                                       tGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCt                                          <  2:2226612/100‑1 (MQ=255)
                                                       tGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCt                                          <  1:1394511/100‑1 (MQ=255)
                                                       tGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCt                                          >  2:719807/1‑100 (MQ=255)
                                                        gCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCtt                                         >  1:1342196/1‑100 (MQ=255)
                                                         cGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTa                                        >  2:1793725/1‑100 (MQ=255)
                                                         cGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTa                                        <  1:261771/100‑1 (MQ=255)
                                                           tagtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATc                                      <  2:951426/100‑1 (MQ=255)
                                                           tagtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATc                                      >  1:1085948/1‑100 (MQ=255)
                                                             gtaAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAg                                    >  2:735854/1‑100 (MQ=255)
                                                                aaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGTGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGcc                                 <  1:1874959/100‑1 (MQ=255)
                                                                aaGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGcc                                 <  2:1335886/100‑1 (MQ=255)
                                                                  gggATCGAGTAAGTGATTATGCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGGCTa                               >  2:1646387/1‑100 (MQ=255)
                                                                  gggATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCCGGCCTa                               >  1:528748/1‑100 (MQ=255)
                                                                  gggATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTa                               >  1:726430/1‑100 (MQ=255)
                                                                   ggATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTAc                              >  2:951840/1‑100 (MQ=255)
                                                                   ggATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTAc                              >  2:946469/1‑100 (MQ=255)
                                                                   ggATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTAc                              >  1:2656051/1‑100 (MQ=255)
                                                                     aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGc                            <  1:197317/100‑1 (MQ=255)
                                                                     aTCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGc                            >  1:2131797/1‑100 (MQ=255)
                                                                      tCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCa                           <  2:757607/100‑1 (MQ=255)
                                                                        gAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGt                         <  2:2325141/100‑1 (MQ=255)
                                                                         aGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTa                        <  2:1872007/100‑1 (MQ=255)
                                                                          gTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAc                       >  2:1948674/1‑100 (MQ=255)
                                                                           tAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTAcc                      <  1:677541/100‑1 (MQ=255)
                                                                             aGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTg                    <  1:2318425/100‑1 (MQ=255)
                                                                              gTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGc                   <  1:234073/100‑1 (MQ=255)
                                                                              gTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGc                   <  1:2484894/100‑1 (MQ=255)
                                                                              gTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGc                   <  1:1874615/100‑1 (MQ=255)
                                                                               tGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTATGCAGTACCTGCa                  <  2:1134158/100‑1 (MQ=255)
                                                                                gATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCaa                 >  1:2271371/1‑100 (MQ=255)
                                                                                 aTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAat                >  2:1102189/1‑100 (MQ=255)
                                                                                 aTTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAat                <  1:1537025/100‑1 (MQ=255)
                                                                                  gccgccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAata               >  2:1096843/5‑100 (MQ=25)
                                                                                  gccgccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAata               >  2:2054499/5‑100 (MQ=25)
                                                                                  gccgccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAata               >  1:705591/5‑100 (MQ=25)
                                                                                  gccgccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCACTACGCTTATCAGGCCTACGCAGTACCTGCAata               >  2:1421122/5‑100 (MQ=25)
                                                                                     tCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat            <  2:2600822/100‑1 (MQ=255)
                                                                                     tCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat            <  2:971149/100‑1 (MQ=255)
                                                                                     tCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat            <  1:2501899/100‑1 (MQ=255)
                                                                                     gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat            >  2:897618/2‑100 (MQ=25)
                                                                                     gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatatat            >  1:406424/2‑100 (MQ=25)
                                                                                     gccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAatat              >  2:2552816/2‑98 (MQ=25)
                                                                                      ccGTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATAtt           >  1:2407105/1‑100 (MQ=25)
                                                                                       cGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTg          >  1:251234/1‑100 (MQ=255)
                                                                                        gTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa         >  1:2601365/1‑100 (MQ=25)
                                                                                        gTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa         >  1:2413064/1‑100 (MQ=25)
                                                                                        gTCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGa         >  2:2495716/1‑100 (MQ=25)
                                                                                         tCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGaa        >  1:2610162/1‑100 (MQ=25)
                                                                                         tCACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGaa        >  1:1098728/1‑100 (MQ=25)
                                                                                          cACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAAt       <  2:1059317/100‑1 (MQ=255)
                                                                                          cACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAAt       <  1:1784669/100‑1 (MQ=255)
                                                                                          cACCCCGACAGTATCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAAt       >  1:1487530/1‑100 (MQ=25)
                                                                                                gATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCa  <  2:2379436/99‑1 (MQ=255)
                                                                                                  |                                                                                                
CTGGAAATCACCCCGGAACGCGTGATCGACTTTAGCGGTAATTACGAAGATTACCTGCGTAGTAAAGGGATCGAGTAAGTGATTATCCGTCACCCCGATACTGTCGGGGTGACGGCGGCTGATGACAAATGCAAAATAGCCTGATGCGCTACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCA  >  NC_000913/857478‑857672

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑