Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,590 C→A 20.6% intergenic (+72/‑33) 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,5900CA20.6% 262.0 / 57.5 116intergenic (+72/‑33)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base C (43/49);  new base A (10/14);  total (53/63)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 118 total aligned reads displayed.

CCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCT  >  NC_000913/3156491‑3156689
                                                                                                   |                                                                                                   
ccgccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCACACTTGCGGTGTCGGGCAGTTCGTCCAGACGGAAGTTCACAACACCTCACCGGAGcc                                                                                                     <  1:1181044/100‑1 (MQ=255)
  gccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTg                                                                                                   >  1:1092470/1‑100 (MQ=255)
  gccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTg                                                                                                   >  2:2090298/1‑100 (MQ=255)
    cgTATATACGAAGCCGCCCGCTAAGCTTGTTACGCATCAAACTTTCGTTATCGGGCATTTCGTCCAAACTTAAGTGCACAACACCTCACCGGAGCCTGCt                                                                                                 <  1:2338672/100‑1 (MQ=255)
     gTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTc                                                                                                >  2:1596076/1‑100 (MQ=255)
         ataCGAAGCCGCCCGCTAAGCTTTTTACGTCTCAAACTTTCGTTTTCGGACCTTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                            <  1:1605794/100‑1 (MQ=255)
         ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGt                                                                                            <  2:1386444/100‑1 (MQ=255)
          taCGAGGCCGCCCGCTAAGCTTTTTACGCCTCAACCCTTCGTCTTCGGCCATTTCGACCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                           <  1:2550296/100‑1 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTg                                                                                           <  2:1730838/100‑1 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTg                                                                                           <  2:2493542/100‑1 (MQ=255)
           aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGa                                                                                          <  1:42939/100‑1 (MQ=255)
           aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGa                                                                                          <  1:1614841/100‑1 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                         >  1:1767003/1‑100 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                         <  1:1517238/100‑1 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                         <  2:2143623/100‑1 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                         <  2:413777/100‑1 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                         <  1:2015943/100‑1 (MQ=255)
             gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAgg                                                                                        <  2:1899360/100‑2 (MQ=255)
             gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAgg                                                                                        <  2:23307/100‑2 (MQ=255)
             gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAgg                                                                                        <  1:336711/100‑2 (MQ=255)
             gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTGCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAgg                                                                                        <  1:2169500/100‑2 (MQ=255)
                gccgccTGCTAAGCTTTTTACGCCTTAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                            <  1:2452296/93‑1 (MQ=255)
                gccgccCGCTAAGCTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCat                                                                                    >  2:1001749/1‑100 (MQ=255)
                 ccgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCat                                                                                    >  1:2128227/1‑100 (MQ=255)
                 ccgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCat                                                                                    >  2:1731473/1‑100 (MQ=255)
                 ccgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCat                                                                                    <  2:48093/100‑1 (MQ=255)
                    cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                 <  1:2090284/100‑1 (MQ=255)
                    cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                 <  1:548078/100‑1 (MQ=255)
                    cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACAACTCACCGGAGCAGGCTCCGGTGAggtgttgtg                                                                                 <  1:1926478/100‑9 (MQ=37)
                     ccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACCTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                 >  1:1884470/1‑99 (MQ=255)
                       gcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAag                                                                              >  2:900083/1‑100 (MQ=255)
                        cTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagg                                                                             <  2:2091904/100‑1 (MQ=255)
                          aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                           <  2:215076/100‑1 (MQ=255)
                           aGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggagg                                                                          >  1:935727/1‑100 (MQ=255)
                           aGCTTTTCACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggagg                                                                          >  1:1732568/1‑100 (MQ=255)
                            gCTTTTTACGCCTGAAACTTGCGATTTCGGGCATTTCGTCCAGACTTAAGTTCACAGCACCTCACCGGAGCCTTCGCCGGTGAGTTCATATAAagga                                                                            >  1:1576347/1‑97 (MQ=255)
                              ttgttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                       <  1:2204071/97‑1 (MQ=255)
                               ttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACg                                                                      <  1:869204/100‑1 (MQ=255)
                               ttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACg                                                                      <  1:1731464/100‑1 (MQ=255)
                                tttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGt                                                                     >  1:1313059/1‑100 (MQ=255)
                                   aCGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATg                                                                  <  2:309165/100‑1 (MQ=255)
                                   aCGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATg                                                                  <  1:785977/100‑1 (MQ=255)
                                    cGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATgg                                                                 >  1:1350666/1‑100 (MQ=255)
                                     gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                                >  2:1767160/1‑100 (MQ=255)
                                      ccTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCt                                                               >  1:2019960/1‑100 (MQ=255)
                                      ccTCAAACTTTCGTTTTCGGGCATTTCGCCCAGACTTAAGTTCACAACACCTCACCGGAGGCTGCTCCGGTGAGTTCATATACAGTAGGAACGTATGGc                                                                <  1:1598902/99‑1 (MQ=255)
                                       cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTa                                                              >  2:1512227/1‑100 (MQ=255)
                                          aaaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATc                                                           >  1:510087/1‑100 (MQ=255)
                                          aaaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATc                                                           >  2:1251047/1‑100 (MQ=255)
                                           aaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATcc                                                          <  2:432395/100‑1 (MQ=255)
                                              tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACATCACCTCACCGGAGCCTGCTTCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAc                                                       <  1:1413381/100‑1 (MQ=255)
                                              tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAc                                                       <  2:2502563/100‑1 (MQ=255)
                                               ttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAcc                                                      <  2:541513/100‑1 (MQ=255)
                                               ttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAcc                                                      <  2:1191610/100‑1 (MQ=255)
                                                  gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                   >  1:2436851/1‑100 (MQ=255)
                                                   ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtta                                                  >  2:267284/1‑100 (MQ=255)
                                                    tttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttat                                                 <  2:812731/100‑1 (MQ=255)
                                                    tttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttat                                                 <  2:2101977/100‑1 (MQ=255)
                                                        gggCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAg                                             <  2:935730/100‑1 (MQ=255)
                                                        gagCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGACTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAg                                             >  2:1747977/3‑100 (MQ=255)
                                                         ggCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGc                                            <  2:2121120/100‑1 (MQ=255)
                                                         ggCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATAGAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAg                                             >  2:414645/1‑99 (MQ=255)
                                                             tttCGTCCAGACTTAAGTTCACAACCCCTCACCGAAGCCTGCTCCGGTGAGTCCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTAc                                         <  1:1055200/99‑1 (MQ=255)
                                                             tttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACa                                        <  2:1092474/100‑1 (MQ=255)
                                                                 gTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                    >  2:545998/1‑100 (MQ=255)
                                                                 gTCCAGACATAAGAACACAACACCTCACAGGAGCCTGCTCCGGTGAGTCCACAAAAAGGAGGAACGTATGGCCAATCCAACAGTTATTAAGCTACAGGAt                                    >  1:1398707/1‑100 (MQ=255)
                                                                      gACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCaa                               <  2:22941/100‑1 (MQ=255)
                                                                        cTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGCACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATg                             >  1:953352/1‑100 (MQ=255)
                                                                            aGTTGACAGCACCTCGCCGGAGCCTGCTCCGGGGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTGATTAAGCTACAGGATGGCAATGTc                           <  1:2413625/98‑1 (MQ=255)
                                                                            aGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCAt                         >  2:1045931/1‑100 (MQ=255)
                                                                             gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                        >  2:468412/1‑100 (MQ=255)
                                                                             gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                        <  2:1217605/100‑1 (MQ=255)
                                                                             gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                        <  1:662910/100‑1 (MQ=255)
                                                                             gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGACAATGTCATg                        >  2:1948739/1‑100 (MQ=255)
                                                                             gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAAGCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                        <  2:2371771/100‑1 (MQ=255)
                                                                              ttCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGc                       >  2:1524079/1‑100 (MQ=255)
                                                                               tCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcc                      <  2:70579/100‑1 (MQ=255)
                                                                               tCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcc                      <  2:2054391/100‑1 (MQ=255)
                                                                                cacaACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCCAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGccc                     >  1:1588281/1‑100 (MQ=255)
                                                                                 acaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGCGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcccc                    <  1:1657727/100‑1 (MQ=255)
                                                                                      aacTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTg               >  2:180102/3‑100 (MQ=255)
                                                                                       ccTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTg               >  1:209978/1‑99 (MQ=255)
                                                                                       ccTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTg               >  2:342923/1‑99 (MQ=255)
                                                                                       acTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg              >  1:558026/2‑100 (MQ=255)
                                                                                        cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTg               <  1:1051709/98‑1 (MQ=255)
                                                                                        cTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg             >  1:2374559/1‑100 (MQ=255)
                                                                                        cTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg             >  2:2477330/1‑100 (MQ=255)
                                                                                        cTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg             >  1:1110112/1‑100 (MQ=255)
                                                                                        atctccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg             <  2:865972/96‑1 (MQ=255)
                                                                                         tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGAAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa            >  1:1375642/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAATGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  2:2458587/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  1:2290950/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  1:2371281/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  2:474117/1‑100 (MQ=255)
                                                                                            ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGagat         >  1:865969/1‑97 (MQ=255)
                                                                                            ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTg         >  2:1145621/1‑100 (MQ=255)
                                                                                            ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTg         <  2:1597778/100‑1 (MQ=255)
                                                                                             cGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTgg        <  2:1554263/100‑1 (MQ=255)
                                                                                                aGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCg     <  1:1145617/100‑1 (MQ=255)
                                                                                                 gCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGt    >  2:289123/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCT  >  NC_000913/3156491‑3156689

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: