Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,587,542	A→G	14.3%	G198G (GGA→GGG)	mrr →	methylated adenine and cytosine restriction protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,587,542	0	A	G	14.3%	193.4 / 3.1	78	G198G (GGA→GGG)	mrr	methylated adenine and cytosine restriction protein
Reads supporting (aligned to +/- strand):  ref base A (13/52);  new base G (8/5);  total (21/57)
Fisher's exact test for biased strand distribution p-value = 4.44e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGGTGTTGGTTGATGGGG  >  NC_000913/4587307‑4587771                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           ccctCGTTTTAGATCCGAAACCCGTCTTCTCCGGTTATGAAGAACTTGCAAAGAGCAGCCCGGACGGTCGGTTAGAACAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGAt                                                                                                                                                                                                                                         <  2:547755/236‑1 (MQ=255)     cGATTTAGAGCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCACGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGAtggt                                                                                                                                                                                                                                      <  1:966230/238‑1 (MQ=255)            gATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              <  1:560915/239‑1 (MQ=255)              tCCCAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGCCGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              <  1:6516/237‑1 (MQ=255)              gccGAAAGCCGACTCGCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              <  1:887284/236‑1 (MQ=255)                  aaaGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGACGTCATTGTTCTGGCTG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGTTGATTTGCATCGTGTTGGCGGTACTGGGGAGGGGGGCACCGgtggtgg                                                                                                                                                                                                                        >  1:1297053/1‑238 (MQ=255)                  aaaGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCCCCTTGTTGATTTGCAGCGGTTTTGCGGTACTGTAGATGGTGGGACCGCTggggg                                                                                                                                                                                                                        >  1:550188/1‑236 (MQ=255)                        gaatctcCCGAACATGAAGAACTTGCAAAGGGCCGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATc                                                                                                                                                                                                                  <  1:781587/236‑1 (MQ=255)                          ctctCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCCCCGGGATGATTTTCAGCGGGTTGGCGGTACTGGGGATGGTGGCATCGCTGGGGTGaattc                                                                                                                                                                                                                  >  1:1116275/1‑233 (MQ=255)                          accccccggcgaTGGAGGACTTGCAAGAGGTAGCCCGGAGCATCTGTGAGATCAGGGTGTAAATGGTCTCCGTGCGGGTGGGGCTGATGATGTTCTGGAAAACTTATGTCAGGTTTCTCCTTTGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCGCCGTGATGATTTGGAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGc                                                                                                                                                                                                                <  2:729183/228‑1 (MQ=255)                              cccGATCATGAAGAACTTGCAAAGGGCCGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTTATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGa                                                                                                                                                                                                            <  2:283800/239‑1 (MQ=255)                                   tCATTAAGAACTTGCAAAGAGCAGCCCGGACGATCGGGTTGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACGTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAAc                                                                                                                                                                                                       <  2:845614/239‑1 (MQ=255)                                       gaagaaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGCCTCGATGGTGTGATATCGCTTGATAAACttt                                                                                                                                                                                                    >  1:454329/1‑237 (MQ=255)                                        aagtaCTTGCAAAGGGCCGTCCGGACGATGGGTTAGATCGGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGc                                                                                                                                                                                                  <  2:732877/239‑1 (MQ=255)                                        aaggaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGGTAAAAGAGCGTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGc                                                                                                                                                                                                  <  1:214822/239‑1 (MQ=255)                                         agaaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTTGCGGTACTGGAGATGGGGGCATCGATGGTGGGGTATGGCTTGATAAAATTGGcc                                                                                                                                                                                                 >  1:156397/1‑239 (MQ=255)                                         agaaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGcc                                                                                                                                                                                                 <  1:659652/239‑1 (MQ=255)                                         agaaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGc                                                                                                                                                                                                  <  1:605578/238‑1 (MQ=255)                                             cTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCCCCGTGATGATTTTCAGCGGGTTGGCGGTACTGGAGATGGTGGCAGCGATGGGGTGATATCGCTTTGTAAACTTGTGCTggg                                                                                                                                                                                             >  1:114992/1‑238 (MQ=255)                                             cTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCAACGGAGGTGTGAGATCGCTTGTTAAACTTGGCCCggg                                                                                                                                                                                             >  1:611709/1‑238 (MQ=255)                                             attGCAAAGGGGAGCCGGGTGGGTTGGGTAGATCAGGAGCTAAAAGGGCTTCGTGATGCGGGGGCTGATGAGGTTCTGGAAAACTTATTTCAGGTTTCCTCGTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTGGCATCGCCTGGGGTAGGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGga                                                                                                                                                                                             <  2:626841/238‑1 (MQ=255)                                               tGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGACTCGATGGGGGGATATCGCTTGGTAAACTTGGCCTGggga                                                                                                                                                                                           >  1:555542/1‑239 (MQ=255)                                                gCAAAGGGCAGCCCGGACGATCGGTTGGATCAGGCGCTAAAAGTGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGaa                                                                                                                                                                                          <  1:950007/239‑1 (MQ=255)                                                          gcCCGGGCTTTCGGTTAGGTCAGGTGGTAAAAGTGCTGCGTGTGGTGTGGGCTGATGTGTTTCTGGAAACCTTATTGCAGGTTTCTCCTTCGTGCATGGAAGGCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGtt                                                                                                                                                                               <  2:446124/239‑1 (MQ=255)                                                           cccGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGGGTTGGCGGTAACGGGGATGGGGGCATCGGTGGGGTGTATTCGCTTGGTAAACTTGGGCTGGGGAAA‑GTTTAttttt                                                                                                                                                                              >  1:748034/1‑238 (MQ=255)                                                            ccGGACGAGCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATTCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGtt                                                                                                                                                                               <  1:669781/237‑1 (MQ=255)                                                                      ggTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGt                                                                                                                                                                   <  1:322365/239‑1 (MQ=255)                                                                      ggTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGt                                                                                                                                                                   <  1:1113843/239‑1 (MQ=255)                                                                       gTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGATGGTGTGATATCg                                                                                                                                                                                                                 <  1:441920/193‑1 (MQ=255)                                                                         ttgATTAGGGGCTAAAAGGGCTTCGTTAGGCGGTGGCTGATGGTGTTCTGGAAAACTTATTGCAGGTTTCTCCTTGGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTgg                                                                                                                                                                <  2:1270284/237‑1 (MQ=255)                                                                          aGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCAACGTGATGATTTGCAGCGTGTTGGAGGGACTGGAGATGGTGGCATCCATGGTGTGATATCGATTGTTTAACTTGGCCTGGGGAAA‑GTTTATTGTCAGGGaaagagatggc                                                                                                                                                               >  1:262356/1‑231 (MQ=255)                                                                           ggtcgggCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCa                                                                                                                                                              <  1:61357/234‑1 (MQ=255)                                                                             tCAGGCGCTAAAAGGGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            <  1:1047488/239‑1 (MQ=255)                                                                             tCAGGAGCTAAAGGAGCTTCGTGATGCCGTGGCTGCGGAGGTTCGGGAAAACTTATTGCAGGTTTCTCCTTCGCCCTGTGAAGTCATTGTTCTGTATG‑GTTTACATCGCCTGGGGTATGGCGGGCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            <  1:251036/239‑1 (MQ=255)                                                                              cAGGGGCTTACAGGGTTTCGTGGTGGGGGGGGGGATGGGGTTCGTGAAACTTGTTTGCTGGTTTTCCCTTCGTGCTTTTAAGCCATTTTTCTGGAGG‑TTTGGCATCGCCTGGGGTATTGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            <  2:921915/238‑1 (MQ=255)                                                                                       aaaGATATTCGTGATGCGGTGGCTGAGGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCGGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTagg                                                                                                                                                 <  1:763572/239‑1 (MQ=255)                                                                                        aaGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTaggc                                                                                                                                                <  1:879715/239‑1 (MQ=255)                                                                                            gCTTCGGTATGGTGTGGTTGATGGGGTTGTGGAATACTTATTGCAGGTTTCTCCGTCGGTCTTTGAACTCATTGTCTTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTaggcaggc                                                                                                                                            <  2:609057/239‑1 (MQ=255)                                                                                             cTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGTGGGCATCGATGGTGGGGTATCGCTTGATAAACTTGGCCTGGGGAAA‑GTTTATGTTCAGGCAAACAGTTGGCAGAGA‑ACTGTAGGCAGGcc                                                                                                                                           >  1:839014/1‑239 (MQ=255)                                                                                              ttCGTGATGGGGTGGCTTATGGTGTTCTGGAAAACTTTTTGCAGGTTTCACCTTCGAGTTTGTGAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGcc                                                                                                                                           <  2:395170/238‑1 (MQ=255)                                                                                                                ttgtggTTTTGGAAAACTTTTTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTAc                                                                                                                       <  2:839014/235‑1 (MQ=255)                                                                                                                      ttCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              >  1:34950/1‑133 (MQ=255)                                                                                                                      ttCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              <  2:34950/133‑1 (MQ=255)                                                                                                                       tttGGAAAACGTATTGCCGGGTTCTCCTTCGCGCTTTGAAGTCATTGTTCGGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGGGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCa                                                                                                                 <  1:1249943/236‑1 (MQ=255)                                                                                                                       tgtgagaagtttATTGGCGGTTTCGCTTTCGGGCTTTGACGGTATTGTTCTGGATG‑GTTTGCATGGGCTGGGTTATGGCGGCAGCCGTGATGAGTTGCAGCGTGTTGACGGTACTGGAGATGGTGGCATCGATGGTGTGATCTCGCTTGATAACCTTGTCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTGGGCAGAAT‑ACTGGAGGCAGGCAAGAA‑TTACAGGCATTTTACGGCGCa                                                                                                                 <  2:790845/228‑1 (MQ=255)                                                                                                                             aaaCTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCActggct                                                                                                           <  1:1049650/238‑1 (MQ=255)                                                                                                                                gttttgtcccttttCTCCTTGGCCCTTTGGATTCATTGTT‑TGGTGGTTTTTGGATCGGCTGGGGAAGTGCGGCGCCCGTGATGATGTGCAGCGGGTTGGCGGTAGTGGGGATGGTGGCATCGATGGTGGGATGTCGCTTGATAAACTTGGCCTGGAGAAA‑GGTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTTCGGAGCTCTGGCTGGGc                                                                                                       <  2:1203954/228‑1 (MQ=255)                                                                                                                                          gTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGATGGGGTGATATCGCTTGATAAACTTGGGCTGGGGAAA‑GGTTATGTTCAGGCAAAACGTTTGCCGAAT‑ACTTTAGGCAGGCCCAAA‑TTTAAGGCCTTTTACGGCGCCCTGGCCGGGCaaaaaaggaa                                                                                             >  1:981761/1‑234 (MQ=255)                                                                                                                                           tttCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGCTTTGCAGGGGGTTGGCGGTACTGGAGAGGGTGGGATCGATGGGGTGATATCGCTTGGTAAACTTGGCCTGGAGAAA‑GGTTATGTTTGGGGAAAACGTTTGGAGAAA‑AATGTAGGGAGGGCCGAATTTAGAGGCTTTTTAGGGCGC‑CTGGGTGGGGAAAAAACGaaa                                                                                            >  2:975280/1‑239 (MQ=255)                                                                                                                                               gccTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGTCGGCCACCGTGATGATGTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTg                                                                                       <  1:842937/238‑1 (MQ=255)                                                                                                                                                  ttCGGGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTgggg                                                                                    <  1:975280/239‑1 (MQ=255)                                                                                                                                                  ttCGCTTTTGGAAGTCTTTGTTTTGGATG‑TTTTGGATCGCTTGGGGTAGGGGGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGGTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTgggg                                                                                    <  2:454329/239‑1 (MQ=255)                                                                                                                                                      cgcTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGTAGCGTGTTGGCGGTACTGGAGAGGGTGGCACCGATGGGGTGATATCGCTTGATAAAATTGGCCTGGAGGAA‑GTTTTTGTTCAGGGAAAACGTTGGCAGAAT‑ACCTTGGGCAGGGCCAAA‑TTACCGGGATTTTATGGCGCACAGGTTTGGGAAAAAAGGAAAC‑AGgggggggt                                                                                >  2:1080649/1‑235 (MQ=255)                                                                                                                                                                 cATTGTTCTGGATG‑GTTTGCATCGACGGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACtt                                                                     <  1:1185873/239‑1 (MQ=255)                                                                                                                                                                     gTTCTGGATT‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTgg                                                                 <  1:1229833/239‑1 (MQ=255)                                                                                                                                                                       tCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAAATTGGCCTGGAGAAA‑GGTTATGTTCAGGCAAAAAGTTGGCAGGAT‑ACTGTATGCCGGCCAGAA‑TTACAGGCATTTTACGGGGCCCTGGCTGGGCAAAAAACGCAAC‑GGGGGGGGTTTTTTCTCCCTTCTgggt                                                               >  1:144665/1‑237 (MQ=255)                                                                                                                                                                         tGGATA‑TTTTGCGTCGCCT‑GGGTATGTCAGCCACCGTGATGATTTGCAGCGTGTGGGCGGTGCTGGAGATGGTGGCATCGATGGAGTGATATCGATTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGAAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTa                                                            <  1:728280/239‑1 (MQ=255)                                                                                                                                                                                 tttGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTTCTGCGGATGGTGGCCGCCATTGGGTGATATCGCTTGATAAACTTGGCCTTGGGAAAGGTTTAGGTTCGGGCAAAAC‑TTGGGAGAATAACTG‑AGGCCGGGCCGGA‑TATCAGGGCTTATTCGGCGGACTGGGGGGGAGCAAAGCAAAACACGGGGGTGTTT‑TTTCCCCTTCTGGTTTTatgttcc                                                      >  2:914026/1‑233 (MQ=255)                                                                                                                                                                                                 tATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGTGCACTGGCTGGGCAAAAAGCGAAAC‑GTg                                                                                       <  1:23430/190‑1 (MQ=255)                                                                                                                                                                                                       gccccccGTGGTGATTGGCGGCGGGGTGGGGGTACTGGGGATGGTGGGAGCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTAGGTTCAGGCAAAACGGTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACGTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAAt                                 <  2:144665/233‑1 (MQ=255)                                                                                                                                                                                                           aCCGTGATGATTTGCAGGGTGTTGGCGGTACTGGAGATGGTGGCCTCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTc                            <  1:1157002/239‑1 (MQ=255)                                                                                                                                                                                                             cgggctgatTTGCAGCGTGTTGGCGTTACTGGAGATGGTGGGTTCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTAGGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑Ga                         <  1:370865/234‑1 (MQ=255)                                                                                                                                                                                                             cGTTATGTTTTGCAGCGTTTTGGCGGTACTGGAGATGGTTGCATCGGTGGTGTGATATCGCTTGATAAACTTGGGCTGGAGAAA‑GTTTATGTTCAGGCAAAACGGTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑Ga                         <  2:981761/239‑1 (MQ=255)                                                                                                                                                                                                             cGTGATTCTTTGCAGCGTGTTGGCGGTATTGGGGATGGTGTCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCTGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑Ga                         <  2:262356/239‑1 (MQ=255)                                                                                                                                                                                                                 aagatTTGCAGCCAGCTGCCGGAACGTGGGATGGTGGCATCGATGGCGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACCGGCATTTTACCGCCCACTGGCTGGGCAAAAAGCGAAAC‑GAGGGGGAGTTCTTACCACTTCCGGAGTTACTTCTCAGGCGCGGGAATGTTCCCAAACCCTC‑Ggggga                     >  2:70861/3‑238 (MQ=255)                                                                                                                                                                                                                  tgatTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGAATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTa                    <  1:871297/239‑1 (MQ=255)                                                                                                                                                                                                                   gatTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTTGGATAGCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTa                    <  1:1284461/238‑1 (MQ=255)                                                                                                                                                                                                                       ttcagtgggtTGGCGTTACTGGAGATGGTGGCATTGATGGTGTGAAATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGgtg               <  1:1143635/231‑1 (MQ=255)                                                                                                                                                                                                                       tggaggatgtTTGCGGTGCTGGAAATGGTGGCATCGATGGTTTGATATCGCTTGATAAACGTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGgtg               <  1:673356/232‑1 (MQ=255)                                                                                                                                                                                                                           gCGTGTTGGCGGTACTGGAGATTGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            >  2:616496/1‑98 (MQ=255)                                                                                                                                                                                                                           gCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            <  1:616496/98‑1 (MQ=255)                                                                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGGGCACTGGCTGGGCAAAAAGGGAAAC‑GTGGGGTGTTTTTTACCACTCCTGGATTTACTTCTCCGGGGCGGGGCTTTGGCCAATCCGCCTGGGGGTAGGGTgtgggtt         >  2:963941/1‑239 (MQ=255)                                                                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTggg                                                                                     <  1:226663/163‑1 (MQ=255)                                                                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTggg                                                                                     >  2:226663/1‑163 (MQ=255)                                                                                                                                                                                                                               tttGGCGGTTCTGGGGATGGGGGCATGGATGGGTGGTTATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGAAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGGTgttggttg        <  2:114992/237‑1 (MQ=255)                                                                                                                                                                                                                                     ggTATTGGAGGTGGGGGCATCGATGGTGTGATACCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGGTGTTGGTTGATgggg  <  2:156397/238‑1 (MQ=255)                                                                                                                                                                                                                                     ggTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGGTGTTGGTTGATgggg  <  1:70861/238‑1 (MQ=255)                                                                                                                                                                                                                                       tACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGGTgttgg           <  1:168593/227‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           CCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGGTGTTGGTTGATGGGG  >  NC_000913/4587307‑4587771

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