Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,712,306:1 (T)7→8 100% 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% 310.8 / 9.8 97intergenic (+141/‑470)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base T (40/52);  total (42/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 113 total aligned reads displayed.

TTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑TTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTT  >  NC_000913/1712200‑1712393
                                                                                                    |                                                                                                
ttATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTg                                                                                                   <  2:2056826/100‑1 (MQ=255)
ttATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTg                                                                                                   <  1:1594180/100‑1 (MQ=255)
   tCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑ttt                                                                                               >  1:3200548/1‑100 (MQ=255)
    ccTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCATGTAATCCGTTTGCGGGTAACCAAAAAGGCAGGTTAATTG‑tttt                                                                                              >  2:1385239/1‑100 (MQ=255)
       ctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑ttttt‑tt                                                                                          >  2:1049740/1‑100 (MQ=255)
        ttttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑ttttt‑ttt                                                                                         >  1:3636514/1‑99 (MQ=255)
        ttttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑ttttt‑ttt                                                                                         <  2:1092596/100‑2 (MQ=255)
          ttAATACATCTTTTTTACGTCAATGAGGCCATTGCTTAGCTTTATCATCAGGTAATCTGTTTGCGGATAACCAAAAATGCAGGGTAATTG‑ttttt‑ttt                                                                                         <  2:360212/98‑2 (MQ=255)
          ttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑ttttt‑tttaa                                                                                       <  1:2325260/100‑4 (MQ=255)
           tAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCATGTAATCCGTTTGCAGATAACCAAAAATGCAGGTTAATTG‑ttttt‑tttaat                                                                                      <  2:1618387/100‑5 (MQ=255)
           tAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑ttttt‑tttaat                                                                                      <  1:2013198/100‑5 (MQ=255)
           aaaaaCATCTTTTTAACGTCAATGATGCCTTTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGACAACCAAAAATGCAGGTTAATTG‑ttttt‑tttaat                                                                                      <  2:2500313/99‑5 (MQ=255)
            aaaaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑ttttt‑tttaata                                                                                     >  1:3490414/1‑95 (MQ=255)
             aaaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAg                                                                                    >  1:1614513/1‑100 (MQ=255)
             aaaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAg                                                                                    >  2:3600709/1‑100 (MQ=255)
             aaaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAg                                                                                    <  1:3436756/100‑1 (MQ=255)
                cATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCg                                                                                  >  2:3409125/1‑99 (MQ=255)
                cATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGa                                                                                 >  1:2298442/1‑100 (MQ=255)
                 aTCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGaa                                                                                >  1:4122234/1‑100 (MQ=255)
                  tCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGaaa                                                                               <  2:3258516/100‑1 (MQ=255)
                       ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTc                                                                          <  1:1920724/100‑1 (MQ=255)
                        tAACGTCATTGATGCCATTGCTTAGCTTTATCACCGCGTAATCCGTTTGCGGATAACCAAAAATGCATGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTc                                                                          <  2:1299560/99‑1 (MQ=255)
                        tAACGTCAATGATGCCATTGCTTAGCGTTATCGTCAGGTGATTCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCt                                                                         <  2:2924720/100‑1 (MQ=255)
                        tAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCt                                                                         >  2:441306/1‑100 (MQ=255)
                        tAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCt                                                                         <  1:2341165/100‑1 (MQ=255)
                         aaCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCTAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCt                                                                         <  2:2272692/99‑1 (MQ=255)
                          aCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTat                                                                       >  1:430052/1‑100 (MQ=255)
                           cGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTatt                                                                      <  1:2276882/100‑1 (MQ=255)
                           cGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTatt                                                                      <  1:3161677/100‑1 (MQ=255)
                           cGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTatt                                                                      >  1:969919/1‑100 (MQ=255)
                           cGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTatt                                                                      >  1:2943383/1‑100 (MQ=255)
                           cGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTatt                                                                      <  2:2808582/100‑1 (MQ=255)
                              cAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTattcat                                                                   >  1:1247195/1‑100 (MQ=255)
                               aaTGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTattcatt                                                                  <  2:1248556/100‑1 (MQ=255)
                                atgatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTattcattc                                                                 >  1:4026529/1‑100 (MQ=255)
                                 tgatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCg                                                                <  2:3577366/100‑1 (MQ=255)
                                 tgatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCg                                                                <  1:1664240/100‑1 (MQ=255)
                                  gatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGt                                                               <  2:1619701/100‑1 (MQ=255)
                                  gatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGt                                                               <  1:1009484/100‑1 (MQ=255)
                                   atgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGAGGATAAACAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTg                                                              <  2:1462735/100‑1 (MQ=255)
                                    tgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGa                                                             >  2:1807140/1‑100 (MQ=255)
                                    tgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGa                                                             <  1:2235533/100‑1 (MQ=255)
                                     gCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGAt                                                            >  1:72102/1‑100 (MQ=255)
                                     gCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGAt                                                            <  2:2643588/100‑1 (MQ=255)
                                     gCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGAt                                                            <  2:2978901/100‑1 (MQ=255)
                                     gCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGAt                                                            <  1:3075010/100‑1 (MQ=255)
                                        aTTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCaa                                                         <  1:750795/100‑1 (MQ=255)
                                           gattAGCTTTATCATCAGGTAATCCGTTTGGTGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGAt                                                      <  2:4165235/98‑1 (MQ=255)
                                            cTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATcttgatc                                                     <  1:2523809/100‑7 (MQ=255)
                                             ttAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATc                                                     <  2:64574/99‑1 (MQ=255)
                                             ttAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCa                                                    <  1:405110/100‑1 (MQ=255)
                                              tAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCTTGATCAc                                                   <  1:1188811/100‑1 (MQ=255)
                                               aGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACg                                                  >  2:1970567/1‑100 (MQ=255)
                                               aGCGTTATCATCAGGTAATCCGTTTGCGGACAAGCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGAGTAAGATCACg                                                  >  2:2738904/1‑100 (MQ=255)
                                                gCGTTATCATCAGGTAATCCGTTTGCGGATAACCAACAATGCAGGTTAATTGTTTTTT‑TT‑ATTAGCGATTTTTTCTATTCATTCGTGATTATGATCAc                                                   >  2:2681788/1‑98 (MQ=255)
                                                gCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGa                                                 >  2:508669/1‑100 (MQ=255)
                                                gCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGa                                                 <  1:1857948/100‑1 (MQ=255)
                                                 cGTTATCATCAGGTAATCCGTTTGCGGATAACAAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGAt                                                >  1:484995/1‑100 (MQ=255)
                                                   ttATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCTTGATCACGaatg                                              <  1:2381274/100‑4 (MQ=255)
                                                   ttATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGc                                              <  2:1539370/100‑1 (MQ=255)
                                                   ttATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGc                                              <  2:3645225/100‑1 (MQ=255)
                                                     atcatcGGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTg                                            >  1:3114270/1‑100 (MQ=255)
                                                      tcatcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGt                                           <  1:1213445/100‑1 (MQ=255)
                                                        atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                         >  2:1046681/1‑100 (MQ=255)
                                                        atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                         >  2:104573/1‑100 (MQ=255)
                                                         tcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                         <  1:226669/99‑1 (MQ=255)
                                                         tcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGc                                        <  2:3391516/100‑1 (MQ=255)
                                                           aGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTTGTGATCAAGATCACGATGCTGTc                                          >  2:2365877/1‑96 (MQ=255)
                                                           aGGTAATCCGTTTGCGGATAACCAAAAATGCAGCTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACTATGTTGTCGCa                                       <  2:2188613/99‑1 (MQ=255)
                                                            ggTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAAt                                     >  1:3676969/1‑100 (MQ=255)
                                                            ggTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAAt                                     >  1:2779065/1‑100 (MQ=255)
                                                             gTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATg                                    <  1:2368313/100‑1 (MQ=255)
                                                              tattctgTTTGTGGATAATCAAAAATGTAGTTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTTGTGATCAAGATCACGATGCTGTCGTAATGt                                   <  2:1421586/94‑1 (MQ=255)
                                                              tAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGt                                   <  2:1425044/100‑1 (MQ=255)
                                                              tAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑TTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTa                                  >  1:548674/1‑100 (MQ=255)
                                                               aaTCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTa                                  <  1:55771/100‑1 (MQ=255)
                                                                aTCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaa                                 >  2:2002096/1‑100 (MQ=255)
                                                                 tacGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑TTTTT‑TTGAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaa                                >  1:2810347/3‑100 (MQ=255)
                                                                 tCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaa                                >  1:4070746/1‑100 (MQ=255)
                                                                  ccGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaaaa                               >  1:2625840/1‑100 (MQ=255)
                                                                      ttGCGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAt                           <  2:1403608/100‑1 (MQ=255)
                                                                         cGGATAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGc                        >  2:3249128/1‑100 (MQ=255)
                                                                         cGGATAACCAAAAATGCAAGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTg                         >  2:854442/1‑99 (MQ=255)
                                                                            aTAACCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAg                     <  1:1472302/100‑1 (MQ=255)
                                                                             tAACCAAAAAGGCATGTTAATTTTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGa                    <  1:1385247/100‑1 (MQ=255)
                                                                              aaCCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGaa                   >  1:2763997/1‑100 (MQ=255)
                                                                              aaCCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGAGCACGATGCTGTCGCAATGTAAAAACATTGCAAGaa                   <  2:3820673/100‑1 (MQ=255)
                                                                               aCCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                  >  1:2716824/1‑100 (MQ=255)
                                                                               aCCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                  >  1:121682/1‑100 (MQ=255)
                                                                               aCCAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTCTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                  >  1:3819544/1‑100 (MQ=255)
                                                                                 caaatgtGCAGGTTAATTG‑TTTTTGTG‑AATAGCGATAGTTTCGATTCATTCGTTATCAAGATCACGATGCTGTCGCAATGTAAAAGCATTGCaa                      <  2:3880459/88‑1 (MQ=255)
                                                                                 cAAAAATGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCa                <  2:1083107/100‑1 (MQ=255)
                                                                                    aaaTGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGt             <  1:1138278/100‑1 (MQ=255)
                                                                                    aaaTGCAGGTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGt             <  1:777700/100‑1 (MQ=255)
                                                                                        gCAGGTTAATTG‑TTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTaaaa        <  1:2983884/100‑1 (MQ=255)
                                                                                           ggTTAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATc      >  1:2650393/1‑100 (MQ=255)
                                                                                             ttAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCa     <  2:880790/99‑1 (MQ=255)
                                                                                              tAATTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAtt   <  1:1967310/100‑1 (MQ=255)
                                                                                              tAATTG‑TTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAttt  <  1:3987073/100‑1 (MQ=255)
                                                                                               aaTTGTTTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAttt  <  2:442897/100‑1 (MQ=255)
                                                                                                    |                                                                                                
TTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑TTTTT‑TT‑AATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTT  >  NC_000913/1712200‑1712393

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