Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,712,518 (C)6→5 100% intergenic (+360/‑251) 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,5130C.92.5% 382.7 / 18.4 93intergenic (+355/‑256)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base . (46/40);  total (51/42)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

TTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTA  >  NC_000913/1712416‑1712612
                                                                                                 |                                                                                                   
ttGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccc                                                                                                   <  1:1419217/100‑1 (MQ=255)
 tGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAcccc                                                                                                  <  2:2412421/100‑1 (MQ=255)
    ccTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccctc                                                                                               >  2:1904771/1‑98 (MQ=255)
     cTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccctct                                                                                              >  2:3017148/1‑97 (MQ=255)
      tACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCAccccctctc                                                                                             <  2:238306/100‑5 (MQ=255)
        cAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTcag                                                                                          <  2:403691/100‑1 (MQ=255)
        cAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTcag                                                                                          <  2:3132372/100‑1 (MQ=255)
        cAAAACATTACACTGGCTATTTTTCAGAAACTGGACTACCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTcag                                                                                          <  1:1387787/100‑1 (MQ=255)
         aaaaCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTcagt                                                                                         >  1:905639/1‑100 (MQ=255)
          aaaCATTACACTGGCTATTTTTCAGAAACGGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTcagtg                                                                                        <  1:2011661/100‑1 (MQ=255)
           aaCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTcagtgc                                                                                       <  1:1845020/100‑1 (MQ=255)
            aCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTcagtgca                                                                                      <  2:761806/100‑1 (MQ=255)
             cATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTcagtgcag                                                                                     <  2:2893474/100‑1 (MQ=255)
             cATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTcagtgcag                                                                                     <  2:1675172/100‑1 (MQ=255)
              aTTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTcagtgcagt                                                                                    <  1:2330392/100‑1 (MQ=255)
              aTTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTcagtgcagt                                                                                    >  1:2614807/1‑100 (MQ=255)
              aTTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTcagtgcagt                                                                                    <  1:1063664/100‑1 (MQ=255)
                  cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGaaa                                                                                >  1:2322129/1‑100 (MQ=255)
                  cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGaaa                                                                                >  1:2646579/1‑100 (MQ=255)
                  cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGaaa                                                                                >  2:2188445/1‑100 (MQ=255)
                  cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGaaa                                                                                >  2:2486742/1‑100 (MQ=255)
                      ggCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAAt                                                                            >  1:3069700/1‑100 (MQ=255)
                         tATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTg                                                                         <  2:1251860/100‑1 (MQ=255)
                          aTTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGc                                                                        >  2:116247/1‑100 (MQ=255)
                           tttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGcc                                                                       <  2:2266045/100‑1 (MQ=255)
                            ttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCg                                                                      >  2:1736894/1‑100 (MQ=255)
                            ttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCg                                                                      >  1:130802/1‑100 (MQ=255)
                                cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTAc                                                                  >  2:1825700/1‑100 (MQ=255)
                                cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTAc                                                                  >  2:2550438/1‑100 (MQ=255)
                                cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTAc                                                                  <  2:2569542/100‑1 (MQ=255)
                                cAGAAACCGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTGAAGCCATCACTCCCTCTCAGTGCAGTGAAAAAATCTGCCg                                                                      <  1:2131206/96‑1 (MQ=255)
                                 aGAAACTGTACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCACAGTGCAGTGAAAAAATCTGCCGTTAc                                                                  >  1:2195998/1‑99 (MQ=255)
                                 aGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACg                                                                 >  1:24371/1‑100 (MQ=255)
                                      cTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGttttt                                                            <  2:2260711/100‑1 (MQ=255)
                                      cTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGttttt                                                            >  1:2657996/1‑100 (MQ=255)
                                      cTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGttttt                                                            >  1:929833/1‑100 (MQ=255)
                                        ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttttt                                                           >  1:2395029/1‑100 (MQ=255)
                                        ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTg                                                          <  2:3269008/100‑1 (MQ=255)
                                        ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTg                                                          >  2:1484284/1‑100 (MQ=255)
                                        ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTg                                                          >  2:49186/1‑100 (MQ=255)
                                        ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTg                                                          <  1:2188430/100‑1 (MQ=255)
                                        ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCGCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGttttttt                                                          >  1:2163006/1‑99 (MQ=255)
                                         gACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTg                                                          <  1:2664595/99‑1 (MQ=255)
                                         gACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGa                                                         >  1:286830/1‑100 (MQ=255)
                                          aCTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                        <  1:834001/100‑1 (MQ=255)
                                             aTCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaaaaa                                                     >  1:3374389/1‑100 (MQ=255)
                                             aTCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaaaaa                                                     >  2:826992/1‑100 (MQ=255)
                                              tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAt                                                    >  2:2608362/1‑100 (MQ=255)
                                                tcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAttt                                                  <  2:929836/100‑1 (MQ=255)
                                                 cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAttt                                                  >  2:642442/1‑100 (MQ=255)
                                                 cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                                 >  2:2844675/1‑100 (MQ=255)
                                                 cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                                 >  2:635779/1‑100 (MQ=255)
                                                 cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                                 >  1:1440455/1‑100 (MQ=255)
                                                 cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                                 >  1:1687843/1‑100 (MQ=255)
                                                     aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGc                                             >  1:2253736/1‑100 (MQ=255)
                                                      gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCt                                            >  1:2294293/1‑100 (MQ=255)
                                                       tAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCt                                            <  1:1131320/99‑1 (MQ=255)
                                                       tAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                           <  2:33232/100‑1 (MQ=255)
                                                       tAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                           >  2:1977831/1‑100 (MQ=255)
                                                       tAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                           >  1:809501/1‑100 (MQ=255)
                                                        aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                           <  1:564590/100‑1 (MQ=255)
                                                        aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTGTTGAAAAATTTAACGCTgg                                          >  1:3225630/1‑100 (MQ=255)
                                                         aaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                         <  2:3194358/100‑1 (MQ=255)
                                                           aaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATa                                       <  1:1825689/100‑1 (MQ=255)
                                                            aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                      >  1:285304/1‑100 (MQ=255)
                                                             cGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAAc                                     >  2:160270/1‑100 (MQ=255)
                                                              gCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACa                                    <  1:161365/100‑1 (MQ=255)
                                                              gCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACa                                    <  2:618098/100‑1 (MQ=255)
                                                                    tAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTccc                              >  2:1886689/1‑100 (MQ=255)
                                                                    tAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTccc                              <  2:2034899/100‑1 (MQ=255)
                                                                    tAGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCg                             >  2:1634411/1‑100 (MQ=255)
                                                                     aGCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCg                             <  1:1484278/100‑1 (MQ=255)
                                                                      gCAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg                            <  1:1551804/100‑1 (MQ=255)
                                                                       cAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGa                           <  1:1176891/100‑1 (MQ=255)
                                                                       cAGAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGa                           <  1:139123/100‑1 (MQ=255)
                                                                        aGAATATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                          >  2:3279135/1‑100 (MQ=255)
                                                                         gAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAAt                         >  2:1996172/1‑100 (MQ=255)
                                                                         gAACATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAAt                         >  1:1881514/1‑100 (MQ=255)
                                                                          aaCATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATg                        <  2:1687852/100‑1 (MQ=255)
                                                                          aaCATATCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATg                        <  1:1737988/100‑1 (MQ=255)
                                                                              tatCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                    >  2:232853/1‑100 (MQ=255)
                                                                              tatCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                    <  2:809504/100‑1 (MQ=255)
                                                                               atCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGa                   >  1:775798/1‑100 (MQ=255)
                                                                               atCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGa                   >  2:2542799/1‑100 (MQ=255)
                                                                                tCGCCGTTAAGCCATCACCCCCTCTCTGTGCAGTGAAAATATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGaa                  <  2:324478/100‑1 (MQ=255)
                                                                                tCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTTTTTTTTTTTACACTAGATAACATTTCCCGGAATGGTTGaa                  >  1:3324159/1‑100 (MQ=255)
                                                                                tCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGaa                  >  2:1415124/1‑100 (MQ=255)
                                                                                tCGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGaa                  <  1:1736883/100‑1 (MQ=255)
                                                                                tCGCCGTTAAGCCATCACCCCCTCTCAGGGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGaa                  <  2:3029801/100‑1 (MQ=255)
                                                                                tCGCCGTTAAGCCATAACCCCCTCTCACGGCAGTGAAAAAATCTGCAGTTACGATTTCTGACAAATTTAACGATGGATAACATTTCACGGTAt                         >  1:1025902/1‑93 (MQ=255)
                                                                                 cGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAt                 >  2:3082596/1‑100 (MQ=255)
                                                                                 cGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAt                 <  2:126559/100‑1 (MQ=255)
                                                                                 cGCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAt                 <  2:2457014/100‑1 (MQ=255)
                                                                                  gCCGTTAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAt                 <  2:2996915/99‑1 (MQ=255)
                                                                                      ttAAGCCATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTcccc            >  1:1074266/1‑100 (MQ=255)
                                                                                      ttAAGCCATCTCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTcccc            <  1:2325583/100‑1 (MQ=255)
                                                                                           ccATCACCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTc       >  2:2312957/1‑100 (MQ=255)
                                                                                                aCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt   >  1:1512783/1‑100 (MQ=255)
                                                                                                 ccccccTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTta  <  1:642439/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
TTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTA  >  NC_000913/1712416‑1712612

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: