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)
        cAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcag                                                                                          <  2:403691/100‑1 (MQ=255)
        cAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcag                                                                                          <  2:3132372/100‑1 (MQ=255)
        cAAAACATTACACTGGCTATTTTTCAGAAACTGGACTACCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcag                                                                                          <  1:1387787/100‑1 (MQ=255)
         aaaaCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagt                                                                                         >  1:905639/1‑100 (MQ=255)
          aaaCATTACACTGGCTATTTTTCAGAAACGGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtg                                                                                        <  1:2011661/100‑1 (MQ=255)
           aaCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtgc                                                                                       <  1:1845020/100‑1 (MQ=255)
            aCATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtgca                                                                                      <  2:761806/100‑1 (MQ=255)
             cATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtgcag                                                                                     <  2:2893474/100‑1 (MQ=255)
             cATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtgcag                                                                                     <  2:1675172/100‑1 (MQ=255)
              aTTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtgcagt                                                                                    <  1:2330392/100‑1 (MQ=255)
              aTTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtgcagt                                                                                    >  1:2614807/1‑100 (MQ=255)
              aTTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTcagtgcagt                                                                                    <  1:1063664/100‑1 (MQ=255)
                  cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGaaa                                                                                >  1:2322129/1‑100 (MQ=255)
                  cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGaaa                                                                                >  1:2646579/1‑100 (MQ=255)
                  cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGaaa                                                                                >  2:2188445/1‑100 (MQ=255)
                  cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGaaa                                                                                >  2:2486742/1‑100 (MQ=255)
                      ggCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAAt                                                                            >  1:3069700/1‑100 (MQ=255)
                         tATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTg                                                                         <  2:1251860/100‑1 (MQ=255)
                          aTTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGc                                                                        >  2:116247/1‑100 (MQ=255)
                           tttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGcc                                                                       <  2:2266045/100‑1 (MQ=255)
                            ttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCg                                                                      >  2:1736894/1‑100 (MQ=255)
                            ttttCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCg                                                                      >  1:130802/1‑100 (MQ=255)
                                cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTAc                                                                  >  2:1825700/1‑100 (MQ=255)
                                cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTAc                                                                  >  2:2550438/1‑100 (MQ=255)
                                cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTAc                                                                  <  2:2569542/100‑1 (MQ=255)
                                cAGAAACCGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTGAAGCCATCA‑CTCCCTCTCAGTGCAGTGAAAAAATCTGCCg                                                                      <  1:2131206/96‑1 (MQ=255)
                                 aGAAACTGTACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCACAGTGCAGTGAAAAAATCTGCCGTTAc                                                                  >  1:2195998/1‑99 (MQ=255)
                                 aGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACg                                                                 >  1:24371/1‑100 (MQ=255)
                                      cTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGttttt                                                            <  2:2260711/100‑1 (MQ=255)
                                      cTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGttttt                                                            >  1:2657996/1‑100 (MQ=255)
                                      cTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGttttt                                                            >  1:929833/1‑100 (MQ=255)
                                        ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttttt                                                           >  1:2395029/1‑100 (MQ=255)
                                        ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTg                                                          <  2:3269008/100‑1 (MQ=255)
                                        ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTg                                                          >  2:1484284/1‑100 (MQ=255)
                                        ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTg                                                          >  2:49186/1‑100 (MQ=255)
                                        ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTg                                                          <  1:2188430/100‑1 (MQ=255)
                                        ggACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATC‑GCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGttttttt                                                          >  1:2163006/1‑99 (MQ=255)
                                         gACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTg                                                          <  1:2664595/99‑1 (MQ=255)
                                         gACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGa                                                         >  1:286830/1‑100 (MQ=255)
                                          aCTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                        <  1:834001/100‑1 (MQ=255)
                                             aTCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaaaaa                                                     >  1:3374389/1‑100 (MQ=255)
                                             aTCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaaaaa                                                     >  2:826992/1‑100 (MQ=255)
                                              tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAt                                                    >  2:2608362/1‑100 (MQ=255)
                                                tcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAttt                                                  <  2:929836/100‑1 (MQ=255)
                                                 cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAAttt                                                  >  2:642442/1‑100 (MQ=255)
                                                 cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                                 >  2:2844675/1‑100 (MQ=255)
                                                 cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                                 >  2:635779/1‑100 (MQ=255)
                                                 cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                                 >  1:1440455/1‑100 (MQ=255)
                                                 cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                                 >  1:1687843/1‑100 (MQ=255)
                                                     aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGc                                             >  1:2253736/1‑100 (MQ=255)
                                                      gTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCt                                            >  1:2294293/1‑100 (MQ=255)
                                                       tAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCt                                            <  1:1131320/99‑1 (MQ=255)
                                                       tAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                           <  2:33232/100‑1 (MQ=255)
                                                       tAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                           >  2:1977831/1‑100 (MQ=255)
                                                       tAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                           >  1:809501/1‑100 (MQ=255)
                                                        aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                           <  1:564590/100‑1 (MQ=255)
                                                        aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTGTTGAAAAATTTAACGCTgg                                          >  1:3225630/1‑100 (MQ=255)
                                                         aaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGa                                         <  2:3194358/100‑1 (MQ=255)
                                                           aaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATa                                       <  1:1825689/100‑1 (MQ=255)
                                                            aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                      >  1:285304/1‑100 (MQ=255)
                                                             cGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAAc                                     >  2:160270/1‑100 (MQ=255)
                                                              gCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACa                                    <  1:161365/100‑1 (MQ=255)
                                                              gCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACa                                    <  2:618098/100‑1 (MQ=255)
                                                                    tAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTccc                              >  2:1886689/1‑100 (MQ=255)
                                                                    tAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTccc                              <  2:2034899/100‑1 (MQ=255)
                                                                    tAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTG‑AAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCg                             >  2:1634411/1‑100 (MQ=255)
                                                                     aGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCg                             <  1:1484278/100‑1 (MQ=255)
                                                                      gCAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCgg                            <  1:1551804/100‑1 (MQ=255)
                                                                       cAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGa                           <  1:1176891/100‑1 (MQ=255)
                                                                       cAGAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGa                           <  1:139123/100‑1 (MQ=255)
                                                                        aGAATATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                          >  2:3279135/1‑100 (MQ=255)
                                                                         gAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAAt                         >  2:1996172/1‑100 (MQ=255)
                                                                         gAACATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAAt                         >  1:1881514/1‑100 (MQ=255)
                                                                          aaCATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATg                        <  2:1687852/100‑1 (MQ=255)
                                                                          aaCATATCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATg                        <  1:1737988/100‑1 (MQ=255)
                                                                              tatCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                    >  2:232853/1‑100 (MQ=255)
                                                                              tatCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                    <  2:809504/100‑1 (MQ=255)
                                                                               atCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGa                   >  1:775798/1‑100 (MQ=255)
                                                                               atCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGa                   >  2:2542799/1‑100 (MQ=255)
                                                                                tCGCCGTTAAGCCATCA‑CCCCCTCTCTGTGCAGTGAAAATATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGaa                  <  2:324478/100‑1 (MQ=255)
                                                                                tCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTTTTTTTTTTTACACTAGATAACATTTCCCGGAATGGTTGaa                  >  1:3324159/1‑100 (MQ=255)
                                                                                tCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGaa                  >  2:1415124/1‑100 (MQ=255)
                                                                                tCGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGaa                  <  1:1736883/100‑1 (MQ=255)
                                                                                tCGCCGTTAAGCCATCA‑CCCCCTCTCAGGGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGaa                  <  2:3029801/100‑1 (MQ=255)
                                                                                tCGCCGTTAAGCCAT‑AACCCCCTCTCACGGCAGTGAAAAAATCTGCAGTTACGATTTCTGACAAATTTAACGATGGATAACATTTCACGGTAt                         >  1:1025902/1‑93 (MQ=255)
                                                                                 cGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAt                 >  2:3082596/1‑100 (MQ=255)
                                                                                 cGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAt                 <  2:126559/100‑1 (MQ=255)
                                                                                 cGCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAt                 <  2:2457014/100‑1 (MQ=255)
                                                                                  gCCGTTAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAt                 <  2:2996915/99‑1 (MQ=255)
                                                                                      ttAAGCCATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTcccc            >  1:1074266/1‑100 (MQ=255)
                                                                                      ttAAGCCATC‑TCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTcccc            <  1:2325583/100‑1 (MQ=255)
                                                                                           ccATCA‑CCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTc       >  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: ‑