Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,712,306 (T)7→8 intergenic (+148/‑463) nth → / → dtpA DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,712,2991.T96.5% 298.6 / 3.8 85intergenic (+141/‑470)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base T (36/46);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTT  >  NC_000913/1712201‑1712393
                                                                                                   |                                                                                               
tATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGt                                                                                                >  2:753902/1‑100 (MQ=255)
  tCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttt                                                                                             <  1:1790962/100‑1 (MQ=255)
   ccTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttgt                                                                                             >  2:1236857/1‑98 (MQ=255)
   ccTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttgt                                                                                             >  2:1237248/1‑98 (MQ=255)
   ccTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttt                                                                                            >  2:1233596/1‑100 (MQ=255)
    ctctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttttt                                                                                           >  2:190800/1‑100 (MQ=255)
       ttttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttttt                                                                                        >  2:145017/1‑99 (MQ=255)
       ttttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttttt                                                                                        >  2:1163781/1‑99 (MQ=255)
        tttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttttta                                                                                       <  1:1025754/100‑3 (MQ=255)
            aaaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAg                                                                                    >  1:602416/1‑100 (MQ=255)
            aaaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAg                                                                                    <  1:948108/100‑1 (MQ=255)
             aaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGc                                                                                   <  1:673192/100‑1 (MQ=255)
              aCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCg                                                                                  >  2:311836/1‑100 (MQ=255)
               cATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGa                                                                                 >  1:1369092/1‑100 (MQ=255)
               cATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGa                                                                                 >  1:777980/1‑100 (MQ=255)
                 tCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGaaa                                                                               >  1:536126/1‑100 (MQ=255)
                   tttttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAAtt                                                                             >  2:1712996/1‑100 (MQ=255)
                    ttttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAAttt                                                                            <  2:801718/100‑1 (MQ=255)
                      ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTc                                                                          <  2:1235042/100‑1 (MQ=255)
                      ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTc                                                                          <  2:393307/100‑1 (MQ=255)
                      ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTc                                                                          <  2:1794897/100‑1 (MQ=255)
                      ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTc                                                                          <  1:1974534/100‑1 (MQ=255)
                      ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTc                                                                          <  1:1926748/100‑1 (MQ=255)
                      ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTc                                                                          <  1:1182182/100‑1 (MQ=255)
                      ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTc                                                                          <  1:1182376/100‑1 (MQ=255)
                      ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTc                                                                          <  2:1699738/100‑1 (MQ=255)
                      ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTc                                                                          <  1:552491/100‑1 (MQ=255)
                       tAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCt                                                                         <  1:1279365/100‑1 (MQ=255)
                         aCGTTAATGATTCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTa                                                                        <  2:1606256/99‑1 (MQ=255)
                         aCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTat                                                                       >  2:1625862/1‑100 (MQ=255)
                         aCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTat                                                                       >  2:1622826/1‑100 (MQ=255)
                         aCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTat                                                                       <  1:1048803/100‑1 (MQ=255)
                         aCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTa                                                                        <  1:190164/99‑1 (MQ=255)
                          cGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTatt                                                                      >  1:81241/1‑100 (MQ=255)
                          cGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTatt                                                                      >  1:1687854/1‑100 (MQ=255)
                          cGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTatt                                                                      <  1:1831973/100‑1 (MQ=255)
                            tCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattca                                                                    <  1:238876/100‑1 (MQ=255)
                             cAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattcat                                                                   >  1:762563/1‑100 (MQ=255)
                              aaTGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTACTTGTTTTTTTTAATAGCGAAATTTTCTattcatt                                                                  >  1:1113899/1‑100 (MQ=255)
                              aaTGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattcatt                                                                  <  1:1854181/100‑1 (MQ=255)
                                 gatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGt                                                               >  1:1681766/1‑100 (MQ=255)
                                 gatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGt                                                               >  2:898850/1‑100 (MQ=255)
                                     ccATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATc                                                           >  2:1502126/1‑100 (MQ=255)
                                      cATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCa                                                          <  2:256801/100‑1 (MQ=255)
                                          gCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGAt                                                      <  2:1624620/100‑1 (MQ=255)
                                          gCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTCTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGAt                                                      >  2:1865104/1‑100 (MQ=255)
                                            ttAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAGGATCa                                                    <  1:753905/100‑1 (MQ=255)
                                            ttAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCa                                                    <  1:1163787/100‑1 (MQ=255)
                                             tAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTTTTTTTTTAATTTTCTATTCATTCGTGATCAAGAt                                                      >  2:1563061/1‑97 (MQ=255)
                                              aGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACg                                                  <  1:1168923/100‑1 (MQ=255)
                                              aGCGTTATCATCAGGTAATCCGTTTGCGGATAAACAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACg                                                  <  2:602414/100‑1 (MQ=255)
                                               gCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGa                                                 >  2:865723/1‑100 (MQ=255)
                                               gCGTTATCATCAGGTAATCCATTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGa                                                 >  2:1251086/1‑100 (MQ=255)
                                                cGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGAt                                                >  1:1531912/1‑100 (MQ=255)
                                                  ttATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGc                                              <  1:475873/100‑1 (MQ=255)
                                                  ttATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGc                                              <  1:1713009/100‑1 (MQ=255)
                                                  gatctatcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGc                                              <  1:1317208/95‑1 (MQ=255)
                                                       atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGcagatc                                         >  2:1317200/1‑95 (MQ=255)
                                                         caGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCa                                       <  1:1251092/100‑1 (MQ=255)
                                                          aGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCaa                                      >  1:1546412/1‑100 (MQ=255)
                                                            gTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAAtn                                    <  2:1687841/100‑2 (MQ=255)
                                                            gTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATg                                    >  2:1302234/1‑100 (MQ=255)
                                                             tAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTTTTTTTTAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGt                                   >  2:1601603/1‑100 (MQ=255)
                                                             tAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGt                                   >  1:1624919/1‑100 (MQ=255)
                                                                tCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaa                                >  1:498242/1‑100 (MQ=255)
                                                                  cGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaaaa                              <  1:1064230/100‑1 (MQ=255)
                                                                   gTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAAc                             <  2:1051763/100‑1 (MQ=255)
                                                                   gTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCATGATCAAGATCACAATGCTGTCGCAATGTAAAAAc                             >  1:167416/1‑100 (MQ=255)
                                                                    tttGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACa                            <  2:1051272/100‑1 (MQ=255)
                                                                      tGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAtt                          >  1:1995270/1‑100 (MQ=255)
                                                                      tGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAtt                          >  1:2008232/1‑100 (MQ=255)
                                                                       gCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTg                         >  2:920824/1‑100 (MQ=255)
                                                                        cGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGc                        <  1:1622839/100‑1 (MQ=255)
                                                                        cGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGc                        <  1:1625875/100‑1 (MQ=255)
                                                                          gATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATATCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCaa                      >  1:1976702/1‑100 (MQ=255)
                                                                          gATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTg                         <  2:1012186/97‑1 (MQ=255)
                                                                           aTAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCaa                      <  2:23464/99‑1 (MQ=255)
                                                                             aaCCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGaa                   <  2:589713/100‑1 (MQ=255)
                                                                              aCCAAAAATGCAGGTTAATTGTTTTTTTTAATATCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                  >  2:74630/1‑100 (MQ=255)
                                                                               ccAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCa                >  2:1036661/1‑100 (MQ=255)
                                                                                cAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCa                <  2:777977/100‑1 (MQ=255)
                                                                                  aaaaTGCTTGTTAATTGTTTTTTTTAATAGCGAAATTTTATATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAgg              <  1:145018/100‑1 (MQ=255)
                                                                                  aaaaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAgg              >  2:1300286/1‑100 (MQ=255)
                                                                                   aaaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGt             <  1:525379/100‑1 (MQ=255)
                                                                                   aaaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGt             <  2:1457755/100‑1 (MQ=255)
                                                                                    aaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGtt            >  1:518860/1‑100 (MQ=255)
                                                                                        cAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCCGTCGCAATGTAAAAACATTGCAAGAATCAGGTTaaaa        <  2:1098136/100‑1 (MQ=255)
                                                                                        cAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAAt       <  1:1036666/100‑1 (MQ=255)
                                                                                          ggTTAATTGTTTTTTTTAATAGCGAATTTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATc      >  1:588656/1‑100 (MQ=255)
                                                                                           gTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCa     <  1:454291/100‑1 (MQ=255)
                                                                                            ttAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAt    <  2:1681753/100‑1 (MQ=255)
                                                                                             tAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAtt   >  1:1599889/1‑100 (MQ=255)
                                                                                              aaTTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAttt  >  1:983290/1‑100 (MQ=255)
                                                                                              aaTTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTAATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAttt  <  2:762560/100‑1 (MQ=255)
                                                                                                   |                                                                                               
TATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTT  >  NC_000913/1712201‑1712393

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: