Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,590 C→A 7.4% intergenic (+72/‑33) yqhD → / → dkgA NADPH‑dependent aldehyde reductase YqhD/methylglyoxal reductase DkgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,156,5900CA 7.4% 156.9 / 4.5 62intergenic (+72/‑33)yqhD/dkgANADPH‑dependent aldehyde reductase YqhD/methylglyoxal reductase DkgA
Reads supporting (aligned to +/- strand):  ref base C (33/24);  new base A (2/3);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

CCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAA  >  NC_000913/3156452‑3156729
                                                                                                                                          |                                                                                                                                           
ccAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGcc                                                                                                                                         >  2:84698/1‑142 (MQ=255)
ccAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTCTCGTTTGCGGGCCTTTCGTCCAGACTGAAGATCACAACACCTCACCGGAGCCAGca                                                                                                                                         >  2:345443/1‑142 (MQ=255)
   agtGGGCGAAATTCATGACATTATTTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTGTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCg                                                                                                                                      <  1:829024/141‑1 (MQ=255)
   aCTTGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCg                                                                                                                                      <  1:22661/143‑1 (MQ=255)
     tGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACAACTCACCGGAGCAGGCTCCGGt                                                                                                                                    <  1:649771/143‑1 (MQ=255)
       ggCGACAATCATGCCGTTACGTCGGATGTCCGCCGCCGTATACTCTCTGTCGCCCGCTAAGCTTTTTACGCCTCAAACCTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGa                                                                                                                                  <  2:1021028/143‑1 (MQ=255)
        gCGAAAATCATGACATTACGGTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTGCGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                                                                 <  2:855597/143‑1 (MQ=255)
            aactCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAggtgt                                                                                                                             <  1:483086/140‑5 (MQ=255)
               tCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGCGCTCatac                                                                                                                          >  2:384461/1‑142 (MQ=255)
                cATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGGGAGTTCAtata                                                                                                                         >  1:951109/1‑143 (MQ=255)
                  tGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaaa                                                                                                                       <  2:532380/143‑1 (MQ=255)
                  tGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaaa                                                                                                                       <  1:232730/143‑1 (MQ=255)
                     cATTACGTTGTATGTCAGCCTCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCTTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagga                                                                                                                    <  1:884630/143‑1 (MQ=255)
                     cATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagga                                                                                                                    <  1:1003775/143‑1 (MQ=255)
                      aTTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                                                                   >  1:519900/1‑143 (MQ=255)
                           gTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACg                                                                                                              >  1:587038/1‑143 (MQ=255)
                                atttCAGCCGCCGTACATACGAGGCCGCCCGTTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCAGTTCGTCCAGAGTTAACTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATgg                                                                                                         <  2:451975/140‑1 (MQ=255)
                                   tCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTa                                                                                                      >  1:890121/1‑143 (MQ=255)
                                    cAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTaa                                                                                                     <  2:890829/143‑1 (MQ=255)
                                       ccgccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGTCATTTCGTCCGGACCTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGGGAGGTCATATAAAGGAGGAACGTAGGGCTGATca                                                                                                  >  2:778267/1‑142 (MQ=255)
                                               tataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtta                                                                                          >  1:998351/1‑143 (MQ=255)
                                                 taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCGAATCCAACCGttatt                                                                                        >  1:225480/1‑143 (MQ=255)
                                                  aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATACAGGAGGAACGTATGGCTAATCCGACCGttatta                                                                                       >  1:779863/1‑143 (MQ=255)
                                                   cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTaa                                                                                      >  1:1017708/1‑143 (MQ=255)
                                                        ccgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCCTTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATCAAGGAGGAACGTATGGCTAATCCAACCGTTATTCAGCtgc                                                                                 >  2:589968/1‑141 (MQ=255)
                                                         cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACa                                                                                >  1:720825/1‑143 (MQ=255)
                                                          gcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAg                                                                               >  1:13711/1‑143 (MQ=255)
                                                              gcTAAGCTTTTTACGCCTCAAACTTTCGTTGTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATg                                                                           <  1:179112/143‑1 (MQ=255)
                                                              gcTAAGCTCTTTCCGCCTCAAACTCTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGAGAGTTCATATAACGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATg                                                                           <  1:983424/143‑1 (MQ=255)
                                                                tAAGCTTTTTACGCCTCAAACTTTCGGTTTCGGGCATTTCGTCCAGACTTAACTTCACAACACCTCACCAGAGCCTGCTCCGGTGAGTTCATATAAACGAGGAACGTATGGCTAATACACACGTACTTGAGCTACAGGaaagc                                                                         >  2:42108/1‑139 (MQ=255)
                                                                   gCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCCTATAAAGCAGGCACGTATGGGTAATCCACCCGCTATTAGACTGGAGGATGGCGAt                                                                      >  2:501785/1‑143 (MQ=255)
                                                                   gCTTTTGACGCCTCAAACTTTCGTTTTCGGGCATGTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATATAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAAt                                                                      <  2:548183/143‑1 (MQ=255)
                                                                       tttACGCCTCAAGCTTTGGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCa                                                                  <  2:711410/143‑1 (MQ=255)
                                                                         tACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                                                                <  1:647802/143‑1 (MQ=255)
                                                                              cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGCGAGTTCATATAAAGGCGGCACGTACGGCTCATCCAACCGTTCTTAAGCCCCAGGCTCGCACTGCCATGCCCCa                                                           >  2:656855/1‑143 (MQ=255)
                                                                              cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACATCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTGCAGGATGGCAATGTCATGCcaca                                                           >  2:51306/1‑143 (MQ=255)
                                                                               tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCACTGTCATGCCCCAg                                                          >  1:183356/1‑143 (MQ=255)
                                                                                ccaaCTCTCGTTTTCGGGCATTTCGCCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGc                                                         <  1:781243/141‑1 (MQ=255)
                                                                                cAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGc                                                         >  1:1023272/1‑143 (MQ=255)
                                                                                     tttCGTGTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa                                                    <  1:408401/143‑1 (MQ=255)
                                                                                      ttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc                                                   >  2:280454/1‑143 (MQ=255)
                                                                                              cGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGAGACTGGGCggc                                           >  1:795255/1‑141 (MQ=255)
                                                                                                    tttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaa                                     <  1:540188/143‑1 (MQ=255)
                                                                                                      tCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagc                                   <  1:946826/143‑1 (MQ=255)
                                                                                                          ccAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGt                               >  2:375961/1‑143 (MQ=255)
                                                                                                               cgtAAGTTCACACCACCTCACCGGGGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCGCCAGCTGGGACTGGGCGGCTGGCAAGCAAGTAATGa                          <  2:873532/141‑1 (MQ=255)
                                                                                                                ttAAGTTCACAACACCTCCGCTGAGCCTGCCCCGATCAGATCAGAGAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAg                         <  2:741380/143‑1 (MQ=255)
                                                                                                                ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAg                         >  1:638820/1‑143 (MQ=255)
                                                                                                                    gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAg                     <  1:84698/143‑1 (MQ=255)
                                                                                                                          aacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTACTCCAACCGTTATTAAGCTACAGGATGGCAATGTCGTGCCCCAGCTGGGACTGGGCGTCTGGCAAGGAAGTAATGAGGAAGTAATCa               >  1:380314/1‑143 (MQ=255)
                                                                                                                          aacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCa               >  2:784199/1‑143 (MQ=255)
                                                                                                                                tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCa         >  1:1429/1‑143 (MQ=255)
                                                                                                                                tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGGAATCCACGccc         >  2:804946/1‑142 (MQ=255)
                                                                                                                                 cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGCACTGGGCGTCTGGCAAGCAAGTCATGAGGAAGGAATCACCGCCAt        >  2:501962/1‑143 (MQ=255)
                                                                                                                                  aCCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTCCTGGATTGCAATGTCATGCCCCAGCTGGGACGGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCAtt       <  1:345443/143‑1 (MQ=255)
                                                                                                                                  aCCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCAtt       >  1:451957/1‑143 (MQ=255)
                                                                                                                                  aCCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCAtt       >  1:110055/1‑143 (MQ=255)
                                                                                                                                  aCCGGAGCAGGCTCCGGTGAGTTCCTATAAAGGAAGAACGTATGGCTAATCCAACCGTTCTTAAGCTCCCGCATGGCAATGCCATGCCCCAGCTGGGACCGGGCGTCTGGCAAGCCAGTAATGAGGCAGTAATCACCGCCagt       >  2:198669/1‑141 (MQ=255)
                                                                                                                                   ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTAAATGATTGCAATGTCATGCCCCAGCTGGGACTAGGCGTCTGACAAGCACGTAATGAGGAAGGAAACACCGCCCtta      >  2:264474/1‑142 (MQ=255)
                                                                                                                                      gAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaa   <  2:779863/143‑1 (MQ=255)
                                                                                                                                       aGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCaaaa  <  1:784199/143‑1 (MQ=255)
                                                                                                                                       aGCCTGCTCCGGAGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCCAGTAATGAGGAAGTAATCACCACCATTCaaac  >  2:534408/1‑142 (MQ=255)
                                                                                                                                          |                                                                                                                                           
CCAACTGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAA  >  NC_000913/3156452‑3156729

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