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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,156,5910TG16.5% 447.3 / 87.2 194intergenic (+73/‑32)yqhD/dkgAaldehyde reductase, NADPH‑dependent/2,5‑diketo‑D‑gluconate reductase A
Reads supporting (aligned to +/- strand):  ref base T (75/87);  new base G (16/16);  total (91/103)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-01

Only 100 of 200 total aligned reads displayed.

CCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTG  >  NC_000913/3156494‑3156690
                                                                                                 |                                                                                                   
 cgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCt                                                                                                  <  1:983108/100‑1 (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)
     ttccgatctccgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCgg                                                                                              <  2:3117448/91‑1 (MQ=255)
     tataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCgg                                                                                              <  1:2125377/100‑1 (MQ=255)
     tataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTAGAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCgg                                                                                              >  1:1746102/1‑100 (MQ=255)
       taCGAAGCCGCTTGCTAAGCTTTTTACTCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTg                                                                                            <  1:1942216/100‑1 (MQ=255)
       taCGAAGCCGCGCGTTAAGCTTTTTACGCCTCAAACTTTAGTTTTCGGGCATTTTGGCGAGACTTAAGGTGAGAAGACCTTACTGGAGGCTGCTCCGGt                                                                                             >  2:777600/1‑99 (MQ=255)
       taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTg                                                                                            <  1:3737350/100‑1 (MQ=255)
       taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGGCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGtn                                                                                            <  1:2909737/100‑2 (MQ=255)
       taCGAAGCCGCCCGCTAAGCTTTTTACACCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTg                                                                                            <  2:3707256/100‑1 (MQ=255)
       gacgttgccgccCGCTAAGCGTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTg                                                                                            <  2:304757/94‑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:3330935/100‑1 (MQ=255)
         cGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAg                                                                                          <  1:3744307/100‑1 (MQ=255)
          gAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg                                                                                          >  2:430818/1‑99 (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:2390857/100‑7 (MQ=38)
               cgcccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACAACTCACCGGAGCAGGCTCCGGTGAggtgttg                                                                                    <  2:1800499/100‑7 (MQ=37)
                 cctgcTCAGCTTTTTACGCCTCAAACGTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACTCCTCACCGGAGCCTGGTCCGGTGAGTTCAtata                                                                                  <  2:1037520/97‑1 (MQ=255)
                 cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                  >  1:3397305/1‑100 (MQ=255)
                 cccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                  <  1:2403586/100‑1 (MQ=255)
                  ccgcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGTCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCAtata                                                                                  <  1:326089/99‑1 (MQ=255)
                    gcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAag                                                                               >  2:1054518/1‑100 (MQ=255)
                    gcTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAag                                                                               <  1:525704/100‑1 (MQ=255)
                       aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                            >  2:502406/1‑100 (MQ=255)
                       aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                            <  1:3299766/100‑1 (MQ=255)
                       aaGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAaggag                                                                            <  2:934757/100‑1 (MQ=255)
                          cTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGaa                                                                         <  2:52910/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                                                                    <  2:2449852/100‑1 (MQ=255)
                                  gCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGCGGAACGTATGGc                                                                 >  2:2255369/1‑100 (MQ=255)
                                   ccTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCt                                                                >  1:1244590/1‑100 (MQ=255)
                                   ccTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCt                                                                >  1:1881290/1‑100 (MQ=255)
                                    cTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTa                                                               >  1:2080596/1‑100 (MQ=255)
                                      cAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAAt                                                             <  2:850243/100‑1 (MQ=255)
                                         aCTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCa                                                          <  2:3173716/100‑1 (MQ=255)
                                          cTTTCGATTTCGAGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACATATGGCTAATCCaa                                                         >  2:1203202/1‑100 (MQ=255)
                                           tttCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAAc                                                        >  2:2558664/1‑100 (MQ=255)
                                             tCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCg                                                      <  2:280950/100‑1 (MQ=255)
                                              cGTTTTCGGGCATTTCGTCCAGCCTTAAGTTCACAACCCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGt                                                     <  2:2641264/100‑1 (MQ=255)
                                               gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                    >  2:1552474/1‑100 (MQ=255)
                                               gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                    >  1:2794733/1‑100 (MQ=255)
                                               gTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtt                                                    >  2:3884436/1‑100 (MQ=255)
                                                 tttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGATGAACGTATGGCTAATCCAACCGttat                                                  <  2:2153089/100‑1 (MQ=255)
                                                 tttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttat                                                  <  1:1639202/100‑1 (MQ=255)
                                                   tCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttatta                                                >  2:402972/1‑100 (MQ=255)
                                                      ggCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGACCCTGCTCCGGTGAGTTCATATAAAGGAGGCACGTATGGCTAATCCAACAGTTATTAAGc                                             >  2:2109779/1‑100 (MQ=255)
                                                          tttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATGTAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACa                                         <  2:2991080/100‑1 (MQ=255)
                                                           ttCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAg                                        <  1:1851627/100‑1 (MQ=255)
                                                              gTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                     >  1:2467634/1‑100 (MQ=255)
                                                              gTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGAt                                     >  2:1977398/1‑100 (MQ=255)
                                                              gTCCAGACTTAAGTGCACAACACCTCACCTGGGCCTGCTCCGGTGAGTTCATATAGAGGGGGTACGTATTGCTAATCTGACCGTTAGTAAGCTACAGGAt                                     >  2:172621/1‑100 (MQ=255)
                                                                ccAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATgg                                   <  2:3784186/100‑1 (MQ=255)
                                                                 cAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCATCCGTTATTAAGCTACAGGATGGc                                  <  2:254680/100‑1 (MQ=255)
                                                                  aGACTTGAGTTCACAACACCTCGCCGGTGTCTGCTGCGGTGAGTTCAGATAAATGCGTAACGTATGGCTAAGCCAACCGTTATTAAGCTACAGGATGGc                                  <  2:1922723/99‑1 (MQ=255)
                                                                  aGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCa                                 <  1:155577/100‑1 (MQ=255)
                                                                   gACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCaa                                >  2:2517322/1‑100 (MQ=255)
                                                                    aCTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAAt                               <  2:2620798/100‑1 (MQ=255)
                                                                         aGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCAt                          <  1:498264/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)
                                                                              acaacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGccag                     >  1:1319259/1‑98 (MQ=255)
                                                                               caacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCa                    >  2:1020856/1‑100 (MQ=255)
                                                                               caacaCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCa                    <  1:2517345/100‑1 (MQ=255)
                                                                                aacaCCTCACCGGAGCCTGCTCCTGTTAGGTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCa                    <  2:3609517/99‑1 (MQ=255)
                                                                                   aacTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTg                >  1:1794733/3‑100 (MQ=255)
                                                                                    ccTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAACTACAAGATGGCAATGTCATGCCCCAGCTgg               >  1:2139819/1‑100 (MQ=255)
                                                                                     cTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg              >  1:3911842/1‑100 (MQ=255)
                                                                                     cTCACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTggg              >  1:2935799/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)
                                                                                       cACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc            <  1:1579494/100‑1 (MQ=255)
                                                                                       cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGTAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAc            >  1:2753204/1‑100 (MQ=255)
                                                                                       cACCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATACAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTgg               >  2:1137862/1‑97 (MQ=255)
                                                                                         ccGGCGCCTGCTCCGGCGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGa             >  2:3383019/1‑97 (MQ=255)
                                                                                         ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACt           >  2:1352555/1‑99 (MQ=255)
                                                                                         ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTg          >  1:946639/1‑100 (MQ=255)
                                                                                         ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTg          >  1:3104002/1‑100 (MQ=255)
                                                                                         ccGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGAATGACAATGTCATGCCCCAGCTGGGACt           >  1:1497022/1‑99 (MQ=255)
                                                                                            gAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc       <  1:3520877/100‑1 (MQ=255)
                                                                                            gAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGc       >  1:823537/1‑100 (MQ=255)
                                                                                             aGCCTGCTCCGGTGAGTTCATATAAATGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCg      <  2:1518351/100‑1 (MQ=255)
                                                                                                 tGCTCCGGTGAGTTCATATAAATGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTg  <  2:861463/100‑1 (MQ=255)
                                                                                                 tGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTg  <  1:3505790/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
CCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTG  >  NC_000913/3156494‑3156690

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: