Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,591 T→G 15.9% intergenic (+73/‑32) yqhD → / → dkgA aldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,156,5910TG15.9% 376.6 / 66.2 151intergenic (+73/‑32)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base T (57/70);  new base G (16/8);  total (73/78)
Fisher's exact test for biased strand distribution p-value = 7.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

Only 100 of 153 total aligned reads displayed.

CGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTG  >  NC_000913/3156492‑3156690
                                                                                                   |                                                                                                   
cgccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCt                                                                                                     <  2:521981/100‑1 (MQ=255)
  ccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGc                                                                                                   <  1:758466/100‑1 (MQ=255)
     tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTcn                                                                                                <  2:359312/100‑2 (MQ=255)
     tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTcc                                                                                                >  2:2933340/1‑100 (MQ=255)
     tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTcc                                                                                                <  2:2700649/100‑1 (MQ=255)
     tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTcc                                                                                                <  2:2396879/100‑1 (MQ=255)
      atataCGCAGCCGCCCGTGCAGCTTTTTACGCCTCAATCTTTCGTTATCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTcc                                                                                                <  2:2215629/99‑1 (MQ=255)
       tataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCgg                                                                                              <  2:1998413/100‑1 (MQ=255)
       tataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCgg                                                                                              <  1:928567/100‑1 (MQ=255)
       tataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCgg                                                                                              <  1:1590784/100‑1 (MQ=255)
        ataCGAAGCCGCCCTCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGt                                                                                             <  2:3452157/100‑1 (MQ=255)
        ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                             >  1:1069372/1‑100 (MQ=255)
        ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGt                                                                                             <  2:1146911/100‑1 (MQ=255)
        ataCGAAGCCGCACGCTAAGCTTTTTACGACTCAAACTTTCGTTTTCGGGCAGTTCGTCCAGACTTAAGTTAACAACACCTCACCGGAGCCTGCTCCGGt                                                                                             <  2:1894610/100‑1 (MQ=255)
         taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                            <  1:1655678/100‑1 (MQ=255)
         taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                            <  1:142971/100‑1 (MQ=255)
         taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                            <  2:1363325/100‑1 (MQ=255)
         taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                            >  2:1188231/1‑100 (MQ=255)
          aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGa                                                                                           <  1:3167806/100‑1 (MQ=255)
           cGAAGCCGCCCGCTAAGCTTTTTACGCCTCGAACTTTCGTTTTCGGACATTTCGTCCAGACTTAAGTTCACAACACCGCACCGGAGCCTGCTCCGGTGAg                                                                                          <  2:736194/100‑1 (MQ=255)
           cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTGGTTTTCGGGCATTTCCTCCAGACGTAATTTCACGACACCTCACCGGAGCAGGCTCCGGTGa                                                                                           <  1:494442/99‑1 (MQ=255)
           cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                          >  1:2152526/1‑100 (MQ=255)
           cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                          >  1:3902068/1‑100 (MQ=255)
           cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                          <  1:1512160/100‑1 (MQ=255)
            gAAGCCGCCCGCTAAGCTTTTTACGCCTCAATCTGTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAGCACCTCACCGGAGCCTGCTCCGGTGAgn                                                                                         <  1:1123097/100‑2 (MQ=255)
            gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAgg                                                                                         <  2:2842594/100‑2 (MQ=255)
             aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGtt                                                                                        <  1:384058/100‑1 (MQ=255)
               gccgccCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAggtgt                                                                                      <  2:3907726/100‑5 (MQ=39)
                 cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCata                                                                                    <  2:333447/100‑1 (MQ=255)
                   cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                  <  1:306950/100‑1 (MQ=255)
                     cgcTAATCTTTTTACGCCTCGAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACACCACCTCACCGGAGCATGCTCCGGTGAggtgttgtgaa                                                                                <  2:1373739/100‑11 (MQ=37)
                     cgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaaa                                                                                >  1:1300307/1‑100 (MQ=255)
                       cTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagg                                                                              <  1:2368430/100‑1 (MQ=255)
                        tAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagga                                                                             >  2:1109688/1‑100 (MQ=255)
                         aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                            <  2:1133142/100‑1 (MQ=255)
                         aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                            <  2:3592874/100‑1 (MQ=255)
                           gCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAggagga                                                                          >  1:3581319/1‑100 (MQ=255)
                           gCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAggagga                                                                          >  2:1540357/1‑100 (MQ=255)
                             tttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                        <  1:496030/100‑1 (MQ=255)
                                  aCGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATg                                                                   <  2:3501067/100‑1 (MQ=255)
                                  aCGCCTCAAACTTTCGTTTTCGGGCATCTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATg                                                                   <  1:391522/100‑1 (MQ=255)
                                    gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                                 <  1:1601187/100‑1 (MQ=255)
                                    gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                                 <  2:3577685/100‑1 (MQ=255)
                                     ccTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCt                                                                >  1:3908091/1‑100 (MQ=255)
                                     ccTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCt                                                                <  1:2429740/100‑1 (MQ=255)
                                     ccTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCt                                                                >  1:3423860/1‑100 (MQ=255)
                                       tCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTaa                                                              >  2:2330318/1‑100 (MQ=255)
                                         aaaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATc                                                            <  1:1418143/100‑1 (MQ=255)
                                          aaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATcc                                                           >  1:3875666/1‑100 (MQ=255)
                                          aaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATcc                                                           >  1:3624199/1‑100 (MQ=255)
                                           aCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCa                                                          >  1:2217966/1‑100 (MQ=255)
                                           aCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCa                                                          <  2:1080604/100‑1 (MQ=255)
                                           aCTTTCGGTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGATCCTGCTCCGGTTAGTTCATATAAAGGAGGAACGTATGGCTAATCCa                                                          <  2:2324265/100‑1 (MQ=255)
                                            cTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCaa                                                         >  2:3866183/1‑100 (MQ=255)
                                              ttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAcc                                                       <  1:3431884/100‑1 (MQ=255)
                                               tCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCg                                                      <  2:1778715/100‑1 (MQ=255)
                                                 gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCGCCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                    <  2:283738/100‑1 (MQ=255)
                                                 gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                    >  1:1159602/1‑100 (MQ=255)
                                                  ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtta                                                   <  1:3898024/100‑1 (MQ=255)
                                                  ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtta                                                   >  2:2504002/1‑100 (MQ=255)
                                                     tCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGTAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatta                                                >  2:3755948/1‑100 (MQ=255)
                                                          cATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTa                                           >  1:2206857/1‑100 (MQ=255)
                                                                  ccAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATgg                                   <  2:2143039/100‑1 (MQ=255)
                                                                     gACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCaa                                >  2:2983943/1‑100 (MQ=255)
                                                                     gACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCaa                                <  2:2002321/100‑1 (MQ=255)
                                                                      aCTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCaa                                <  2:1174814/99‑1 (MQ=255)
                                                                        ttAAGTTCACAACACCACACCGGAGCCTGCTCCGGTGAGATCATATAAATGAAGCACGTAAGGCGAAATAAACCGTTATTAAGCAACAAGAAGGCaaag                              >  2:3832155/1‑97 (MQ=255)
                                                                           aGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCAt                          >  2:390106/1‑100 (MQ=255)
                                                                            gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                         <  1:2095180/100‑1 (MQ=255)
                                                                            gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                         >  1:2583220/1‑100 (MQ=255)
                                                                             ttCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGc                        <  1:330392/100‑1 (MQ=255)
                                                                               cacaACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGccc                      >  2:3757929/1‑100 (MQ=255)
                                                                               cacaACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGccc                      <  1:1515651/100‑1 (MQ=255)
                                                                                 caacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAACGTTATTAAGCTACAGGATGGCAATGTCATGCCCCa                    >  2:1547760/1‑100 (MQ=255)
                                                                                     acCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGCGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTg                <  2:3425531/100‑1 (MQ=255)
                                                                                      ccTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg               >  1:168133/1‑100 (MQ=255)
                                                                                      ccTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg               <  1:3076236/100‑1 (MQ=255)
                                                                                      ccTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg               >  1:3750960/1‑100 (MQ=255)
                                                                                      acTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg               >  1:2382464/2‑100 (MQ=255)
                                                                                       cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg              >  1:2594880/1‑100 (MQ=255)
                                                                                       cTCACCGGAGCAGGCTCCGGTGAGTTCATATAAGGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg              >  1:2973296/1‑100 (MQ=255)
                                                                                       cTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg              >  1:267208/1‑100 (MQ=255)
                                                                                        tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa             <  1:3209333/100‑1 (MQ=255)
                                                                                        tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa             >  2:2896430/1‑100 (MQ=255)
                                                                                         cACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc            >  1:715592/1‑100 (MQ=255)
                                                                                         cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGCACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCCTGCCCCAGCTGGGAc            >  1:1234887/1‑100 (MQ=255)
                                                                                         cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc            >  2:3126282/1‑100 (MQ=255)
                                                                                         cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc            >  1:2971457/1‑100 (MQ=255)
                                                                                         cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc            >  1:1721544/1‑100 (MQ=255)
                                                                                         cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc            >  2:1762366/1‑100 (MQ=255)
                                                                                           ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTg          <  1:1858130/100‑1 (MQ=255)
                                                                                           ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTg          <  2:1028233/100‑1 (MQ=255)
                                                                                           ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTg          >  1:1021732/1‑100 (MQ=255)
                                                                                           ccGGAGCCTGCTCCGGTGAGATCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGCGACt           >  2:1680935/1‑99 (MQ=255)
                                                                                            cGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTgg         <  1:1804801/100‑1 (MQ=255)
                                                                                             ggAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTggg        <  2:1651019/100‑1 (MQ=255)
                                                                                              gAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc       >  2:1046209/1‑100 (MQ=255)
                                                                                              gAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc       <  1:3757957/100‑1 (MQ=255)
                                                                                                  cTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCt   >  2:702829/1‑100 (MQ=255)
                                                                                                   tGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTg  >  2:27251/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTG  >  NC_000913/3156492‑3156690

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: