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.T95.6% 321.6 / 6.8 91intergenic (+141/‑470)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base T (46/41);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 108 total aligned reads displayed.

TATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAT  >  NC_000913/1712201‑1712391
                                                                                                   |                                                                                            
 aTCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtt                                                                                            >  1:1670072/1‑100 (MQ=255)
  tCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttt                                                                                           >  1:271050/1‑100 (MQ=255)
    ntctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttttt                                                                                         >  1:2308409/2‑100 (MQ=255)
    ctctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttttt                                                                                         >  1:631277/1‑100 (MQ=255)
    ctctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttttt                                                                                         <  2:1414461/100‑1 (MQ=255)
     tctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttt                                                                                        <  1:1812874/100‑1 (MQ=255)
      ctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACAAAAAATGCAGGTTAATTGttttttt                                                                                       <  2:29417/100‑1 (MQ=255)
       ttttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttttt                                                                                      <  1:2335411/100‑2 (MQ=255)
       ttttAAAACATCTTTTTAAAGTCAATTATGCCATTGCTTAGCGTTATTATCAGGTAATCCGTTTGCGCATAACCAAAAATGCAGGTTAAttttttttttt                                                                                      <  1:350231/100‑2 (MQ=255)
        tttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttttta                                                                                     <  2:1935660/100‑3 (MQ=255)
        tttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttttta                                                                                     <  1:2268277/100‑3 (MQ=255)
        tttAAAACATCTTTTTAACGTCAATGATGCCATTACTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttttta                                                                                     <  2:1935322/100‑3 (MQ=255)
         ttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttttttttaa                                                                                    <  1:306152/100‑4 (MQ=255)
          tAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttttttttaat                                                                                   <  1:946566/100‑5 (MQ=255)
          tAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTAAAAAAAAAAAAAAATGCAGGTTAATTGttttttttaan                                                                                   <  2:1448450/100‑5 (MQ=255)
              aCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCg                                                                                <  2:1004262/100‑1 (MQ=255)
               cATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGa                                                                               >  1:877590/1‑100 (MQ=255)
               cATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGa                                                                               >  1:182177/1‑100 (MQ=255)
               cATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGa                                                                               >  2:2268310/1‑100 (MQ=255)
                 tCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGaaa                                                                             <  2:1460272/100‑1 (MQ=255)
                  cTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAgatcgga                                                                            >  1:544631/1‑94 (MQ=255)
                      ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAAttttn                                                                        <  1:861358/100‑2 (MQ=255)
                       tAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCt                                                                       <  2:182175/100‑1 (MQ=255)
                       tAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCt                                                                       <  1:1181410/100‑1 (MQ=255)
                        atcGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTa                                                                      <  1:1211936/98‑1 (MQ=255)
                        aaCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTat                                                                     <  2:1903549/100‑1 (MQ=255)
                         aCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTat                                                                     >  2:1880479/1‑100 (MQ=255)
                           gTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattc                                                                   >  2:1118729/1‑100 (MQ=255)
                           gTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattc                                                                   >  2:1778447/1‑100 (MQ=255)
                           gTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattc                                                                   >  2:2077895/1‑100 (MQ=255)
                            tCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattca                                                                  >  2:902192/1‑100 (MQ=255)
                             cAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattcat                                                                 <  2:1908235/100‑1 (MQ=255)
                                tgatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCg                                                              <  2:707480/100‑1 (MQ=255)
                                 gatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGt                                                             <  1:710584/100‑1 (MQ=255)
                                  atgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTg                                                            <  2:588151/100‑1 (MQ=255)
                                   tgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGa                                                           <  2:2236347/100‑1 (MQ=255)
                                     ccATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATc                                                         >  1:1006899/1‑100 (MQ=255)
                                      cATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCa                                                        <  2:2257351/100‑1 (MQ=255)
                                      cATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCa                                                        >  2:720387/1‑100 (MQ=255)
                                      cATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCa                                                        <  2:1125964/100‑1 (MQ=255)
                                          gCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGAt                                                    >  1:606829/1‑100 (MQ=255)
                                          gCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGAt                                                    >  2:1769851/1‑100 (MQ=255)
                                          gCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGAt                                                    >  2:823798/1‑100 (MQ=255)
                                           cTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATc                                                   >  2:1972419/1‑100 (MQ=255)
                                           cTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATc                                                   >  2:794123/1‑100 (MQ=255)
                                             tAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCAc                                                 <  1:349175/100‑1 (MQ=255)
                                             tAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCAc                                                 <  1:1769856/100‑1 (MQ=255)
                                               gCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGa                                               >  2:1635467/1‑100 (MQ=255)
                                                   tATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCt                                           >  1:628537/1‑100 (MQ=255)
                                                      catcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTc                                        >  1:615446/1‑100 (MQ=255)
                                                      catcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTc                                        <  1:2077902/100‑1 (MQ=255)
                                                       atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                       <  2:537850/100‑1 (MQ=255)
                                                       atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                       <  2:631276/100‑1 (MQ=255)
                                                       atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGc                                      >  1:194356/1‑100 (MQ=255)
                                                        tcaGGTAATCCGTTTGCGGATAATCAACAATGCAGGTTAATTTTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                       <  1:610954/99‑1 (MQ=255)
                                                        tcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGc                                      >  2:1177076/1‑100 (MQ=255)
                                                         caGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGc                                      <  2:1338944/99‑1 (MQ=255)
                                                         caGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCa                                     <  1:2132869/100‑1 (MQ=255)
                                                         caGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGCCGCa                                     >  1:737219/1‑100 (MQ=255)
                                                           ggTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAAt                                   >  1:697368/1‑100 (MQ=255)
                                                             tAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGt                                 >  2:439532/1‑100 (MQ=255)
                                                             tAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGt                                 >  2:25235/1‑100 (MQ=255)
                                                              aaTCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTa                                <  1:1442878/100‑1 (MQ=255)
                                                               aTCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaa                               <  1:1956727/100‑1 (MQ=255)
                                                                 ccGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaaa                             >  2:147870/1‑100 (MQ=255)
                                                                 ccGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGttttta                            <  2:194354/100‑6 (MQ=255)
                                                                  cGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaaaa                            <  1:1635472/100‑1 (MQ=255)
                                                                   gTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAAc                           >  2:690600/1‑100 (MQ=255)
                                                                   gTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAAc                           >  2:1542006/1‑100 (MQ=255)
                                                                   gTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAAc                           <  1:439534/100‑1 (MQ=255)
                                                                      tGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAtt                        >  1:1773309/1‑100 (MQ=255)
                                                                       gCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTg                       >  2:452784/1‑100 (MQ=255)
                                                                        cGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGc                      >  1:1293172/1‑100 (MQ=255)
                                                                        cGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGc                      >  2:897443/1‑100 (MQ=255)
                                                                          gATAACCAAAAATGCTGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCaa                    <  1:1485038/100‑1 (MQ=255)
                                                                          gATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGc                      <  2:753905/98‑1 (MQ=255)
                                                                             aaCCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGaa                 <  1:147872/100‑1 (MQ=255)
                                                                             aaCCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGaa                 >  2:1978253/1‑100 (MQ=255)
                                                                              aCCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                >  1:395805/1‑100 (MQ=255)
                                                                              aCCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                >  1:2197258/1‑100 (MQ=255)
                                                                               ccAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc               >  1:924324/1‑100 (MQ=255)
                                                                               ccAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc               <  2:50020/100‑1 (MQ=255)
                                                                               ccAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc               <  2:1293169/100‑1 (MQ=255)
                                                                               ccAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc               >  1:2248350/1‑100 (MQ=255)
                                                                               ccAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc               >  1:771429/1‑100 (MQ=255)
                                                                               ccAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc               <  2:661710/100‑1 (MQ=255)
                                                                               ccAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGAGCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc               >  2:414595/1‑100 (MQ=255)
                                                                               ccAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGAACAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc               >  1:688610/1‑100 (MQ=255)
                                                                                cAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCa              <  2:1128691/100‑1 (MQ=255)
                                                                                 aaaaaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATTAAGATTACGATGCTGTCGCAATGTAAAAACATTGTAAGAATTAg             >  2:434354/1‑100 (MQ=255)
                                                                                 aaaaaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAg             >  1:548479/1‑100 (MQ=255)
                                                                                  aaaaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAgg            <  1:33636/100‑1 (MQ=255)
                                                                                   aaaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAggn           <  2:648918/100‑2 (MQ=255)
                                                                                   aaaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGt           <  2:548099/100‑1 (MQ=255)
                                                                                      tGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTaa        <  2:526974/100‑1 (MQ=255)
                                                                                      tGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTaa        <  1:897445/100‑1 (MQ=255)
                                                                                           gTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCa   >  1:717013/1‑100 (MQ=255)
                                                                                           gTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCa   >  1:1213583/1‑100 (MQ=255)
                                                                                            ttAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAt  <  1:720388/100‑1 (MQ=255)
                                                                                            ttAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAt  >  1:1491525/1‑100 (MQ=255)
                                                                                                   |                                                                                            
TATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAT  >  NC_000913/1712201‑1712391

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: