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.T94.8% 382.3 / 12.4 115intergenic (+141/‑470)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base T (58/51);  total (62/53)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

TTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTT  >  NC_000913/1712200‑1712393
                                                                                                    |                                                                                                
ttATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTg                                                                                                   <  1:3329565/100‑1 (MQ=255)
 tATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGt                                                                                                 >  2:2401855/1‑100 (MQ=255)
  aTCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtt                                                                                               >  2:1735974/1‑100 (MQ=255)
   tCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttt                                                                                             <  2:1977459/100‑1 (MQ=255)
   tCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttt                                                                                             <  1:652289/100‑1 (MQ=255)
    ccTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttt                                                                                            >  2:1491769/1‑100 (MQ=255)
     ctctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttttt                                                                                           <  2:766207/100‑1 (MQ=255)
      tctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttt                                                                                          <  2:3314172/100‑1 (MQ=255)
       ctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttttttt                                                                                         >  2:1012811/1‑100 (MQ=255)
        ttttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttttt                                                                                        <  1:223387/100‑2 (MQ=255)
        ntttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttttt                                                                                        >  2:1833619/2‑99 (MQ=255)
         tttAAAACATCTTTTTAACGTCAATGTTGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttttta                                                                                       <  2:3173929/100‑3 (MQ=255)
         tttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATTCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttttta                                                                                       >  1:2445257/1‑98 (MQ=255)
         tttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGtttttttta                                                                                       >  1:120069/1‑98 (MQ=255)
           tAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttttttttaat                                                                                     <  2:420031/100‑5 (MQ=255)
            aaaaCATCTTTTTAACGTCAATGATGCCATTGCTTACCGTTATCATCAGGTAATCCGTTTGCGGCTAAAAAAAAATCCAGCTTACtttttttttttaat                                                                                     >  1:682650/1‑95 (MQ=255)
             aaaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAg                                                                                    <  1:988407/100‑1 (MQ=255)
             aaaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGttttttttaatan                                                                                   <  1:1735950/100‑7 (MQ=255)
              aaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTTATCCGTTTAAAAAAAAAAAAAAATGCAGGTTAATTGTTTTTTTTAATAGc                                                                                   <  1:1174665/100‑1 (MQ=255)
                cATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGa                                                                                 >  2:828498/1‑100 (MQ=255)
                   cTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAAt                                                                              >  1:2839486/1‑100 (MQ=255)
                   cTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAAt                                                                              >  2:3370994/1‑100 (MQ=255)
                   cTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAAt                                                                              >  2:1306632/1‑100 (MQ=255)
                      tttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCAGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAAtttt                                                                           >  2:2005688/1‑100 (MQ=255)
                       ttAACGTCAATGATGCCATTTCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAAACAAAAATGCAGGTTAATTGTTTTTTTTTATAGTTAAAtttt                                                                           >  1:3123375/1‑99 (MQ=255)
                       ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTc                                                                          <  1:2754696/100‑1 (MQ=255)
                          aCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTat                                                                       >  2:2943312/1‑100 (MQ=255)
                          aCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTat                                                                       >  1:2553347/1‑100 (MQ=255)
                            gTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattn                                                                     <  1:1833594/100‑2 (MQ=255)
                            gTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattc                                                                     >  2:1314618/1‑100 (MQ=255)
                             tCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattca                                                                    >  1:237948/1‑100 (MQ=255)
                              cAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTATTAGCTAAATTTTCTattcat                                                                   >  1:1185038/1‑100 (MQ=255)
                               aaTGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattcatt                                                                  >  1:2979721/1‑100 (MQ=255)
                                 tgatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCg                                                                >  2:1932997/1‑100 (MQ=255)
                                 tgatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCg                                                                <  1:2884333/100‑1 (MQ=255)
                                  gatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGt                                                               <  1:2401822/100‑1 (MQ=255)
                                    tgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGa                                                             >  2:1915353/1‑100 (MQ=255)
                                     gCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCTAAATTTTCTATTCATTCGTGAt                                                            >  1:1294481/1‑100 (MQ=255)
                                      ccATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATc                                                           >  2:263820/1‑100 (MQ=255)
                                      ccATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATc                                                           >  2:623300/1‑100 (MQ=255)
                                      ccATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATc                                                           <  2:902369/100‑1 (MQ=255)
                                      ccATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATATCGAAATTTTCTATTCATTCGTGATCa                                                          >  1:1012505/1‑100 (MQ=255)
                                        aTTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCaa                                                         >  2:3163075/1‑100 (MQ=255)
                                          tGCTTAGCGTTATCATCAGGTAATCCGTTTGCTGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGa                                                       <  1:2156130/100‑1 (MQ=255)
                                          tGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCaaca                                                       <  2:2353021/100‑3 (MQ=255)
                                          tGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGa                                                       <  2:2953573/100‑1 (MQ=255)
                                          tGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGa                                                       <  2:2631301/100‑1 (MQ=255)
                                          tGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGa                                                       <  1:566077/100‑1 (MQ=255)
                                          tGCTTAGCGTTATAATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGa                                                       <  1:2462397/100‑1 (MQ=255)
                                           gCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGAt                                                      >  1:102842/1‑100 (MQ=255)
                                           gCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGAt                                                      >  1:2811267/1‑100 (MQ=255)
                                           gCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGAt                                                      >  1:1849526/1‑100 (MQ=255)
                                             ttAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCa                                                    >  2:2428517/1‑100 (MQ=255)
                                             ttAGCGTTATCATCAGGTAATACGTTTAAAAAAAAAAAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATcttgatca                                                    <  1:1259721/100‑8 (MQ=255)
                                              tAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCAc                                                   >  1:2191514/1‑100 (MQ=255)
                                              tAGCGTTATCATCACCTAATCCCTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGTGAAATTTTCTATTCATTCGTGATCAAGATCAc                                                   <  2:3299267/100‑1 (MQ=255)
                                               aGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACg                                                  <  1:734118/100‑1 (MQ=255)
                                               aGCGTTATCATCAGGTAATCCGTTAAAAAAAAAAAAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCAc                                                   <  1:256518/99‑1 (MQ=255)
                                                gCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGa                                                 >  1:1425884/1‑100 (MQ=255)
                                                  gTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCTTGATCACGaat                                               <  2:2793279/100‑3 (MQ=255)
                                                   ttATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGc                                              >  2:2221981/1‑100 (MQ=255)
                                                   ttATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGc                                              <  1:2159163/100‑1 (MQ=255)
                                                    tATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCt                                             >  2:2211591/1‑100 (MQ=255)
                                                     atcatcAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTg                                            <  2:2506648/100‑1 (MQ=255)
                                                      tcatcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCTTGATCACGaatgaat                                           <  1:2495756/100‑7 (MQ=255)
                                                       catcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTCTTTTTTTTTTTTTTTTTTTTCTATTCATTCGTGATCAAGATCACGATGCTGt                                           >  1:1929921/1‑99 (MQ=255)
                                                        atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTTATCAAGATCACGATGCTGTCg                                         >  1:449616/1‑100 (MQ=255)
                                                        atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                         >  2:3101364/1‑100 (MQ=255)
                                                        atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                         <  2:237951/100‑1 (MQ=255)
                                                        atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACCATGCTGTc                                          >  1:1662873/1‑99 (MQ=255)
                                                         tcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGc                                        >  2:2888799/1‑100 (MQ=255)
                                                          caGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGCTGCTGTCGCa                                       >  2:1878579/1‑100 (MQ=255)
                                                          caGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCa                                       >  2:1087219/1‑100 (MQ=255)
                                                              tAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGt                                   >  1:3188634/1‑100 (MQ=255)
                                                               aaTCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTa                                  >  2:2711286/1‑100 (MQ=255)
                                                                   cGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaaa                               <  1:2894451/99‑1 (MQ=255)
                                                                    gTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAAc                             >  2:1487205/1‑100 (MQ=255)
                                                                    gTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAAc                             <  2:681007/100‑1 (MQ=255)
                                                                     tttGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACa                            >  2:3112591/1‑100 (MQ=255)
                                                                          ttaTAACCATAAAGGCAGGTTAATTGTGGTTTTTAATAGCTAAAGGTTCTATGCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCa                       <  2:286681/98‑1 (MQ=255)
                                                                          ggATAACCAAAAATGCAGGTTAATTGTTTTTTTTTTTTTCTAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCa                       >  1:3068681/1‑100 (MQ=255)
                                                                          ggATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCa                       <  2:73731/100‑1 (MQ=255)
                                                                             tAACCAAAAATGCAGGTTAATTGTTGTTTTTAATAGCGAAATTTTCTATTCATTCGTGATTAAGATCAGGATGCTGTCGCAATGTAAAAACATTGCaa                      <  1:2964088/98‑1 (MQ=255)
                                                                               aCCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                  >  1:1003865/1‑100 (MQ=255)
                                                                               aCCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                  >  2:295094/1‑100 (MQ=255)
                                                                                ccAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc                 <  1:2491027/100‑1 (MQ=255)
                                                                                ccAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc                 >  1:2591002/1‑100 (MQ=255)
                                                                                ccAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc                 >  2:1401853/1‑100 (MQ=255)
                                                                                 cAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCa                >  2:1175518/1‑100 (MQ=255)
                                                                                 cAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTCCAAGAATCa                >  1:2750079/1‑100 (MQ=255)
                                                                                 cAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAg               >  1:3453166/1‑100 (MQ=255)
                                                                                   cgatctcAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAgg              <  2:393447/94‑1 (MQ=255)
                                                                                   aaaaTGCAGTTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAgg              <  1:1400216/100‑1 (MQ=255)
                                                                                    aaaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGt             <  1:945848/100‑1 (MQ=255)
                                                                                       tGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGttn           <  1:2831138/99‑2 (MQ=255)
                                                                                         cAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAggagatcg        >  1:393442/1‑94 (MQ=255)
                                                                                          aGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACAATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAAt       >  1:2226124/1‑100 (MQ=255)
                                                                                           ggTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTAGCAATGTAAAAACATTACATTAATCAGGTTaaaa        >  1:2056499/1‑98 (MQ=255)
                                                                                            gTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCa     >  2:1636940/1‑100 (MQ=255)
                                                                                              tAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAtt   <  2:2071609/100‑1 (MQ=255)
                                                                                                    |                                                                                                
TTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTT  >  NC_000913/1712200‑1712393

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: