Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,549,601 T→G 8.1% intergenic (+195/+45) 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,6010TG 8.1% 183.1 / 4.9 61intergenic (+195/+45)yfeW/yfeXpenicillin binding protein PBP4B; weak DD‑carboxypeptidase activity/porphyrinogen oxidase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base G (3/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CCGATCCGCAAAAGAATCCCAATATGTTCGAAAGCGGTCAGTTGCCGATTGCAACTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGG  >  NC_000913/2549303‑2549896
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                       
ccGATCCGCAAAAGAATCCCAATATGTTCGAAAGCGGTCAGTTGCCGATTGCAACTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCcagtcaa                                                                                                                                                                                                                                                                                                      >  2:97011/1‑300 (MQ=255)
              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)
                            cGAAAGCGGTCAGTTGCCGATTGCAACTTCTGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACACAAGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCAGGATAAGGCGTTCACGACGCATCCGGCGTGTACCACGCGCACCTGGTCAACAAACAGAGAACAGTCAGAAAtacaaaagccaggcagttaaccaa                                                                                                                                                                                                                                                                           >  2:337375/1‑278 (MQ=37)
                            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)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                       
CCGATCCGCAAAAGAATCCCAATATGTTCGAAAGCGGTCAGTTGCCGATTGCAACTTATGGTTGGGTAGTTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGTCCAGCGACGGTGCGAAATAATAGCCGCCGGTTACCGGTTTGGTGAAACGCAACATCGCATCACGCTTACCATCGGTATCGCCAAACATGCTCAGCAGTTGCTGCTCAATGTTATGCAGACGCGCGCAGTAGGCGCAGAAGTACAGACCGTGAGTGCCACTGGCAGTGCCGTACGGCAGGCTCTGGCGAACAATCTTCAGCCCTTTGCCATCTTCTTTCAGATCAACGCGGG  >  NC_000913/2549303‑2549896

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: