Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,590 C→A 20.1% 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.1% 341.0 / 74.5 138intergenic (+72/‑33)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base C (53/57);  new base A (18/10);  total (71/67)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 148 total aligned reads displayed.

CCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCT  >  NC_000913/3156491‑3156689
                                                                                                   |                                                                                                   
ccgccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGcc                                                                                                     <  1:2371701/100‑1 (MQ=255)
ccgccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGcc                                                                                                     <  1:549879/100‑1 (MQ=255)
ccgccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGcc                                                                                                     <  1:1867713/100‑1 (MQ=255)
  gccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTg                                                                                                   <  1:461493/100‑1 (MQ=255)
      tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTcc                                                                                               <  1:2412721/100‑1 (MQ=255)
      tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTcc                                                                                               <  1:2845372/100‑1 (MQ=255)
         ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACTCCTCACCGGAGCAGGCTCCGGt                                                                                            <  2:2114123/100‑1 (MQ=255)
         ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                            <  2:1023097/100‑1 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                           <  1:1492470/100‑1 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                           >  2:2848310/1‑100 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                           >  2:1294768/1‑100 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTg                                                                                           <  1:485979/100‑1 (MQ=255)
           aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGa                                                                                          <  1:1085450/100‑1 (MQ=255)
            cGAAGCGGGCCGCTAAGCTTTTTCCGCCTCAAACTTCCGTTTTCGGGCATGTCGTCCAGACTTAAGTTCACAACACCTCGCCGGAGCCTGCTcc                                                                                               <  2:3395389/94‑1 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                         >  1:1679660/1‑100 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                         >  1:2748352/1‑100 (MQ=255)
             gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGt                                                                                        <  2:2124476/100‑1 (MQ=255)
             gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGt                                                                                        <  1:690135/100‑1 (MQ=255)
             gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAgg                                                                                        <  1:1766964/100‑2 (MQ=255)
             gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAgg                                                                                        <  2:1520523/100‑2 (MQ=255)
             gAAGCCGCCCGCTAAGCTCTTTACTCATCAAACTTTCTCTCTCGGGCCTTTCGTCCAGACTTAACTTCCCCACACCTCACCGGCTCCTGCTCCGGTg                                                                                           <  2:2245522/97‑1 (MQ=255)
              aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGTATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAggt                                                                                       <  2:485054/100‑3 (MQ=39)
              aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCTCCGGAGCCTGCTCCGGTGAGtt                                                                                       <  1:658088/100‑1 (MQ=255)
              aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGtt                                                                                       <  1:1615813/100‑1 (MQ=255)
              aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGtt                                                                                       <  1:2640310/100‑1 (MQ=255)
               agcagcCCGCTAAGCTTTTGACGCCTTAAACTTTCGTTTTCGGGCATTTCGTCGAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAggtg                                                                                      <  2:2020799/100‑4 (MQ=37)
                 ccgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCa                                                                                     <  2:3091821/99‑1 (MQ=255)
                        cTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagg                                                                             <  1:28852/100‑1 (MQ=255)
                        cTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagg                                                                             <  2:852803/100‑1 (MQ=255)
                        cTAAGCGTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagg                                                                             <  1:2702045/100‑1 (MQ=255)
                         tAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagga                                                                            >  1:198707/1‑100 (MQ=255)
                          aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                           <  2:1767009/100‑1 (MQ=255)
                           aGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGGCTTAAGTTCACAACACCTCACCGGCGCCTGCTCCGGTGAGTTCATATAAaggagg                                                                          >  2:1624239/1‑100 (MQ=255)
                            gCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAggagga                                                                         >  1:3430835/1‑100 (MQ=255)
                              tttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                       >  1:2022933/1‑100 (MQ=255)
                              tttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                       <  1:3349835/100‑1 (MQ=255)
                              tttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                       <  1:3236542/100‑1 (MQ=255)
                                 ttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTa                                                                    >  1:1365432/1‑100 (MQ=255)
                                     gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                                >  2:2283641/1‑100 (MQ=255)
                                     gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                                <  1:2589836/100‑1 (MQ=255)
                                     gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                                <  1:2421166/100‑1 (MQ=255)
                                           aaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATcc                                                          <  1:1586978/100‑1 (MQ=255)
                                            aCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATcc                                                          <  2:1385972/99‑1 (MQ=255)
                                             cgttCGTCGTCGGTCATTTCGTCCAGACTTAACTTCACCACACCCCCCCGGCACCTACGGCGGTGAGCTCATAAAAACCCGGAACGTATGCCTAATCCaa                                                        >  2:595975/3‑100 (MQ=255)
                                             cTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATcc                                                          >  2:2214215/1‑98 (MQ=255)
                                             cTTTCGTTTTAGGGAATTTCGTATAGACTTAAGTTCACAATACCTTATCTGAGCCTTCTACGGTGAGTTCATATAAAGGATGAATTTATGGTTAAt                                                            >  2:242113/1‑96 (MQ=255)
                                                tCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCg                                                     <  2:1740241/100‑1 (MQ=255)
                                                  gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                   >  1:384289/1‑100 (MQ=255)
                                                  gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                   <  2:2648336/100‑1 (MQ=255)
                                                    tttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtta                                                  <  2:1092000/99‑1 (MQ=255)
                                                      tCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCGCCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatta                                               >  2:1955231/1‑100 (MQ=255)
                                                       cGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTaa                                              >  2:3541610/1‑100 (MQ=255)
                                                         ggCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGc                                            <  1:751854/100‑1 (MQ=255)
                                                         ggCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGc                                            <  1:1134979/100‑1 (MQ=255)
                                                              ttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAg                                       >  1:3076009/1‑100 (MQ=255)
                                                              ttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAg                                       <  2:3403403/100‑1 (MQ=255)
                                                                  tCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATg                                   >  2:279778/1‑100 (MQ=255)
                                                                   ccAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATgg                                  >  2:2955748/1‑100 (MQ=255)
                                                                   ccAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATgg                                  >  1:1632455/1‑100 (MQ=255)
                                                                   ccAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATgg                                  >  1:1884219/1‑100 (MQ=255)
                                                                    cAGACTTAAGTCCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGc                                 >  2:1612238/1‑100 (MQ=255)
                                                                     aGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCa                                <  1:3606046/100‑1 (MQ=255)
                                                                        cTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATg                             >  1:539795/1‑100 (MQ=255)
                                                                        cTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATg                             <  1:2283657/100‑1 (MQ=255)
                                                                          tAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGATAATCCAACCGTTATTAAGCTACAGGATGGCAATGt                            <  2:1473503/99‑1 (MQ=255)
                                                                            aGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATACAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCa                          <  2:2331689/99‑1 (MQ=255)
                                                                             gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                        >  2:2194528/1‑100 (MQ=255)
                                                                             gTTCACAACACCTCACCGGAGCCTGCACCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                        >  1:2222811/1‑100 (MQ=255)
                                                                               tCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATGCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcc                      <  2:1949541/100‑1 (MQ=255)
                                                                               tCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcc                      >  1:540054/1‑100 (MQ=255)
                                                                               tCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcc                      <  1:2657273/100‑1 (MQ=255)
                                                                               tCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcc                      >  2:1300628/1‑100 (MQ=255)
                                                                                 acaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcccc                    <  1:1304147/100‑1 (MQ=255)
                                                                                  tatgacCTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGCACGTATGGCTCGTCCATCCGTTCCTAACATGCCCGATGGCCCTGTCCTGccc                     >  2:2564419/5‑98 (MQ=12)
                                                                                   aacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAg                  <  2:388522/100‑1 (MQ=255)
                                                                                       ccTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg              >  1:1919559/1‑100 (MQ=255)
                                                                                       ccTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg              >  2:1157264/1‑100 (MQ=255)
                                                                                       ccTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg              >  1:3126796/1‑100 (MQ=255)
                                                                                        cTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg             >  1:2864815/1‑100 (MQ=255)
                                                                                        cTCACCGGAGCAGGCTCCGGTGAGATCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg             >  1:1501111/1‑100 (MQ=255)
                                                                                         tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGCGGAACGTCTGGCTAATACAACCGTTATTAAGCTACAGGTTGGCAATGTCATAACCCAGCTggg             >  2:1194186/1‑99 (MQ=255)
                                                                                         tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa            >  1:2586768/1‑100 (MQ=255)
                                                                                         tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa            >  1:1998956/1‑100 (MQ=255)
                                                                                         tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa            >  1:2842471/1‑100 (MQ=255)
                                                                                         tCACCGGAGCAGGCTCCGGTGAAGTCGTAAAACGCAGCCACGTAGGCCTAATGCAATCGTCATTAAGCTACAGGATGGCAAGGACATTCTCCAACTGGGa            >  2:2956306/1‑100 (MQ=255)
                                                                                          cACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  2:1412550/1‑100 (MQ=255)
                                                                                          cACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           <  1:2720728/100‑1 (MQ=255)
                                                                                          cACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  1:2784475/1‑100 (MQ=255)
                                                                                          cACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  1:3520587/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  1:1818078/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  2:1124556/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc           >  1:915827/1‑100 (MQ=255)
                                                                                           aCCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACt          <  1:1717064/100‑1 (MQ=255)
                                                                                             cGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTgg        <  2:3098608/100‑1 (MQ=255)
                                                                                              ggAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTggg       >  1:1024492/1‑100 (MQ=255)
                                                                                              ggAGCCTGCGCCGGTGAGTTCATATAAAGGAGTAACGTATGGCTAATCCAACCGTTATTAAGCGACAGGATGGCAATGTCATGCCCCAGCTGGGACTggg       <  1:1971937/100‑1 (MQ=255)
                                                                                                 gCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGt    >  1:2283506/1‑100 (MQ=255)
                                                                                                  ccTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTc   >  2:808609/1‑100 (MQ=255)
                                                                                                  ccTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATCTCATGCCCCAGCTGGGACTGGGCGTc   <  1:2214229/100‑1 (MQ=255)
                                                                                                   cTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCt  >  1:558116/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: