Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,156,590 C→A 15.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,5900CA15.1% 372.7 / 52.6 133intergenic (+72/‑33)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base C (53/60);  new base A (9/11);  total (62/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 137 total aligned reads displayed.

CCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTC  >  NC_000913/3156491‑3156688
                                                                                                   |                                                                                                  
  gccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTg                                                                                                  <  2:2831332/100‑1 (MQ=255)
      tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTc                                                                                               <  1:1243378/99‑1 (MQ=255)
        tataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCgg                                                                                            >  1:177643/1‑100 (MQ=255)
         ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                           >  1:3216878/1‑100 (MQ=255)
         ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                           <  2:2370989/100‑1 (MQ=255)
         ataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGt                                                                                           <  2:1400627/100‑1 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGt                                                                                           >  1:85314/1‑99 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                          <  1:343380/100‑1 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                          >  2:343382/1‑100 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                          >  1:85252/1‑100 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                          >  2:3308861/1‑100 (MQ=255)
          taCGAAGCCGCCCGCTAAGCTTTTTACCCCTCAAACTTTCGTTTTCTGGCATTTCGTCCAGACTTAAGTTCACAACAACTCACCGGAGCAGGCTCCGGt                                                                                           <  1:3304558/99‑1 (MQ=255)
           aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGa                                                                                         <  1:601737/100‑1 (MQ=255)
           aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGa                                                                                         <  2:1631085/100‑1 (MQ=255)
            cGAAGCCGCCCGCTCACCTTTTTACGCCTCAAACTTTCGTTTCCGGGCATTTCGTCCAGACCTCAGTTCACAACACCCCACCGGCGCAGGCTCCGGTGAg                                                                                        <  1:937592/100‑1 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                        >  2:2189811/1‑100 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                        <  2:1080960/100‑1 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                        <  2:1047643/100‑1 (MQ=255)
            cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                        <  1:2594816/100‑1 (MQ=255)
               aGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTc                                                                                     >  2:2354249/1‑100 (MQ=255)
                gccgccCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCa                                                                                    <  2:2954255/100‑1 (MQ=255)
                gccgccCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCa                                                                                    >  1:1349073/1‑100 (MQ=255)
                  cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCata                                                                                  <  2:3502152/100‑1 (MQ=255)
                   gcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAggtgttgt                                                                                 <  2:1980934/100‑8 (MQ=38)
                    cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                <  2:462288/100‑1 (MQ=255)
                      cgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaaa                                                                              >  2:1299719/1‑100 (MQ=255)
                      cgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaaa                                                                              >  2:2795319/1‑100 (MQ=255)
                       gcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAag                                                                             <  2:2866623/100‑1 (MQ=255)
                        cTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAagg                                                                            <  2:1774352/100‑1 (MQ=255)
                          aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                          >  1:50155/1‑100 (MQ=255)
                           aGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggagg                                                                         >  1:2625270/1‑100 (MQ=255)
                           aGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggagg                                                                         >  2:3574778/1‑100 (MQ=255)
                            gCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAggagga                                                                        >  2:1879797/1‑100 (MQ=255)
                            gCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAggagga                                                                        >  2:598921/1‑100 (MQ=255)
                              tttttACGCCTCAACCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAAc                                                                      <  2:1121383/100‑1 (MQ=255)
                               ttttACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACg                                                                     <  2:552573/100‑1 (MQ=255)
                                  tACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTAt                                                                  <  1:1723301/100‑1 (MQ=255)
                                  tACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTAt                                                                  <  1:1989108/100‑1 (MQ=255)
                                     gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGc                                                               <  1:3361565/100‑1 (MQ=255)
                                       cgcAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGGTCATATAAAGGAGGAACGTATGGc                                                               <  1:1607867/96‑1 (MQ=255)
                                          aaaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACTTATGGCTAATc                                                          <  1:2977271/100‑1 (MQ=255)
                                          aaaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATc                                                          <  1:2097295/100‑1 (MQ=255)
                                           aaCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATcc                                                         <  2:1434725/100‑1 (MQ=255)
                                             cTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCaa                                                       <  1:1456241/100‑1 (MQ=255)
                                             cTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCaa                                                       >  2:1429503/1‑100 (MQ=255)
                                                 cGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACATCACCGGAGCCTGCTTCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCAAACCGt                                                   >  2:500791/1‑100 (MQ=255)
                                                  gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                  <  2:2852937/100‑1 (MQ=255)
                                                    tttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttat                                                <  1:123879/100‑1 (MQ=255)
                                                    tttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttat                                                <  1:1429495/100‑1 (MQ=255)
                                                     ttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGGGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatt                                               >  2:570934/1‑100 (MQ=255)
                                                      tCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatta                                              <  2:1870873/100‑1 (MQ=255)
                                                       cGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTaa                                             >  1:1667475/1‑100 (MQ=255)
                                                       cGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTaa                                             >  2:1237510/1‑100 (MQ=255)
                                                        gggCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGCGCCTGCTCCCGTGAGTTCATATAAACGAGATACGTATGGCTAATCCAACCGttatta                                              <  1:1185443/98‑1 (MQ=255)
                                                         ggCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGc                                           <  2:2435898/100‑1 (MQ=255)
                                                          gCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCt                                          <  2:1097065/100‑1 (MQ=255)
                                                           cATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTa                                         <  1:530469/100‑1 (MQ=255)
                                                            aTTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTAc                                        <  1:2839912/100‑1 (MQ=255)
                                                            aTTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTAc                                        >  2:1884708/1‑100 (MQ=255)
                                                             tttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACa                                       <  2:3023690/100‑1 (MQ=255)
                                                                 gTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                   >  1:2765323/1‑100 (MQ=255)
                                                                    cAGACTTAGGGTCACAACACCACGCCGGAGACTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCACCCGTTATTAAGCTACAGGATGGc                                <  1:2360630/100‑1 (MQ=255)
                                                                     aGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCa                               >  2:1642950/1‑100 (MQ=255)
                                                                        cTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATg                            <  1:1130365/100‑1 (MQ=255)
                                                                         ttAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATgn                           <  2:2401577/100‑2 (MQ=255)
                                                                          tAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTc                          <  1:3028855/100‑1 (MQ=255)
                                                                          tAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTc                          >  2:164377/1‑100 (MQ=255)
                                                                          tAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTc                          >  1:2009215/1‑100 (MQ=255)
                                                                           aaGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCa                         >  2:1486137/1‑100 (MQ=255)
                                                                            ngTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATGCAACCGTTATTAAGCTACAGGATGGCAATGTCAt                        >  1:2041571/2‑100 (MQ=255)
                                                                            aGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCAt                        >  2:1987309/1‑100 (MQ=255)
                                                                                cacaACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGccc                    >  1:1652528/1‑100 (MQ=255)
                                                                                 acaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcccc                   >  2:2022184/1‑100 (MQ=255)
                                                                                 acaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcccc                   >  1:716805/1‑100 (MQ=255)
                                                                                  caacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCa                  >  2:1459527/1‑100 (MQ=255)
                                                                                  caacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCa                  >  2:316727/1‑100 (MQ=255)
                                                                                    tgaacTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGc                >  2:281785/5‑100 (MQ=39)
                                                                                     cacCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAGTGTCATGCCCCAGCt               >  2:3817053/1‑100 (MQ=255)
                                                                                      aacTCACCGGAGCAGGCTCCAGTGAGTTCATATAAAGGAGGAACCTATGGCTAATCAAACCGGTATTAAGCTACAGGATGGCAAAGTCATGCCCCAGCTg              >  1:2365278/3‑100 (MQ=255)
                                                                                       ccTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg             >  1:2998343/1‑100 (MQ=255)
                                                                                       acTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg             >  2:1209799/2‑100 (MQ=255)
                                                                                        cacACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg            >  2:2626222/3‑100 (MQ=255)
                                                                                        cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg            <  2:485299/100‑1 (MQ=255)
                                                                                        cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg            <  1:1612340/100‑1 (MQ=255)
                                                                                        cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg            <  2:3070986/100‑1 (MQ=255)
                                                                                        cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg            <  2:1447996/100‑1 (MQ=255)
                                                                                        cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg            <  2:2765345/100‑1 (MQ=255)
                                                                                        cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg            <  2:2465377/100‑1 (MQ=255)
                                                                                         tCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa           <  1:817581/100‑1 (MQ=255)
                                                                                         tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGcccc                   >  1:304948/1‑92 (MQ=255)
                                                                                         tCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCaccggg             >  2:1637462/1‑93 (MQ=255)
                                                                                          cACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc          >  2:645766/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc          >  2:3288721/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc          >  2:2168226/1‑100 (MQ=255)
                                                                                          cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc          >  2:576064/1‑100 (MQ=255)
                                                                                              ggAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTggg      <  2:3072805/100‑1 (MQ=255)
                                                                                               gAGCCTGCTCCGGTGAGTTCATATAAAGTAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTggg      <  1:2338641/99‑1 (MQ=255)
                                                                                                  ccTGCTCCGGTGAGTTCATATAACGGAGGAACGTATTGCTAATCCCACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCg    <  1:1035596/98‑1 (MQ=255)
                                                                                                  ccTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTc  <  2:3739551/100‑1 (MQ=255)
                                                                                                  ccTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTc  <  1:1299711/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTC  >  NC_000913/3156491‑3156688

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