Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,549,609 G→A 8.0% intergenic (+203/+37) yfeW → / ← yfeX penicillin binding protein PBP4B; weak DD‑carboxypeptidase activity/porphyrinogen oxidase, cytoplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,549,6090GA 8.0% 180.7 / 4.4 62intergenic (+203/+37)yfeW/yfeXpenicillin binding protein PBP4B; weak DD‑carboxypeptidase activity/porphyrinogen oxidase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base A (3/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

AATCCCAATATGTTCGAAAGCGGTCAGTTGCCGATTGCAACTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTG  >  NC_000913/2549317‑2549903
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                      
aaTCCCTATGTGTTCGAATGCGGTAGGTTGCCGATTGAAACTTATGGTTTGGGTGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGcc                                                                                                                                                                                                                                                                                                <  1:402356/301‑1 (MQ=255)
    cgttcttGTTCGAAAGCGGTCAGTTGCCGATTGCAACTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGGCAGTAACATTACTGACTGGc                                                                                                                                                                                                                                                                                            <  2:310779/295‑1 (MQ=255)
     cAATATGTTCGAAAGCGGTCAGTTGCCGATTGCAACTTATGGTTGGGTAGGTGAGCAGGTGTATGGGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGGCAGTAACATTACTGACTGGcc                                                                                                                                                                                                                                                                                           <  2:569296/301‑1 (MQ=255)
          tGTTCGAAAGCGGTCAGTTGCCGATTGCAACTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTAAGGCCAAGCAGTAACATTACTGACTGGcagaacg                                                                                                                                                                                                                                                                                      >  1:310375/1‑295 (MQ=255)
              cGAAAGCGGTCAGTTGCCGATTGCAACTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCaa                                                                                                                                                                                                                                                                                  >  2:406901/1‑301 (MQ=255)
                aaaGCGGTCAGTTGCCGATTGCAACTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATg                                                                                                                                                                                                                                                                                <  2:107158/301‑1 (MQ=255)
                     ggTCAGTTGCCGATTGCAACTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCa                                                                                                                                                                                                                                                                           <  2:147324/301‑1 (MQ=255)
                        aaGTTGCCGATTGCAATTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGaaa                                                                                                                                                                                                                                                                       <  1:96887/300‑1 (MQ=255)
                             gCCGATTGCAACTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCt                                                                                                                                                                                                                                                                   >  1:456782/1‑301 (MQ=255)
                                     cAACTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCAATTGGTCAACAATATGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTAACAgcgc                                                                                                                                                                                                                                                           >  1:340626/1‑301 (MQ=255)
                                       atttttGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCa                                                                                                                                                                                                                                                         <  2:639551/296‑1 (MQ=255)
                                          tATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCa                                                                                                                                                                                                                                                      >  2:171283/1‑301 (MQ=255)
                                             gattGGGTAGTTGGTCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACt                                                                                                                                                                                                                                                   <  1:55501/299‑1 (MQ=255)
                                                 gggTAGTTGTTCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTc                                                                                                                                                                                                                                               <  1:171086/301‑1 (MQ=255)
                                                  ggTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACtagatc                                                                                                                                                                                                                                              >  2:55569/1‑296 (MQ=255)
                                                              ggTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGGGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATAAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCg                                                                                                                                                                                                                                  >  2:642197/1‑301 (MQ=255)
                                                               gtgtATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGATTACAGCGCCATCAACTTGTCCAGCAACGGTGCgc                                                                                                                                                                                                                                 >  1:288102/1‑300 (MQ=255)
                                                                     gcggcgTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAaataaa                                                                                                                                                                                                                           >  2:168090/1‑300 (MQ=255)
                                                                              aaGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCgg                                                                                                                                                                                                                  <  2:457346/301‑1 (MQ=255)
                                                                                     aGTAATACAAAGGGCATTCAGTAATTTTATGGGCTGGTCTGAGTGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACTTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTTCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCgg                                                                                                                                                                                                           <  1:535276/301‑1 (MQ=255)
                                                                                                 gCCATTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACg                                                                                                                                                                                               <  2:341089/301‑1 (MQ=255)
                                                                                                       cagtAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACaa                                                                                                                                                                                         >  2:138050/1‑300 (MQ=255)
                                                                                                                     ggtgggtggagaGTGCTGGCAAACGCAAAACGTCCTGTGTCGCTACGCTTTTCAGGCCTACGTGGAGCCCACAATTTATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTAc                                                                                                                                                                           <  1:406338/293‑1 (MQ=255)
                                                                                                                      ctggTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTTCCAAATAATAGCCGCCAGATACCGGTTTGGTGAAACCCAACATCGCATCACGCTTAcc                                                                                                                                                                          >  1:414180/1‑301 (MQ=255)
                                                                                                                           ttGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCgg                                                                                                                                                                     <  2:639247/300‑1 (MQ=255)
                                                                                                                                        cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACAt                                                                                                                                                        >  1:120480/1‑301 (MQ=255)
                                                                                                                                        cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCACCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACaa                                                                                                                                                        >  2:472302/1‑300 (MQ=255)
                                                                                                                                        cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACaa                                                                                                                                                        >  2:251791/1‑300 (MQ=255)
                                                                                                                                          aaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGc                                                                                                                                                      <  2:288504/301‑1 (MQ=255)
                                                                                                                                           aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGcc                                                                                                                                                     >  2:544810/1‑300 (MQ=255)
                                                                                                                                                         gATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCACCGGTAACGCCAAACATGCCCAGCAGTtcctgcc                                                                                                                                       >  1:72721/1‑300 (MQ=255)
                                                                                                                                                             ggttACGCTTATCAGTCCTACGTGGTGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAAt                                                                                                                                   <  1:251524/298‑1 (MQ=255)
                                                                                                                                                             cgcTATGCTTATCAGGCCTACGTGGAGCCCACAAGATATTGAATTAGCACGATTTCGTAGGCCGGAGAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGAACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTGGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAAt                                                                                                                                   <  1:544137/301‑1 (MQ=255)
                                                                                                                                                             cgcTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATAGCATCACGCTTACCATCGGTATCGCCAAACATGCTAAGCAGTTGCTGCTCaaa                                                                                                                                   >  2:76514/1‑300 (MQ=255)
                                                                                                                                                                 aCGCTTATCAGGCTTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTa                                                                                                                               <  2:120604/301‑1 (MQ=255)
                                                                                                                                                                  cGCTTATCAGGCCTACGTGGTGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTAt                                                                                                                              <  2:414743/301‑1 (MQ=255)
                                                                                                                                                                                  ggggtgCCCACAATATTTTGCATTATCACGGTTTTGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTa                                                                                                              <  1:137925/296‑1 (MQ=255)
                                                                                                                                                                                    ggAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAgg                                                                                                            >  2:57002/1‑301 (MQ=255)
                                                                                                                                                                                       tcccgtccctgcgcttttccgatctttttCGTAGGCCGGATAATGCTTTCACGCCGCATCCGGCGTGTACAACTCGCACTTTGTCAACAATCTGAGGCCATTCAGTAATGTTACTGCCTGTCCTTTAATCAATGAAGCAGAAACGCTTACAGCGCCATCAACTTTTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGgcgc                                                                                                         <  1:468781/275‑1 (MQ=255)
                                                                                                                                                                                          tctcttccgatctaTTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGa                                                                                                      <  2:156116/288‑1 (MQ=255)
                                                                                                                                                                                          cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGa                                                                                                      <  2:72811/301‑1 (MQ=255)
                                                                                                                                                                                                tatttGAATTTTCACGGTTTCGTTGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACa                                                                                                <  1:641320/298‑1 (MQ=255)
                                                                                                                                                                                                      aaTTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCCATGAATCAGAACCGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCCCGCTTACCATCGGTATCGCCAACCATGCTCAGCAGTTGCTGCTCAATGTTCAGCAGACGCGCACAGTAAGCGCCGAAGTACAGAccag                                                                                          >  1:79385/1‑299 (MQ=255)
                                                                                                                                                                                                       aTTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCGGCTCAATGTTATGCAGACCCGCGCAGTAGGCGCAGAAGAACGGaagagg                                                                                         >  1:155965/1‑296 (MQ=255)
                                                                                                                                                                                                        ttAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCCCAATGTTATGCAGACGCCCGCAGGAGGCGCCGAAGAACAACCCGTgg                                                                                        >  1:221384/1‑300 (MQ=255)
                                                                                                                                                                                                           gcccgttttCTTAGGCCGGTTAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTAATGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTg                                                                                     <  1:167893/295‑1 (MQ=255)
                                                                                                                                                                                                                aTTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGatcggaagagaacacgactgaaa                                                                                >  2:469382/1‑279 (MQ=255)
                                                                                                                                                                                                                    cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAAACCACGAGTGCACCTGGcc                                                                            >  1:411400/1‑300 (MQ=255)
                                                                                                                                                                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACTGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTACCACTGGCAGTGcc                                                                       >  2:446939/1‑301 (MQ=255)
                                                                                                                                                                                                                               tAAGGCGTTCAAGCCTCTTCCGGCGTGTACAACGTGAACTTGGTCAACAATATGAGGCAGGTCAGTAAGGTTACTGCCTGGCCTGTAATCAATGAATCAGAAACGCTTACAGCGCCCTCAACTTGTGCAGGGACGGTGCGAAATAAGAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACgg                                                                 <  1:622479/301‑1 (MQ=255)
                                                                                                                                                                                                                                 aGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCa                                                               >  2:68845/1‑301 (MQ=255)
                                                                                                                                                                                                                                          cgccgcTTCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGc                                                      <  1:68755/301‑1 (MQ=255)
                                                                                                                                                                                                                                                 tCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATTTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAAACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCaaaaaaa                                               >  2:372299/1‑300 (MQ=255)
                                                                                                                                                                                                                                                     gCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTc                                           <  2:221627/301‑1 (MQ=255)
                                                                                                                                                                                                                                                          tACAACGCGCTCTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCTCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGccc                                      >  2:88613/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                 ctctcTTGGTCAACAATCTGAGTCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCACTGACTCAGAACCGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGGCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCa                               <  2:501126/297‑1 (MQ=255)
                                                                                                                                                                                                                                                                             acaaTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCACATGGCAGTTCCGTACGGCAAGCTCTGGCGAACAAACTTCAGCCCTTTGCCCTCTTCtttacga                   >  1:243228/1‑297 (MQ=255)
                                                                                                                                                                                                                                                                                aTCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCa                >  2:55161/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                    aaGGCCAGGCAGTAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAAcgc            >  2:421506/2‑301 (MQ=255)
                                                                                                                                                                                                                                                                                       gCCAGGCAGTAACATTACTGACTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCggg         >  2:156800/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                             cagtAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGggtgagg  >  2:565805/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                              agtAATGTTACTGTCTGGCGTTTTATCAATGAATCAGAAACGGTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACTGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGggtgaggtg  <  1:446375/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                      
AATCCCAATATGTTCGAAAGCGGTCAGTTGCCGATTGCAACTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGGTGAGGTG  >  NC_000913/2549317‑2549903

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