Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,891,874 T→C 6.9% Y346H (TAC→CAC)  ttuD → hydroxypyruvate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,891,8740TC 6.9% 252.9 / 6.9 87Y346H (TAC→CAC) ttuDhydroxypyruvate reductase
Reads supporting (aligned to +/- strand):  ref base T (37/44);  new base C (2/4);  total (39/48)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

GTTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGC  >  NC_002947/4891727‑4892022
                                                                                                                                                   |                                                                                                                                                    
gTTCAACACGGCCAGCCGCTCAACGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCCCCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCgtgt                                                                                                                                                      <  4:235500/148‑1 (MQ=255)
gTTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTa                                                                                                                                                     <  1:63009/149‑1 (MQ=255)
    aaCACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCAGGGCGTACACGcc                                                                                                                                                 >  5:391110/1‑149 (MQ=255)
       acGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTg                                                                                                                                              <  6:24621/149‑1 (MQ=255)
       acGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTg                                                                                                                                              <  3:52458/149‑1 (MQ=255)
       acGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTg                                                                                                                                              <  7:329454/149‑1 (MQ=255)
       acGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCAGGGCGTACACGcccgg                                                                                                                                              <  6:391110/149‑3 (MQ=255)
         ggCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCAGGGCGTACACGcccggca                                                                                                                                            <  8:212191/149‑5 (MQ=255)
          gCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCGCCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTACACGcccggcag                                                                                                                                           <  6:335378/149‑6 (MQ=255)
           ccAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGTCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCg                                                                                                                                          <  1:275580/149‑1 (MQ=255)
           ccAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCg                                                                                                                                          <  7:16187/149‑1 (MQ=255)
              gCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTg                                                                                                                                       <  5:14633/149‑1 (MQ=255)
               ccGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTAc                                                                                                                                                    >  8:197296/1‑135 (MQ=255)
               ccGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTAc                                                                                                                                                    <  7:197296/135‑1 (MQ=255)
               ccGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGa                                                                                                                                      >  3:152583/1‑149 (MQ=255)
                 gCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGaca                                                                                                                                    >  3:91511/1‑149 (MQ=255)
                 gCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGaca                                                                                                                                    <  4:315351/149‑1 (MQ=255)
                 gCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGaca                                                                                                                                    <  8:301693/149‑1 (MQ=255)
                   tCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACAcc                                                                                                                                  >  7:302486/1‑149 (MQ=255)
                        gcgcCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACgg                                                                                                                             <  4:91511/149‑1 (MQ=255)
                           cccTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCAt                                                                                                                          >  1:299709/1‑149 (MQ=255)
                            ccTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATc                                                                                                                         >  4:225147/1‑149 (MQ=255)
                             cTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCg                                                                                                                        >  4:252804/1‑149 (MQ=255)
                                 gTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATgg                                                                                                                    <  4:152583/149‑1 (MQ=255)
                                   gATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCt                                                                                                                  >  2:365189/1‑149 (MQ=255)
                                          tCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCgtgt                                                                                                                                                      >  2:353972/1‑106 (MQ=255)
                                          tCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCgtgt                                                                                                                                                      <  1:353972/106‑1 (MQ=255)
                                          tCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGgaag                                                                                                            <  4:35829/148‑1 (MQ=255)
                                              gTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGaaa                                                                                                         >  4:151577/1‑147 (MQ=255)
                                              gTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGaaa                                                                                                         <  3:151577/147‑1 (MQ=255)
                                                  cGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCg                                                                                                   <  4:16090/149‑1 (MQ=255)
                                                   gAAACCACCGTGACCGTGTGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCgg                                                                                                  <  2:1480/149‑1 (MQ=255)
                                                   gAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCgg                                                                                                  >  3:70478/1‑149 (MQ=255)
                                                      accaccGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGc                                                                                               >  4:374897/1‑149 (MQ=255)
                                                       ccaccGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGT‑ACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCt                                                                                             >  8:287089/1‑149 (MQ=255)
                                                         accgtgaccgtgCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCtt                                                                                            <  4:210280/149‑1 (MQ=255)
                                                          ccgtgaccgtgCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCtt                                                                                            >  5:380227/1‑148 (MQ=255)
                                                                     cgcgGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCgg                                                                                 <  4:70478/148‑1 (MQ=255)
                                                                        ggCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAg                                                                             >  7:186480/1‑149 (MQ=255)
                                                                               gCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGc                                                                       >  1:200896/1‑148 (MQ=255)
                                                                               gCCGTGGCGGGCGTAACGCCGAGTTCCTGCTAAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGc                                                                       >  6:5146/1‑148 (MQ=255)
                                                                                ccGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCa                                                                     <  3:225147/149‑1 (MQ=255)
                                                                                      cgggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCg                                                               <  6:380227/149‑1 (MQ=255)
                                                                                       gggcgTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGa                                                              <  2:299709/149‑1 (MQ=255)
                                                                                           gTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAgc                                                                   >  4:112750/1‑140 (MQ=255)
                                                                                           gTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAgc                                                                   <  3:112750/140‑1 (MQ=255)
                                                                                           gTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGcc                                                          <  8:302486/149‑1 (MQ=255)
                                                                                             aaCGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCtt                                                        >  8:223680/1‑149 (MQ=255)
                                                                                              aCGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTg                                                       <  2:26843/149‑1 (MQ=255)
                                                                                               cGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTgg                                                      >  8:205632/1‑149 (MQ=255)
                                                                                                gCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTggg                                                     >  3:310207/1‑149 (MQ=255)
                                                                                                gCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTggg                                                     >  1:337105/1‑149 (MQ=255)
                                                                                                 ccGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGc                                                    <  1:365189/149‑1 (MQ=255)
                                                                                                 ccGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGc                                                    <  7:287089/149‑1 (MQ=255)
                                                                                                      ttCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTc                                               <  4:310207/149‑1 (MQ=255)
                                                                                                        ccTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCgg                                             >  3:158399/1‑149 (MQ=255)
                                                                                                         cTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCgg                                             <  5:5146/148‑1 (MQ=255)
                                                                                                         atGCTCAGCCTCACCGAAAGACTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCCCTGACACACACGGCAACGATCGCTCGGAAGAAAACGCCGCTCCCTTCAAGACCCCAGCCAGCTAAGCAAGAGCCGAAGCCTTGCCCCTGTCgcc                                            >  8:113475/2‑147 (MQ=255)
                                                                                                            cTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAg                                         <  3:374897/149‑1 (MQ=255)
                                                                                                             tCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGc                                        >  5:135838/1‑149 (MQ=255)
                                                                                                               atccTCACCGAAACCCTTAAAGGCATGACTGGCGTGTACGCCCTGGCCTGTTACACCGACTGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGcg                                       <  6:298814/146‑1 (MQ=255)
                                                                                                                 ccTCACCGAAAGCCTGAAAGGCCTGTCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGcgac                                     >  4:321964/1‑148 (MQ=255)
                                                                                                                 ccTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGcgac                                     >  7:85910/1‑148 (MQ=255)
                                                                                                                   tCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAg                                  >  4:191466/1‑149 (MQ=255)
                                                                                                                   tCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAg                                  >  5:85245/1‑149 (MQ=255)
                                                                                                                   tCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAg                                  >  4:78761/1‑149 (MQ=255)
                                                                                                                      ccGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                               <  7:205632/149‑1 (MQ=255)
                                                                                                                      ccGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                               <  1:347070/149‑1 (MQ=255)
                                                                                                                       cGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTggg                                                     >  2:127306/1‑126 (MQ=255)
                                                                                                                       cGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTggg                                                     <  1:127306/126‑1 (MQ=255)
                                                                                                                       cGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGcgacga                                   <  3:298984/144‑1 (MQ=255)
                                                                                                                       cGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGcgacga                                   >  4:298984/1‑144 (MQ=255)
                                                                                                                        gAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGa                             <  2:200896/149‑1 (MQ=255)
                                                                                                                         aaaGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGac                            <  3:252804/149‑1 (MQ=255)
                                                                                                                           aGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGacaa                          <  2:337105/149‑1 (MQ=255)
                                                                                                                              cTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGacaacaa                       >  6:189164/1‑149 (MQ=255)
                                                                                                                               tGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAg                                                                             >  4:328979/1‑94 (MQ=255)
                                                                                                                               tGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAg                                                                             <  3:328979/94‑1 (MQ=255)
                                                                                                                               tGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                               <  8:302035/140‑1 (MQ=255)
                                                                                                                               tGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                               >  7:302035/1‑140 (MQ=255)
                                                                                                                                 aaaGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTc                                               <  4:104190/122‑1 (MQ=255)
                                                                                                                                 aaaGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTc                                               >  3:104190/1‑122 (MQ=255)
                                                                                                                                     ccggccagggCGTACACGCCC‑GGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAacggct               <  3:237013/142‑1 (MQ=255)
                                                                                                                                     ccggccagggCGTACACGCCC‑GGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAacggc                >  4:237013/8‑148 (MQ=255)
                                                                                                                                       cTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGc                                            >  6:340932/1‑119 (MQ=255)
                                                                                                                                       cTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGc                                            <  5:340932/119‑1 (MQ=255)
                                                                                                                                       cTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAacggcta              >  2:96849/1‑149 (MQ=255)
                                                                                                                                          ccGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAacggctacgg           <  4:158399/149‑1 (MQ=255)
                                                                                                                                             ggCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAAcggctacggcta        <  3:191466/149‑1 (MQ=255)
                                                                                                                                              gCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACggctacggctac       <  5:189164/149‑1 (MQ=255)
                                                                                                                                                   tACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTcgc  <  8:85910/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTACTTCGC  >  NC_002947/4891727‑4892022

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑