Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,970,817:1 +G 5.1% coding (685/1131 nt) wecB → UDP‑N‑acetyl glucosamine‑2‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,970,8171.G 5.1% 876.4 / 40.1 198coding (685/1131 nt)wecBUDP‑N‑acetyl glucosamine‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base . (95/93);  new base G (5/5);  total (100/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-02

Only 100 of 202 total aligned reads displayed.

GTATCACTTCTCTCCAACCGAAACTTCCCGGCAAAACTTGCTGCGTGAAAACGTTGCGGATAGCCGAATCTTCATTACCGGTAATACAGTCATTGATGCACTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTG  >  NC_000913/3970561‑3971083
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                          
gTATCACTTCTCTCCAACCGAAACTTCCCGGCAAAACTTGCTGCGTGAAAACGTTGCGGATAGCCGAATCTTCATTACCGGTAATACAGGCATTGATGCACTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGccacca                                                                                                                                                                                                                                                          <  1:36707/275‑1 (MQ=255)
 tATCACTTCTCTCCAACCGAAACTTCCCGGCAAAACTTGCTGCGTGAAAACGTTGCGGATAGCCGAATCTTCATTACCGGTAATACAGTCATTGATGCACTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGccacca                                                                                                                                                                                                                                                          <  4:69369/274‑1 (MQ=255)
     aCTTCTCTCCAACCGAAACTTCCCGGCAAAACTTGCTGCGTGAAAACGTTGCGGATAGCCGAATCTTCATTACCGGTAATACAGTCATTGATGCACTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAACAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTTCCAACGCCGTGCAGACATCGccacccccccc                                                                                                                                                                                                                                                     >  4:214969/1‑269 (MQ=255)
            tcCAACCGAAACTTCCCGGCAAAACTTGCTGCGTGAAAACGTTGCGGATAGCCGCATCTTCATTACCGGTAATACAGTCATTCCTGCCCTGTTATGGGTGCGTGACCAGGTGATGACCAGCGCCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATACAACGATGCTTCTGGTGACCGGTCACAGGCGTGAAGGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGAGCTGGAAGCCATCGCTACCCCGCCCCAGGAc                                                                                                                                                                                                                                               >  4:75904/1‑274 (MQ=255)
                 ccaaaaCATCCCGGAAAAACTTGCTTAGTTAAAACGTGTCGGATAGCCGAATCTTCAGTACCCGTAATACAGTCATTGATGGAATGTTAGGGGTGCGTGACCAGGTGCTGAGCAGCGACAAGCTGCGTGCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCCATAAAAAGATTATTCTGGTGACCGGGCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCTCCACCACGCACCAGGACATCCAg                                                                                                                                                                                                                                         <  4:519470/272‑1 (MQ=255)
                  cGAAACTTCCCTGCAAAACTTGCTGCGTTAAAACGTTGCGGATTGCCGAATCTTCATGACCGGTAATACCGTCATTGATGCGCTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGa                                                                                                                                                                                                                                        <  2:325670/275‑1 (MQ=255)
                                                                    tCTTCATTAACGGTAATACAGTCATTGATGCACTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCg                                                                                                                                                                                       <  2:75905/274‑1 (MQ=255)
                                                                         aTTACCGGTAATACAGTCATTGATGCACTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGccacca                                                                                                                                                                                                                                                          >  1:26172/1‑202 (MQ=255)
                                                                         aTTACCGGTAATACAGTCATTGATGCACTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCCCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCt                                                                                                                                                                                 >  4:505005/1‑275 (MQ=255)
                                                                             ccGGTAATAAAGTCCTAGATGCACGGTTATGGGGGCGTGACCAGGGGATGAGCAGCGACAAGCNGCGTTCNGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCa                                                                                                                                                                                                                                                      <  3:26172/202‑1 (MQ=255)
                                                                                gTAATACAGTCATTGATGCACTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAAAACCGTTCAATCGCATTCTGGGGCAt                                                                                                                                                                          >  2:23214/1‑275 (MQ=255)
                                                                                       aGTCATTGATGCACTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCCAATTACCCGTTTATCGACCCCGTTCAACAGATGATTCTGGTGACCGGTCCCAGGCGTGACAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCCCGCGCTGGCAGACCTCGCCACCTCGCACCAGGACCTCCCGATTGTCTATCCGGTGCATACCAACCCGAACCTCAGCGAACCGGTCACTCGCATTCCGCGGCCTGCGaaa                                                                                                                                                                    >  4:304421/1‑274 (MQ=255)
                                                                                               aTGCACTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGa                                                                                                                                                                                                                                        <  1:68701/198‑1 (MQ=255)
                                                                                                tggacgttaaagggTGTGTGCCCAGGTGAGGAGCAGCGACAAGCTGGGTTCATGACTGGCGGCAAATTACCCGTTTATCGACCGCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGGAAGTTTCGGTCGTGGCGTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCGCAACCCGAACGTCAGAGAACCGGTAAATCGCATTCGGGCGCATGTGAAAAATGTCAtt                                                                                                                                                           <  4:412343/263‑1 (MQ=255)
                                                                                                   aCTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGAATCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCAAGATTGTCCATCCGGTGACTCCCAACCCGAACGTCAAAGAACCGGGCAATCGCATTCTGGGGCATGGGAAAAATGTCATTCGga                                                                                                                                                       >  4:351203/1‑275 (MQ=255)
                                                                                                               gCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACAt                                                                                                                                                                                                                                             >  2:556657/1‑177 (MQ=255)
                                                                                                               gCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTTGGGCATGTGAAAAATGTCATTCTGATCGATCCCCCgg                                                                                                                                           >  1:97986/1‑275 (MQ=255)
                                                                                                                cGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTgg                                                                                                                                                                                                                                                                         >  6:763627/1‑148 (MQ=255)
                                                                                                                 gTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCTATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAACACTGTCATTCTTAACGATCCCCAGGat                                                                                                                                         >  2:8365/1‑274 (MQ=255)
                                                                                                                 gTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGGTAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAg                                                                                                                                         >  2:107218/1‑275 (MQ=255)
                                                                                                                 gTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                                                                                                                                          >  2:162839/1‑274 (MQ=255)
                                                                                                                 gTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCACTCTGGGGCATGTTAAACATGTCATTCTGACCGATCCCCAGGa                                                                                                                                          >  1:11342/1‑274 (MQ=255)
                                                                                                                   gACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAg                                                                                                                                                                                                                                         <  4:556656/177‑1 (MQ=255)
                                                                                                                           gatgaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGcc                                                                                                                                                                                                                                                              <  8:119078/148‑1 (MQ=255)
                                                                                                                                           cTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAgg                                                                                                                                                                                                                                                 <  8:1009642/145‑1 (MQ=255)
                                                                                                                                           cTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAgg                                                                                                                                                                                                                                                 >  7:1009642/1‑145 (MQ=255)
                                                                                                                                             gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGAc                                                                                                                                                                                                                                                                    <  6:183943/124‑1 (MQ=255)
                                                                                                                                             gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGAc                                                                                                                                                                                                                                                                    >  5:183943/1‑124 (MQ=255)
                                                                                                                                              cGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATcc                                                                                                                                                                                                                                           <  5:1074412/148‑1 (MQ=255)
                                                                                                                                                       cTGGCGGCAAATTACCCGTTGATCGACCCCGATAAAAAGAGGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTTGCTTTGAAGACATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTgg                                                                                                     <  2:467100/273‑1 (MQ=255)
                                                                                                                                                         ggcggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTg                                                                                                                                                                                                                                                                          >  5:963606/1‑106 (MQ=255)
                                                                                                                                                             gcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGa                                                                                                                                                                                                                                                                     >  7:317038/1‑107 (MQ=255)
                                                                                                                                                             gcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGa                                                                                                                                                                                                                                                                     <  8:317038/107‑1 (MQ=255)
                                                                                                                                                                aaTTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCa                                                                                                                                                                                                                                                      <  5:454142/119‑1 (MQ=255)
                                                                                                                                                                aaTTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCa                                                                                                                                                                                                                                                      >  6:454142/1‑119 (MQ=255)
                                                                                                                                                                     cccGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCa                                                                                                                                                                                                                                                   <  6:301678/117‑1 (MQ=255)
                                                                                                                                                                     cccGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCa                                                                                                                                                                                                                                                   >  5:301678/1‑117 (MQ=255)
                                                                                                                                                                       cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTg                                                                                                                                                                                                                          <  5:606509/140‑1 (MQ=255)
                                                                                                                                                                       cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTg                                                                                                                                                                                                                          >  6:606509/1‑140 (MQ=255)
                                                                                                                                                                       cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCaa                                                                                                                                                                                                                  >  8:1782/1‑148 (MQ=255)
                                                                                                                                                                         tttACCGACTCCGTTAAAAAGAGTATTGTGGTGACCGGCCACAGGCGGGGGAGTTTCGGTCGTGGCTTGGATGAAAGCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATTAAGATAGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTATGAGCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAgg                                                                                 <  4:107829/275‑1 (MQ=255)
                                                                                                                                                                          ttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGGGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAgg                                                                                 <  2:351204/274‑1 (MQ=255)
                                                                                                                                                                          ttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCa                                                                                                                                                                                                                                                   >  7:403931/1‑112 (MQ=255)
                                                                                                                                                                          ttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCa                                                                                                                                                                                                                                                   >  7:803321/1‑112 (MQ=255)
                                                                                                                                                                          ttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCa                                                                                                                                                                                                                                                   <  8:803321/112‑1 (MQ=255)
                                                                                                                                                                                aCCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGccacc                                                                                                                                                                                                                                                           <  7:56312/98‑1 (MQ=255)
                                                                                                                                                                                aCCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGccacc                                                                                                                                                                                                                                                           >  8:56312/1‑98 (MQ=255)
                                                                                                                                                                                   ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACa                                                                                                                                                                                                                                              >  8:1036493/1‑108 (MQ=255)
                                                                                                                                                                                   ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                                                                                                                                                     >  8:946646/1‑117 (MQ=255)
                                                                                                                                                                                   ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGttt                                                                                                                             >  2:212033/1‑221 (MQ=255)
                                                                                                                                                                                   ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACCCACCAGGACa                                                                                                                                                                                                                                              <  7:1036493/108‑1 (MQ=255)
                                                                                                                                                                                         aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACAt                                                                                                                                                                                                                                                                  <  6:381735/82‑1 (MQ=255)
                                                                                                                                                                                         aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACAt                                                                                                                                                                                                                                                                  >  5:381735/1‑82 (MQ=255)
                                                                                                                                                                                           aaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGAt                                                                                                                                                                                                                                       <  5:297348/107‑1 (MQ=255)
                                                                                                                                                                                           aaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCa                                                                                                                                                                                                                   >  8:108778/1‑127 (MQ=255)
                                                                                                                                                                                           aaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATAGCCACCACGCACCAGGACATCCAGAt                                                                                                                                                                                                                                       >  6:297348/1‑107 (MQ=255)
                                                                                                                                                                                              atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATc                                                                                                                                                                                                                               >  7:185349/1‑112 (MQ=255)
                                                                                                                                                                                              atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCa                                                                                                                                                                                           >  6:319124/1‑148 (MQ=255)
                                                                                                                                                                                              atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCa                                                                                                                                                                                           >  7:952005/1‑148 (MQ=255)
                                                                                                                                                                                                   tCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTTAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATCGTCTATCCGGTGCATCTCAACCCGACCGTCAGAGAACCGGTCAATAGCATTCAGGGGCATGTGAAAAATGTAATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCGGGCTGATGAACCACGCCTGGCTGATTTTGGATGACCCAGGCGGCATTCAGGAAGCAGCGCCTTCGc                                                       >  4:103875/1‑275 (MQ=255)
                                                                                                                                                                                                   tCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                                                                                                                                                      >  7:114782/1‑148 (MQ=255)
                                                                                                                                                                                                       gTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGGGAAAAATGTCATTCTg                                                                                                                                                        >  4:455217/1‑174 (MQ=255)
                                                                                                                                                                                                         gACCGGTCACAGGCGTGAGAGTTTCGGTCGTAGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCt                                                                                                                                                                                 >  6:1059260/1‑147 (MQ=255)
                                                                                                                                                                                                           ccGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatcg                                                                                                                                                    <  2:455218/174‑1 (MQ=255)
                                                                                                                                                                                                              gTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACAt                                                                                                                                                                                                                                             >  2:294813/1‑82 (MQ=255)
                                                                                                                                                                                                               taacgggCGTGGGTGTTTCNNTCGTGNTGTTGAAAACACCGGCCCCGCGCTGGCAGACACCGCCCCCACTCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACGCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGTTCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAAc                                              <  3:11342/267‑1 (MQ=255)
                                                                                                                                                                                                               tCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGt                                                                                                                                                                                                         <  5:631509/117‑1 (MQ=255)
                                                                                                                                                                                                               tCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGt                                                                                                                                                                                                         >  6:631509/1‑117 (MQ=255)
                                                                                                                                                                                                                  caGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAg                                                                                                                                                                                                                                         <  4:294812/82‑1 (MQ=255)
                                                                                                                                                                                                                        tGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                                                                                                                                                        <  6:729713/141‑1 (MQ=255)
                                                                                                                                                                                                                        tGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                                                                                                                                                        >  5:729713/1‑141 (MQ=255)
                                                                                                                                                                                                                             gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaaa                                                                                                                                                                  <  5:234245/142‑1 (MQ=255)
                                                                                                                                                                                                                             gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaaa                                                                                                                                                                  >  6:234245/1‑142 (MQ=255)
                                                                                                                                                                                                                             gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCAt                                                                                                                                                            <  5:319124/148‑1 (MQ=255)
                                                                                                                                                                                                                             gTTTCGGTCCGGGCTTTGAAGAAATCGGCCACCGGCTGGCAGACATCGCCAGCACGCACCAGGACAACCAGATGGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATgcg                              <  2:517397/274‑1 (MQ=255)
                                                                                                                                                                                                                               ttCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatcga                                                                                                                                                   >  4:272809/1‑155 (MQ=255)
                                                                                                                                                                                                                                         tcttctaagaaATCTGCCACGCGCGGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAACCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCGGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGa                  <  4:162838/268‑1 (MQ=255)
                                                                                                                                                                                                                                            ttGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCTCCa                                                                                                                                             >  8:647307/1‑148 (MQ=255)
                                                                                                                                                                                                                                               aagaaATCTGTCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                                                                                                                                          >  5:837971/1‑148 (MQ=255)
                                                                                                                                                                                                                                               aagaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAc                                                                                                                                                                                                                 >  7:209656/1‑77 (MQ=255)
                                                                                                                                                                                                                                               aagaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAc                                                                                                                                                                                                                 <  8:209656/77‑1 (MQ=255)
                                                                                                                                                                                                                                                agaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCa                                                                                                                                                                           <  6:1126709/114‑1 (MQ=255)
                                                                                                                                                                                                                                                agaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAg                                                                                                                                         >  8:827532/1‑148 (MQ=255)
                                                                                                                                                                                                                                                agaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAg                                                                                                                                         >  5:1127746/1‑148 (MQ=255)
                                                                                                                                                                                                                                                 gaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                                                                                                                                                      <  7:510192/102‑1 (MQ=255)
                                                                                                                                                                                                                                                 gaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                                                                                                                                                      >  8:510192/1‑102 (MQ=255)
                                                                                                                                                                                                                                                 gaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGt                                                                                                                                        >  8:403755/1‑148 (MQ=255)
                                                                                                                                                                                                                                                    aTCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtg                                                                                                                                                                         >  6:925456/1‑112 (MQ=255)
                                                                                                                                                                                                                                                    aTCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAtt                                                                                                                                     <  7:911012/148‑1 (MQ=255)
                                                                                                                                                                                                                                                          cACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGt                                                                                                                               >  7:1076010/1‑148 (MQ=255)
                                                                                                                                                                                                                                                          cACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGt                                                                                                                               >  5:784876/1‑148 (MQ=255)
                                                                                                                                                                                                                                                          cACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCACCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATCCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGt                                                                                                                               <  7:647307/148‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGtt                                                                                                                              >  8:907391/1‑148 (MQ=255)
                                                                                                                                                                                                                                                             cgcgCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCaa                                                                                                                                                                                          >  6:359774/1‑87 (MQ=255)
                                                                                                                                                                                                                                                             cgcgCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCaa                                                                                                                                                                                          <  5:359774/87‑1 (MQ=255)
                                                                                                                                                                                                                                                             cgcgCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAttt                                                                                                                                    <  5:1025974/141‑1 (MQ=255)
                                                                                                                                                                                                                                                              gcgcTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCt                                                      <  1:4660/218‑1 (MQ=255)
                                                                                                                                                                                                                                                               cgcTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTga                                                                                                                                                       <  8:596624/120‑1 (MQ=255)
                                                                                                                                                                                                                                                               cgcTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTga                                                                                                                                                       >  7:596624/1‑120 (MQ=255)
                                                                                                                                                                                                                                                               cgcTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGt                                                                                                                           <  7:827532/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                          
GTATCACTTCTCTCCAACCGAAACTTCCCGGCAAAACTTGCTGCGTGAAAACGTTGCGGATAGCCGAATCTTCATTACCGGTAATACAGTCATTGATGCACTGTTATGGGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTG  >  NC_000913/3970561‑3971083

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: