Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,399,929 A→T L17Q (CTG→CAG)  mreC ← cell shape determining protein MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,399,9290AT97.8% 308.4 / ‑1.5 90L17Q (CTG→CAG) mreCcell shape determining protein MreC
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (50/38);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

TCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACG  >  NC_000913/3399694‑3400175
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       
tCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCt                                                                                                                                                                                                                                          >  2:124712/1‑251 (MQ=255)
   aaGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAACGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCCATAATAATGCCGAGCGCCACCAGCCCCGCCTGAATAAGGCGAATTTGt                                                                                                                                                                                                                                        >  2:209068/1‑250 (MQ=255)
      cATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCg                                                                                                                                                                                                                                    >  2:152303/1‑251 (MQ=255)
         cagcagTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGAc                                                                                                                                                                                                                                  <  1:256463/250‑1 (MQ=255)
             agTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGcc                                                                                                                                                                                                                            >  2:69994/1‑251 (MQ=255)
                tCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACg                                                                                                                                                                                                                         >  1:7475/1‑251 (MQ=255)
                 cACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACgg                                                                                                                                                                                                                        <  1:530879/251‑1 (MQ=255)
                      tttttCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAAGCGAATCTGTAGCGACGGGCCACGGCTaa                                                                                                                                                                                                                    >  1:302463/1‑250 (MQ=255)
                         ttcggCAACAGTTACTGAAGTAACGCCCGGTTTTCAAGTTCTACTTGGTCACGCGAGGCAAGCGTCTGCGATACGCCAGCCAGCAATTCACGAGGAGCATTGGAAACAAAGAAGAAAGGACTGACGACGGTATCCAAATAAGTACGGATTTGACTGAAAGACCCCAGGCGACTGTCGGCAAAAATAATGCCGGTCGCCCCCAGCCCCGCCTGAATAAGAGGAATCTGTAgggccggcccaaagataacatt                                                                                                                                                                                                                 >  2:147238/5‑230 (MQ=255)
                         ttaaGCAACAGGTCCTGAAGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCGGCGATACGACATCCAGCAATTCACGCGGCGCATTGGAAACAAAGGAGAAAGGACTGACGGCGGTATCCATATAAGGACGGATTTGACTGAACGGCCCCAGGAGGATGTCGGCAATAATAAAGCCGAGCGACACCAGCACCACCTGAATAAGGCGAATCTGttgggagaagccacggataaaaa                                                                                                                                                                                                                  >  2:147236/4‑228 (MQ=255)
                           cAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGGTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCCGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGCCGGGCCACGGCTAAAACTTg                                                                                                                                                                                                              >  1:421017/1‑251 (MQ=255)
                              cAACAGATCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATACTAATGCCGAGCGCCACCAGCACCGCCTGAATACGGCGAATCTGTAGCGACGTGCCCCGGCTAAAAATTGGCt                                                                                                                                                                                                           >  2:148319/1‑251 (MQ=255)
                               accAGTGCCTGACGTAACGCCCGGTTTCCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGTGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTAACCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCtt                                                                                                                                                                                                          <  1:581055/249‑1 (MQ=255)
                                  aGTTCCTGACGAAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGCAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCCGTAGCGACGGGCCCCGGCTAAAAATTGGCTTCAt                                                                                                                                                                                                       >  2:272316/1‑251 (MQ=255)
                                      ccTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCCGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGCCGGGCCCCGGCTAAAAATTGGCTTCATAAGt                                                                                                                                                                                                   >  2:576582/1‑251 (MQ=255)
                                       cTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGtt                                                                                                                                                                                                  >  1:462359/1‑251 (MQ=255)
                                       cTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGtt                                                                                                                                                                                                  >  1:83051/1‑251 (MQ=255)
                                         gACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTAt                                                                                                                                                                                                <  1:219824/251‑1 (MQ=255)
                                         gACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTAt                                                                                                                                                                                                >  1:115519/1‑251 (MQ=255)
                                            gtaccGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCAGCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCg                                                                                                                                                                                             <  1:25203/247‑1 (MQ=255)
                                            gTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGCATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCg                                                                                                                                                                                             >  2:237978/1‑251 (MQ=255)
                                                 gCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGttatc                                                                                                                                                                                        >  2:454446/1‑247 (MQ=255)
                                                    cGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAACGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGCCGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCg                                                                                                                                                                                     >  2:220989/1‑251 (MQ=255)
                                                     ggTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGACTGAATAAGGCGAATCTGTAGCGACGGGCCCCGGCTAAAAATTGGCTTCCTAAGTTATGCGTATTCTCGt                                                                                                                                                                                    >  2:570604/1‑251 (MQ=255)
                                                       ttttCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCTTATAAGTCGGGATTGAACTGAACGTCCCCAGGCGGCTGTCGGCAAAAAAAATGCCGACGGCCAACAGCACCGCCTGAATAAGGCGAATTTTAACGACCGGGCCACGCCTAAAAATTGGCTTCAAAAGttctcgttttttccgtat                                                                                                                                                                                  >  2:427845/1‑235 (MQ=255)
                                                       ttttCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTAt                                                                                                                                                                                  >  1:509880/1‑251 (MQ=255)
                                                        tttCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTAt                                                                                                                                                                                  <  1:338416/250‑1 (MQ=255)
                                                          tCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTTAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCa                                                                                                                                                                                <  1:64085/250‑1 (MQ=255)
                                                            aaGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGAc                                                                                                                                                                             >  2:105121/1‑251 (MQ=255)
                                                                   aaTTGGTCACGCGAGGCCAGCGTCGGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACcagg                                                                                                                                                                         <  2:132846/248‑1 (MQ=255)
                                                                    aTTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACcaggcagg                                                                                                                                                                     >  1:177636/1‑251 (MQ=255)
                                                                       ggTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCTTATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGt                                                                                                                                                                   <  1:124712/250‑1 (MQ=255)
                                                                        gTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATACCGAGCGCCACCAGCACCGCCTGAATAAGTCGAATCTGTAGCGCCGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGGATTCTCGTATCAGACCAGGCCgggga                                                                                                                                                                  >  2:198601/1‑248 (MQ=255)
                                                                          accgcgAGGCCAGCGTCGCCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTGCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTaaa                                                                                                                                                                <  2:83051/248‑1 (MQ=255)
                                                                              cgAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGac                                                                                                                                                           >  2:130711/1‑251 (MQ=255)
                                                                               gAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAACTTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGac                                                                                                                                                           >  1:581120/1‑250 (MQ=255)
                                                                               gAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAAACCAGGCAGGGTAAAAAGac                                                                                                                                                           >  2:378594/1‑250 (MQ=255)
                                                                               gAGGCCAGAGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGACGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGACCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAgcc                                                                                                                                                           >  1:494226/1‑248 (MQ=255)
                                                                                  gCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACtt                                                                                                                                                       <  1:4855/251‑1 (MQ=255)
                                                                                    cAGCGTCTGCGACACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTcc                                                                                                                                                     >  2:528084/1‑251 (MQ=255)
                                                                                     aGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTccc                                                                                                                                                    >  1:347684/1‑251 (MQ=255)
                                                                                        gTCTGAGATACGCCATCCAGCAATTCACGAGGAGCATTGGAACCAACGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGCTTTGACTAAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGGGAGCCACCAGCCCCGCCTGAATACGGCGAATCTGTAGCGACGGGCCCCGGCTAAAAATTGGCTTCATAAGTTATGCGTCTTCTCGTAACCAACCAGGCAGGGTAAACAGACACTTACcctt                                                                                                                                                 >  2:342018/1‑250 (MQ=255)
                                                                                          cTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCcctgcc                                                                                                                                               >  2:506832/1‑251 (MQ=255)
                                                                                           tGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCCCCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCGTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGGAACCATCCCCTTCCcctgcct                                                                                                                                              >  2:182602/1‑251 (MQ=255)
                                                                                                  cGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCg                                                                                                                                       <  1:54084/251‑1 (MQ=255)
                                                                                                        ccAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTAct                                                                                                                                 >  1:89761/1‑251 (MQ=255)
                                                                                                        cacgaaatccACGAGGAGAATTGGAAACAAAGTAGAAAGGACGAAGGGCGGTTTCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGCCGGCAATAATAATGCCGAGGGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTAct                                                                                                                                 <  2:464022/242‑1 (MQ=255)
                                                                                                          atcAATCCACGAGGAGAATGGGAAACAAAGTAGAAAGGACGGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGAAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTActct                                                                                                                               <  2:89761/249‑1 (MQ=255)
                                                                                                               ttCACGAGGAGCATTGGAAACAAAGTAGAACGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATATGTAGCGACGGGCCACGGCTCAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAAACCAGGCAGGGTAAACCGACACTTCCCCTGCCTGCATCCGATTACTCTTCgat                                                                                                                          >  2:125529/1‑249 (MQ=255)
                                                                                                                  aCGAGGAGCATTGGAAACAAAGTAGAAAGGACTGAAGGCGGTATCCATATAAGTAAGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGAAACTTCCCCTGCCTGAATCCGATTACTCTTCGCTgc                                                                                                                        >  1:418115/1‑249 (MQ=255)
                                                                                                                    aaggacCATTGAAAACAAATTAGAAAGGACTGACGCCGGTATCAATATAAGTACGCATTTGAATCAACGTCCCCAGGCGGTTTTCGTCAATAATAATACCAATCGCAACCAGCACCGCCTGAATAAGGCAAATCTTTATCTACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCTGACCAGCCAGGTTAAACAAACACTTCCCCTGCCTGCATCCGATTACTCTTCTCTGAACa                                                                                                                     <  2:255647/250‑1 (MQ=255)
                                                                                                                              ggAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCgccgcc                                                                                                           >  1:569564/1‑251 (MQ=255)
                                                                                                                               gAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGccgccg                                                                                                          <  1:59344/251‑1 (MQ=255)
                                                                                                                                aaaCAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGACCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCACACCAGGCAGGGTAAACAGAAACTTCCCCTGCCTGCATCCGATGCCTCTTCGCTGAACAGGTCCCCGCCGt                                                                                                         >  2:391899/1‑251 (MQ=255)
                                                                                                                                 aaCAAAGTAGAAAGGACTGACGGCGGTATCCATATACGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTg                                                                                                        >  2:180327/1‑251 (MQ=255)
                                                                                                                                  aCAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGACCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATACGGCGAATCTGTAGCGACGGGCCACGGCTAACAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGAAACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGc                                                                                                       >  1:531212/1‑251 (MQ=255)
                                                                                                                                    aaaGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCCGGCAGGGTAAACAGACACTTCCACTGCCTGCCTCCGCTTACTCTTCGCTGAACAGGTCGCCGCCGTGCAt                                                                                                     >  2:575688/1‑251 (MQ=255)
                                                                                                                                      aGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGt                                                                                                   <  1:225630/251‑1 (MQ=255)
                                                                                                                                          aaaaGGACGGACGGCGGTTCCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAAGCCCGAGCCCCCCCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCGCGTATCAGACCAGGCAGGGTAAACAGTCACTTCCCCTGCCTGCAGCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGAt                                                                                               <  2:28313/250‑1 (MQ=255)
                                                                                                                                            aaGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATc                                                                                              >  1:142382/1‑250 (MQ=255)
                                                                                                                                             aGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCAt                                                                                            >  2:347749/1‑251 (MQ=255)
                                                                                                                                                   gccggcggTATCCATATAAGTACGGATTCGACTGATCGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGATCGCCCCCAGCACCGCCTGAATAAGGCGAATCTGTAGCGTCGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGAGCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCAGCCGATTACTCTTCGCTGAACAGGTCGCGGCCGTGCATGTCGATCATTTCCAg                                                                                      <  1:299343/249‑1 (MQ=255)
                                                                                                                                                   gACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTcc                                                                                        <  1:392587/249‑1 (MQ=255)
                                                                                                                                                        cggTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCtt                                                                                 <  1:378594/251‑1 (MQ=255)
                                                                                                                                                         ggTAGCCATATAAGTACGGATTTGACGGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACCCTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCttt                                                                                <  1:209068/251‑1 (MQ=255)
                                                                                                                                                             accATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCGGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgccg                                                                            <  2:115519/250‑1 (MQ=255)
                                                                                                                                                              ccATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgccgc                                                                           >  2:55495/1‑251 (MQ=255)
                                                                                                                                                               cATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTgccgcc                                                                          >  2:368802/1‑251 (MQ=255)
                                                                                                                                                                aaaaaaGTACGGATTTGACGGAACGCCCCCAGGCGGCTGTCGGCATGAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAACTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCa                                                                         <  2:7475/247‑1 (MQ=255)
                                                                                                                                                                         cGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGcc                                                                >  1:83935/1‑251 (MQ=255)
                                                                                                                                                                                 aCTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGt                                                        <  1:152303/251‑1 (MQ=255)
                                                                                                                                                                                      ccgcccccAGGGGGCGGTGGGCAAAAATAAGCCCGAGCGCACCAAGCCCCGCCTGAATAAGGCGAATTGGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCg                                                   <  2:138147/247‑1 (MQ=255)
                                                                                                                                                                                        gTCCCCAGGCGGCTGTCGGCAATAAGAATGCCGAGCTCCACCACCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGACTAAAAATTGTCTCCATAAGGTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCTCTGCCTGCATCCGCTTACTCTTCGCTTAACAGGTCGCCTCCGTTCATGTCGATCATTTCCAGTTCTTTGCCGCCACCGCGCGCCACACAGGTCAGCggg                                                 <  1:43647/251‑1 (MQ=255)
                                                                                                                                                                                           cccAGGCGGCGGTCGGAAATAATAATGCCGGGGGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCt                                              <  2:142382/251‑1 (MQ=255)
                                                                                                                                                                                              aGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCa                                           <  1:69994/251‑1 (MQ=255)
                                                                                                                                                                                               ggcggcTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGAt                                                                                                                                     >  2:386761/1‑160 (MQ=255)
                                                                                                                                                                                               ggcggcTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGAt                                                                                                                                     <  1:386761/160‑1 (MQ=255)
                                                                                                                                                                                               ggcggcTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAg                                          <  1:184899/251‑1 (MQ=255)
                                                                                                                                                                                                  ggcTGTGGTCAATAATAATGCCGTGCGCCACCATCACCGCCTGAATAAGGCGAATCTGTAGCTACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTCCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGcaa                                       <  1:320695/251‑1 (MQ=255)
                                                                                                                                                                                                          gcaataatatctcctatctccaccaTCACCGCCTGAATAAGGCGAATCCGTAGGGACGGGCCACTGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCTCCACACAGGTCAGCGGGTCTTCAGCAACAACGACt                               <  1:147238/232‑1 (MQ=255)
                                                                                                                                                                                                          gCAATATTAATGCCGAGCGCCACCAGCACCGCCTTAATAATGCGAATCTGTAGCGACGGGCCTCGGCTAAAAATTGGCTTCATAATTTATGCGTATTCTCGTATCAGACCATGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCATCGCTTTGCCGCCACCGCGCGCCACACATGTCAGCGGGTCTTCAGcaacatctact                               <  1:147236/251‑7 (MQ=255)
                                                                                                                                                                                                          gCAATAAAAATGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACt                               <  1:570604/251‑1 (MQ=255)
                                                                                                                                                                                                            aataataatGCCGAGCGCCCCCAGCACCGCCTGCATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGCTTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGCTCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGGCGGCGGGGCTTCAGCAACAACGCCTgg                             >  2:294740/1‑251 (MQ=255)
                                                                                                                                                                                                              taataatGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGa                            >  2:461066/1‑250 (MQ=255)
                                                                                                                                                                                                              aaataatGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACTCTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGaa                           <  1:105121/250‑1 (MQ=255)
                                                                                                                                                                                                                  aatGCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGcc                       >  1:107295/1‑251 (MQ=255)
                                                                                                                                                                                                                  aaggCCGAGCGCCACCAGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGcc                       <  2:117054/248‑1 (MQ=255)
                                                                                                                                                                                                                                 caGCACCGCCTGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTa         >  1:151660/1‑250 (MQ=255)
                                                                                                                                                                                                                                      cccGCCTGAATAAGGCGAATCTGGAGCGACGGGCCACGGCCAAAAATTGGCTTCATAAGTTATCCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGCCCCTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACaaa    <  2:437149/250‑1 (MQ=255)
                                                                                                                                                                                                                                        cGCCTGAATAAGGCGAATTTGTAGGGACGGGCCACGGCTAAAAATTGGCTGCATAAGTTATGCGTATTCTCGGATAAGACCAGGCAGGGTAAACAGACCCTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACg  <  2:45671/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       
TCCAAGCATCAGCAGTTCACTGTTTTTCAGCAACAGTTCCTGACGTAACGCCCGGTTTTCAAGTTCTAATTGGTCACGCGAGGCCAGCGTCTGCGATACGCCATCCAGCAATTCACGAGGAGCATTGGAAACAAAGTAGAAAGGACTGACGGCGGTATCCATATAAGTACGGATTTGACTGAACGTCCCCAGGCGGCTGTCGGCAATAATAATGCCGAGCGCCACCAGCACCGCCAGAATAAGGCGAATCTGTAGCGACGGGCCACGGCTAAAAATTGGCTTCATAAGTTATGCGTATTCTCGTATCAGACCAGGCAGGGTAAACAGACACTTCCCCTGCCTGCATCCGATTACTCTTCGCTGAACAGGTCGCCGCCGTGCATGTCGATCATTTCCAGCGCTTTGCCGCCACCGCGCGCCACACAGGTCAGCGGGTCTTCAGCAACAACGACTGGAATGCCGGTTTCTTCCATTAACAAACG  >  NC_000913/3399694‑3400175

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