Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,891,873 G→A 6.9% V345V (GTG→GTA)  ttuD → hydroxypyruvate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,891,8730GA 6.9% 263.0 / 6.0 87V345V (GTG→GTA) ttuDhydroxypyruvate reductase
Reads supporting (aligned to +/- strand):  ref base G (37/44);  new base A (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 = 4.06e-01

TCGTTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTAC  >  NC_002947/4891725‑4892017
                                                                                                                                                    |                                                                                                                                                
tCGTTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCgtg                                                                                                                                                  >  3:35829/1‑149 (MQ=255)
  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                                               <  7:287089/149‑1 (MQ=255)
                                                                                                   ccGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGc                                               <  1:365189/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:78761/1‑149 (MQ=255)
                                                                                                                     tCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAg                             >  4:191466/1‑149 (MQ=255)
                                                                                                                     tCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAg                             >  5:85245/1‑149 (MQ=255)
                                                                                                                        ccGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                          <  1:347070/149‑1 (MQ=255)
                                                                                                                        ccGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                          <  7:205632/149‑1 (MQ=255)
                                                                                                                         cGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTggg                                                >  2:127306/1‑126 (MQ=255)
                                                                                                                         cGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTggg                                                <  1:127306/126‑1 (MQ=255)
                                                                                                                         cGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGcgacga                              >  4:298984/1‑144 (MQ=255)
                                                                                                                         cGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGcgacga                              <  3:298984/144‑1 (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                          >  7:302035/1‑140 (MQ=255)
                                                                                                                                 tGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTg                          <  8:302035/140‑1 (MQ=255)
                                                                                                                                   aaaGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTc                                          >  3:104190/1‑122 (MQ=255)
                                                                                                                                   aaaGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTc                                          <  4:104190/122‑1 (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)
                                                                                                                                                    |                                                                                                                                                
TCGTTCAACACGGCCAGCCGCTCAAAGCGCCCTGCGTGATCCTGTCGGGTGGCGAAACCACCGTGACCGTGCGCGGCAATGGCCGTGGCGGGCGTAACGCCGAGTTCCTGCTCAGCCTCACCGAAAGCCTGAAAGGCCTGCCGGGCGTGTACGCCCTGGCCGGTGACACCGACGGCATCGATGGCTCGGAAGAAAACGCCGGTGCCTTCATGACCCCGGCCAGCTACGCCAGCGCCGAAGCCTTGGGCCTGTCGGCCAGCGACGAGCTGGACAACAACAACGGCTACGGCTAC  >  NC_002947/4891725‑4892017

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: ‑