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.T96.3% 268.0 / 3.0 80intergenic (+141/‑470)nth/dtpADNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base T (47/30);  total (47/33)
Fisher's exact test for biased strand distribution p-value = 6.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑TTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTT  >  NC_000913/1712200‑1712393
                                                                                                    |                                                                                                
ttATCCTCTTTTAAAACATCTTTTTAACGTCAGTGATGCGATGGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTg                                                                                                   <  1:2374203/100‑1 (MQ=255)
ttATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTg                                                                                                   >  1:2795366/1‑100 (MQ=255)
  aTCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑tt                                                                                               >  2:1737999/1‑100 (MQ=255)
   tCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑ttt                                                                                             >  1:2311410/1‑100 (MQ=255)
    ccTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑tttt                                                                                            >  2:1845031/1‑100 (MQ=255)
     ctctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑ttttt                                                                                           >  1:1693428/1‑100 (MQ=255)
       ctTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑ttttttt                                                                                         >  2:2923689/1‑100 (MQ=255)
        ttttAAACCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑tttttttt                                                                                        <  2:1906699/100‑2 (MQ=255)
        ttttAAAACATTTTTTTAACGTCAATGATGTCATTGCTTAGTGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑tttttttt                                                                                        <  1:1024041/100‑2 (MQ=255)
        ttttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑tttttttt                                                                                        <  1:2855046/100‑2 (MQ=255)
        ttttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑tttttttt                                                                                        <  2:2017164/100‑2 (MQ=255)
         tttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑tttttttta                                                                                       <  2:273990/100‑3 (MQ=255)
         tttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑tttttttta                                                                                       <  2:85061/100‑3 (MQ=255)
          ttAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑ttttttttaa                                                                                      <  2:2777282/100‑4 (MQ=255)
           tAAAACATCGTTTTAACGTCAATGATGCCATTGCTAAGCGTTATCATCAGTTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑ttttttttaat                                                                                     <  1:1532382/100‑5 (MQ=255)
              aaCATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGc                                                                                   <  2:717209/100‑1 (MQ=255)
                cATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGa                                                                                 >  2:195505/1‑100 (MQ=255)
                    tttttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAAtt                                                                             >  2:139123/1‑100 (MQ=255)
                    tttttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAAtt                                                                             >  2:2604990/1‑100 (MQ=255)
                     ttttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAAttt                                                                            >  1:3294623/1‑100 (MQ=255)
                       ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTc                                                                          <  2:1807041/100‑1 (MQ=255)
                       ttAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATATTc                                                                          >  1:1945028/1‑100 (MQ=255)
                        tAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCt                                                                         <  2:1707490/100‑1 (MQ=255)
                         aaCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTa                                                                        <  2:1645794/100‑1 (MQ=255)
                         aaCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTa                                                                        <  2:1358707/100‑1 (MQ=255)
                          aCGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTat                                                                       >  1:3268980/1‑100 (MQ=255)
                          aCGTCAATGATGCCATTGCTTAGCGTTATACTCAGGTAATCCGTTTGCTGATAACCAAAAATGCATGTTAATTGTTTTTTTTAATAGCGAAATTTTCTa                                                                        <  1:2207066/99‑1 (MQ=255)
                           cgatgatgatgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCGAAAATGCAGGTGAATTG‑‑TTGTTTTTAATAGCGAACTTTTCt                                                                         <  1:2028889/92‑1 (MQ=255)
                           cGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTATTATCGAAATTTTCTat                                                                       >  1:1267295/1‑99 (MQ=255)
                           cGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTatt                                                                      >  1:251750/1‑100 (MQ=255)
                           cGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTatt                                                                      >  2:225983/1‑100 (MQ=255)
                             tCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattca                                                                    >  2:2664617/1‑100 (MQ=255)
                             tCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattca                                                                    <  2:2794787/100‑1 (MQ=255)
                             tCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattca                                                                    >  1:485578/1‑100 (MQ=255)
                              cAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTcagctgtt                                                                   <  2:191815/100‑8 (MQ=255)
                              cAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTattcat                                                                   >  2:1419223/1‑100 (MQ=255)
                                   atgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGt                                                               <  2:1873033/99‑1 (MQ=255)
                                    tgCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGa                                                             >  1:2105337/1‑100 (MQ=255)
                                     gCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGAt                                                            <  2:379380/100‑1 (MQ=255)
                                     gCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGAt                                                            <  2:942627/100‑1 (MQ=255)
                                      ccATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATc                                                           >  2:816385/1‑100 (MQ=255)
                                       cATTGCTTAGCGTTATCATCAGGTAATTCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCa                                                          <  1:1626439/100‑1 (MQ=255)
                                       cATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCa                                                          >  2:1131323/1‑100 (MQ=255)
                                       cATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCa                                                          >  2:2814898/1‑100 (MQ=255)
                                         ttGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAg                                                        <  1:2820609/100‑1 (MQ=255)
                                         ttGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAg                                                        >  1:3214714/1‑100 (MQ=255)
                                           gCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGAt                                                      >  2:2398984/1‑100 (MQ=255)
                                            cTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATGCGTGATCAAGATc                                                     >  2:1969883/1‑100 (MQ=255)
                                             ttAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCa                                                    <  1:2461492/100‑1 (MQ=255)
                                             tctgCGTTATCATCAGGTAATCCGTTTACAAAAAACAAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCa                                                    <  1:813429/97‑1 (MQ=255)
                                              tAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCAc                                                   <  2:655168/100‑1 (MQ=255)
                                                gCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGAtcaagatcaaga                                                 >  2:813432/1‑100 (MQ=255)
                                                gCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGa                                                 >  1:2266029/1‑100 (MQ=255)
                                                 cGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCAcca                                                 <  2:2297746/99‑3 (MQ=255)
                                                 cGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑TTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATg                                               <  1:91422/100‑1 (MQ=255)
                                                  gTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATg                                               >  1:238306/1‑100 (MQ=255)
                                                   ttATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGc                                              <  2:708009/100‑1 (MQ=255)
                                                   ttATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGc                                              <  2:822335/100‑1 (MQ=255)
                                                       catcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATATAGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCt                                             >  1:800373/1‑97 (MQ=255)
                                                       catcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATTCTGTc                                          >  1:2996580/1‑100 (MQ=255)
                                                        atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                         <  1:195505/100‑1 (MQ=255)
                                                        atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                         >  2:3234016/1‑100 (MQ=255)
                                                        atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                         >  2:2131221/1‑100 (MQ=255)
                                                        atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCg                                         <  2:2203668/100‑1 (MQ=255)
                                                        atcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGAAGCTGTc                                          >  1:1442342/1‑99 (MQ=255)
                                                        atcaGGTAATCCGTTTGCGGATAACCAAAAATGAAGGATAATTGTTTTTTTTTAATAGCGCAATTTTATAATAATACGTGTTCCAGATCACGATACTGt                                           >  1:2606059/1‑99 (MQ=255)
                                                         tcaGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTCTTCATTCGTGATCAAGATCACGATGCTGTCGc                                        <  2:3214742/100‑1 (MQ=255)
                                                           aGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCaa                                      >  2:3086448/1‑100 (MQ=255)
                                                            ggTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAAt                                     >  2:1454136/1‑100 (MQ=255)
                                                             gTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATg                                    >  1:553466/1‑100 (MQ=255)
                                                              tAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGt                                   >  1:470037/1‑100 (MQ=255)
                                                               naTCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTa                                  >  1:1562068/2‑100 (MQ=255)
                                                                aTCCGTTTGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTaa                                 <  2:1578249/100‑1 (MQ=255)
                                                                      ttGCGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACAt                           <  1:2610249/100‑1 (MQ=255)
                                                                         cGGATAATCCAAAATCCAGGTTAATTG‑‑TTAATTCTACTATACAAATTCTTTATTCATTCTAAATCAAGATCACGAGGCTGTCGCAATGTAAAAACATTg                         <  1:2826090/99‑1 (MQ=255)
                                                                         cGGATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATATCTAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGc                        >  2:311289/1‑100 (MQ=255)
                                                                          ggATAACCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCa                       >  1:405024/1‑100 (MQ=255)
                                                                              aaCCAAAAATGCAGGTTAATTGTTTTTTTTTTTTTTTTAATCTTCTATTCAGTAGTGATCAAGAGCACGATTATGTTGCAATTTAAAAACAtttaaatat                   >  1:1524902/1‑93 (MQ=255)
                                                                               aCCAAAAATGCAGGTTAATTGTTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGaa                   >  2:3341152/1‑100 (MQ=255)
                                                                               aCCAAAAATGCAGGTTAATTGTTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGaa                   >  2:1312451/1‑100 (MQ=255)
                                                                               aCCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                  >  2:105915/1‑100 (MQ=255)
                                                                               aCCAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAAt                  >  2:2440932/1‑100 (MQ=255)
                                                                                ccAAAAATGCGGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc                 >  1:3112690/1‑100 (MQ=255)
                                                                                ccAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGAACAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATc                 >  2:1048437/1‑100 (MQ=255)
                                                                                 cAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCa                >  2:2097045/1‑100 (MQ=255)
                                                                                 cAAAAATGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCa                >  2:17071/1‑100 (MQ=255)
                                                                                  aaaaaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAg               >  2:1393818/1‑100 (MQ=255)
                                                                                   aaaaTGCAGGTTAATTGTTTTTTTTAATATTGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCa                >  1:1119235/1‑98 (MQ=255)
                                                                                    aaaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGt             <  1:242895/100‑1 (MQ=255)
                                                                                    aaaTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGt             <  1:2796454/100‑1 (MQ=255)
                                                                                      aTGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTa           <  2:1231890/100‑1 (MQ=255)
                                                                                       tGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTaa          <  2:453269/100‑1 (MQ=255)
                                                                                       tGCAGGTTAATTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTaa          >  1:431868/1‑100 (MQ=255)
                                                                                              tAATTGTTTTTTTTAATAGCGAAATTGTATATTCATTCGTTATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAtt   <  1:1083474/100‑1 (MQ=255)
                                                                                               aaTTGTTTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCAttt  <  2:271664/100‑1 (MQ=255)
                                                                                                    |                                                                                                
TTATCCTCTTTTAAAACATCTTTTTAACGTCAATGATGCCATTGCTTAGCGTTATCATCAGGTAATCCGTTTGCGGATAACCAAAAATGCAGGTTAATTG‑‑TTTTTTTAATAGCGAAATTTTCTATTCATTCGTGATCAAGATCACGATGCTGTCGCAATGTAAAAACATTGCAAGAATCAGGTTAAAATCATTT  >  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: