Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,590 C→A 16.9% 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,5900CA16.9% 471.0 / 71.4 193intergenic (+72/‑33)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base C (74/85);  new base A (17/16);  total (91/102)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

Only 100 of 196 total aligned reads displayed.

CCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCG  >  NC_000913/3156491‑3156686
                                                                                                   |                                                                                                
ccgccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGcc                                                                                                  >  1:1770993/1‑100 (MQ=255)
     gTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTACAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTc                                                                                             >  2:2634532/1‑100 (MQ=255)
     gTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGAGCATTTCATCAAAAATAAAGTTCACAAAACCAAAACGGAGCCTGCTc                                                                                             >  1:1096166/1‑100 (MQ=255)
      tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTcc                                                                                            <  1:714509/100‑1 (MQ=255)
      tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTcc                                                                                            <  2:1510376/100‑1 (MQ=255)
      tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACAGGAGCCTGATcc                                                                                            >  1:850253/1‑100 (MQ=255)
        tataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCgg                                                                                          <  1:2125377/100‑1 (MQ=255)
          taCGAAGCCGCTTGCTAAGCTTTTTACTCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                        <  1:1942216/100‑1 (MQ=255)
          taCGAAGCCGCGCGTTAAGCTTTTTACGCCTCAAACTTTAGTTTTCGGGCATTTTGGCGAGACTTAAGGTGAGAAGACCTTACTGGAGGCTGCTCCGGt                                                                                         >  2:777600/1‑99 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                        >  1:3440726/1‑100 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                        >  1:2991107/1‑100 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                        <  1:2787474/100‑1 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTg                                                                                        <  1:3737350/100‑1 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAAACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                        >  1:3330031/1‑100 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACACCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTg                                                                                        <  2:3707256/100‑1 (MQ=255)
           aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGa                                                                                       <  2:1721868/100‑1 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                      <  1:2442531/100‑1 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                      >  2:598736/1‑100 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                      <  1:3178561/100‑1 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                      <  1:3744307/100‑1 (MQ=255)
             gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                      >  2:430818/1‑99 (MQ=255)
             gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGt                                                                                     <  2:1617940/100‑1 (MQ=255)
             gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAgn                                                                                     <  2:1651076/100‑2 (MQ=255)
              aaGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCTTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGtt                                                                                    <  2:704788/100‑1 (MQ=255)
                 ccgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGgagatcggaa                                                                                 >  1:3117478/1‑94 (MQ=255)
                  cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACAACTCACCGGAGCAGGCTCCGGTGAggtgttg                                                                                <  2:1800499/100‑7 (MQ=37)
                    cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                              <  1:2403586/100‑1 (MQ=255)
                       gcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAag                                                                           <  2:3129168/100‑1 (MQ=255)
                       gcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAag                                                                           <  1:525704/100‑1 (MQ=255)
                       gcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAag                                                                           <  2:3214151/100‑1 (MQ=255)
                          aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                        >  2:502406/1‑100 (MQ=255)
                           aGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggagg                                                                       <  2:880052/100‑1 (MQ=255)
                             cTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGaa                                                                     <  2:52910/100‑1 (MQ=255)
                              tttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                    <  1:1208761/100‑1 (MQ=255)
                              tttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                    <  2:2970624/100‑1 (MQ=255)
                               ttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACg                                                                   <  2:673812/100‑1 (MQ=255)
                               ttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACg                                                                   <  2:1281904/100‑1 (MQ=255)
                                 ttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGTCTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTa                                                                 <  1:1781029/100‑1 (MQ=255)
                                  tACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTAt                                                                <  1:3263780/100‑1 (MQ=255)
                                   aCGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATg                                                               <  2:1032867/100‑1 (MQ=255)
                                      ccTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCt                                                            >  1:1881290/1‑100 (MQ=255)
                                      ccTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCt                                                            <  2:2619705/100‑1 (MQ=255)
                                       cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTa                                                           >  1:2080596/1‑100 (MQ=255)
                                       cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAACCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTa                                                           >  1:333098/1‑100 (MQ=255)
                                         cAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAAt                                                         <  2:850243/100‑1 (MQ=255)
                                         cAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAAt                                                         >  2:1529012/1‑100 (MQ=255)
                                             cTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAACCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCaa                                                     >  1:1604936/1‑100 (MQ=255)
                                             cTTTCGATTTCGAGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACATATGGCTAATCCaa                                                     >  2:1203202/1‑100 (MQ=255)
                                                tCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCg                                                  <  2:280950/100‑1 (MQ=255)
                                                 cGTTTTCGGGCATTTCGTCCAGCCTTAAGTTCACAACCCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGt                                                 <  2:2641264/100‑1 (MQ=255)
                                                  gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                >  2:3884436/1‑100 (MQ=255)
                                                  gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                >  1:1712651/1‑100 (MQ=255)
                                                    tttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGATGAACGTATGGCTAATCCAACCGttat                                              <  2:2153089/100‑1 (MQ=255)
                                                      tCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatta                                            >  2:402972/1‑100 (MQ=255)
                                                       cGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTaa                                           <  2:110932/100‑1 (MQ=255)
                                                         ggCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGACCCTGCTCCGGTGAGTTCATATAAAGGAGGCACGTATGGCTAATCCAACAGTTATTAAGc                                         >  2:2109779/1‑100 (MQ=255)
                                                              ttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAg                                    <  1:1851627/100‑1 (MQ=255)
                                                              ttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAg                                    >  2:1010118/1‑100 (MQ=255)
                                                                 gTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                 >  1:948747/1‑100 (MQ=255)
                                                                 gTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                 >  1:2467634/1‑100 (MQ=255)
                                                                 gTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                 >  2:1977398/1‑100 (MQ=255)
                                                                 gTCCAGACTTAAGTGCACAACACCTCACCTGGGCCTGCTCCGGTGAGTTCATATAGAGGGGGTACGTATTGCTAATCTGACCGTTAGTAAGCTACAGGAt                                 >  2:172621/1‑100 (MQ=255)
                                                                  tCCAGACTTAAGTTCACAACACCTCTCTGGTGCCTGCTCTGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATg                                <  1:1579013/100‑1 (MQ=255)
                                                                   ccAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATgg                               <  1:1448314/100‑1 (MQ=255)
                                                                   ccAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATgg                               <  2:3784186/100‑1 (MQ=255)
                                                                     aGACTTGAGTTCACAACACCTCGCCGGTGTCTGCTGCGGTGAGTTCAGATAAATGCGTAACGTATGGCTAAGCCAACCGTTATTAAGCTACAGGATGGc                              <  2:1922723/99‑1 (MQ=255)
                                                                        cTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATg                          >  1:2921054/1‑100 (MQ=255)
                                                                            aGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCAt                      <  1:498264/100‑1 (MQ=255)
                                                                             gTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATg                     <  1:2837616/100‑1 (MQ=255)
                                                                              ttCACAACACCTCACCGGAGCCTGCTCCGGTGTGTTCGTATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGc                    >  2:2084355/1‑100 (MQ=255)
                                                                              ttCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATgt                    <  2:1449679/100‑2 (MQ=255)
                                                                               ctacaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcc                   <  2:1319245/98‑1 (MQ=255)
                                                                                 acaacaCCTCACCGGAGCCTGCTCCTGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcccc                 <  1:624664/100‑1 (MQ=255)
                                                                                  caacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCa                <  1:2517345/100‑1 (MQ=255)
                                                                                   aacaCCTCACCGGAGCCTGCTCCTGTTAGGTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCa                <  2:3609517/99‑1 (MQ=255)
                                                                                    acacCTCACCGGATCCTGCTCCTGTGAGTTCATATAAAGGATGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAgt              <  2:1047871/100‑2 (MQ=255)
                                                                                    acacCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGTAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACATGATGGCAATGTCATGCCCCAGc              <  2:2421949/100‑1 (MQ=255)
                                                                                      aacTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTg            >  1:1794733/3‑100 (MQ=255)
                                                                                       ccTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAACTACAAGATGGCAATGTCATGCCCCAGCTgg           >  1:2139819/1‑100 (MQ=255)
                                                                                       acTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg           >  2:3011785/2‑100 (MQ=255)
                                                                                        cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg          <  2:2467611/100‑1 (MQ=255)
                                                                                        cTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg          >  1:3911842/1‑100 (MQ=255)
                                                                                        cTCACCGGAGCAGGCACCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTCCAGGATGGCAATGTCATGCCCCAGCTggg          >  2:2526575/1‑100 (MQ=255)
                                                                                         tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa         >  1:1934044/1‑100 (MQ=255)
                                                                                         tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa         >  1:261700/1‑100 (MQ=255)
                                                                                         tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGAATGGCAATGTCATACCCCAGCTGGGa         >  1:2389716/1‑100 (MQ=255)
                                                                                         tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgaga         >  2:3741727/1‑100 (MQ=255)
                                                                                         tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa         >  1:2595774/1‑100 (MQ=255)
                                                                                         tCACCGGAGCAGGCTACGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCTCTAGCTGGGa         >  2:54075/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGTAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc        >  1:2753204/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc        >  2:309297/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATACAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg           >  2:1137862/1‑97 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGATCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc        >  2:2383041/1‑100 (MQ=255)
                                                                                           aCCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACt       >  1:1506672/1‑100 (MQ=255)
                                                                                            ccGGCGCCTGCTCCGGCGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa         >  2:3383019/1‑97 (MQ=255)
                                                                                            ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACt       >  2:1352555/1‑99 (MQ=255)
                                                                                              ggAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTggg    >  2:2553764/1‑100 (MQ=255)
                                                                                               gAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc   <  1:3520877/100‑1 (MQ=255)
                                                                                               gAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc   >  1:823537/1‑100 (MQ=255)
                                                                                                aGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCg  <  1:1305800/100‑1 (MQ=255)
                                                                                                   |                                                                                                
CCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCG  >  NC_000913/3156491‑3156686

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