Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	734,952	T→G	14.3%	pseudogene (820/1521 nt)	rhsO →	pseudogene, Rhs family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	734,952	0	T	G	14.3%	181.9 / 4.0	89	pseudogene (820/1521 nt)	rhsO	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base T (28/45);  new base G (9/7);  total (37/52)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGGTTCAGTATATAGATTCAATGGG  >  NC_000913/734715‑735190                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               ggggtggtgTCCCCGCCGGTGCTGGGGCAGATGCTCGACCGGCTGGAAGGTGAAACCCTGGCTGCCCGGGTGAGTGAGGAAACCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGt                                                                                                                                                                                                                                                <  1:259228/235‑1 (MQ=255)   tgttgtGTCCCCCCGGGTTGTGGTGCAGTTGCTCGACAGGCTGGAAGGTGAAATCCCGGCTGACCGGATGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCCTCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg                                                                                                                                                                                                                                             <  2:434282/239‑1 (MQ=11)   tgtgGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGGGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCACCACGGAAGGGGCAACAGAATGGGgc                                                                                                                                                                                                                                              >  1:637605/1‑238 (MQ=255)    gggttgcccccGGCGGTGGTGGGGGAGGTGCCCGACGGGGTGGAAAGTAAAATCCTGGCGGCCGGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg                                                                                                                                                                                                                                             <  2:911409/230‑1 (MQ=11)    ggggtgTTCCCGCGGGGTGTGGGGGAGATGCCCGGCCGGGTGAAAGGTGAAATCCGGGCTGACCGGGGGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg                                                                                                                                                                                                                                             <  1:593539/236‑1 (MQ=11)     tggtgTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGGgcgcg                                                                                                                                                                                                                                           >  1:836612/1‑238 (MQ=37)       gtttccccGCCGGTGCTGGTGCAGATGCTCGCCCGGCTGGAAAGTGAAAGCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAg                                                                                                                                                                                                                                          <  1:1076441/233‑1 (MQ=255)         gtcccggccgGTGCGGGGGGAGATGCTCGACCGGGTGGAAAGGGAAATCCTGGCTGACCGGGGGGGGGGGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGaa                                                                                                                                                                                                                                        <  1:226024/232‑1 (MQ=11)         gtTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCCGATGGACCCGGGGTCCACG‑CCGGCGCGAAAAATCCACCAGGACCACGGCGACCATCGGGGCCTGCCCCCGGGGCTTTTCAACACGGAACGGGCACAAGACTGGGCGGCAGaaa                                                                                                                                                                                                                                       >  2:813905/1‑238 (MQ=11)         gtTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCCCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGGGCGCAGAAt                                                                                                                                                                                                                                       >  1:990779/1‑239 (MQ=37)            cccGCCGGGGTGGGGCCAGATGCCCGACGGGTGTGAAAGGGAAAGCCGGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACg                                                                                                                                                                                                                                    <  1:481366/239‑1 (MQ=11)             ccgccgGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGACTGGTGCGCCGAATACGa                                                                                                                                                                                                                                   >  1:922185/1‑239 (MQ=37)             ccgccgGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGGGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGCCCATCGCGGCCCGCCGCGGCCCCTTGTCAGCACGGAAAGGGCACCAGAATCGGGCGAAGAAAACGa                                                                                                                                                                                                                                   >  2:292720/1‑239 (MQ=11)               gccgGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGCAGGGGCAACAGAATGTGGCGCAGAATACGATg                                                                                                                                                                                                                                 >  1:625992/1‑239 (MQ=255)                  ggTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGAAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAAt                                                                                                                                                                                                                              >  1:365277/1‑239 (MQ=37)                        ggggCAGATCCCCGACCGGCTGGAAAGGGAAATCCTGGCTGACCGGGTGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCa                                                                                                                                                                                                                        <  1:1095692/236‑1 (MQ=255)                           cacaaggccagACCGGGTGGAAAGGGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAAc                                                                                                                                                                                                                      <  1:758404/228‑1 (MQ=255)                              gATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCTTCGTGCGGCCTGACGGGGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGCCCTGCCCCTGGCCCGTGTCACCACCGAACGGGCCACAAAATGGTGCGCAAAATACGa                                                                                                                                                                                                                                   >  2:401927/1‑222 (MQ=11)                              gATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCCTCGCGGCCTGCCCCTGGCGCTTGGCAGCACGGCAGGGGCAACAACACGGGGCGCGCAAAACGATGAATGGGGCAACcagc                                                                                                                                                                                                                  >  2:763357/1‑239 (MQ=11)                                tGCTCGACCGGCTGGAAAGTGAAACCCTGGCTGACCGGGGGAGTGGGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAACCCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACctgc                                                                                                                                                                                                                  <  1:1014149/237‑1 (MQ=255)                                   tggccgggCGGGAAAGTGAAATCCGGGGTGGCCGGGGGGGGGGGGAAAGCCCCCGCGGGCTGGCATCGTGCGGCCTGACGGTGGAGAAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCtgaat                                                                                                                                                                                                             <  2:59289/233‑1 (MQ=11)                                    cGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGGACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCGGCCGCTGGCGCTTGTCAGCCCGGAAGGGGCAACAGAAAGGTGCGCAGAAAACAATAAATGGGGCAACCCGCtgaatg                                                                                                                                                                                                            >  2:649008/1‑239 (MQ=11)                                                     aaaTCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGACGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCAGCCTCa                                                                                                                                                                                            >  1:393660/1‑238 (MQ=255)                                                        tCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCt                                                                                                                                                                                         >  1:278543/1‑238 (MQ=37)                                                        tCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTg                                                                                                                                                                                        >  1:462740/1‑239 (MQ=37)                                                           tggttggcggggtgggtGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcag                                                                                                                                                                                     <  1:649008/230‑1 (MQ=255)                                                           tGGCTGACGGGGGGGGGGGGGAAGCCCGCCGCTGGCTGCAATCGTGGGGCCTGACGGGGGAGAAGAGG‑CAAAACCAGATGGCCCCGGTGTACACG‑CCGGCGCGAAAAATCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcag                                                                                                                                                                                     <  2:462740/239‑1 (MQ=11)                                                           tGGCTGACCGGGTGAGGGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTgcag                                                                                                                                                                                     <  1:763357/239‑1 (MQ=37)                                                                 aCCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAAGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCAGCAGCCCCagaagcgt                                                                                                                                                                               >  2:624457/1‑237 (MQ=255)                                                                       gtagtgggAAAAGCGCCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTTCACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCg                                                                                                                                                                          <  1:907697/236‑1 (MQ=255)                                                                              ggAAAGCCGCCGGTTGCTGGCATGGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCg                                                                                                                                                                   <  1:912701/238‑1 (MQ=255)                                                                               gAACGCCCCCGCTGGTGGGCTCCGGGCGGCTGGACGGGGGGGCA‑ATGCCAAAACGCGGTGGCCCCGGTGTAGCCG‑CCGGCGCGAAAAACCCACCTGTACCACTGCGCCCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCgg                                                                                                                                                                  <  2:747123/238‑1 (MQ=11)                                                                               caatggcgccgtGGGCTGGCATTGTGCGGCCGGCCGGTGTAGAACAGG‑CAGAGCTAGATGGACCCGGTTGCGCCG‑CCGGCGCTAATAATCGACCTGGACACCTGCTCACATGCCGCCCTGCCGCTGGCGCTTGTCAGAACGGAAGGGGCAACAGAAGGGTGCGCAGAATACGATGAATGGGGCAACCTTCTGAATGAAGAGAACCCGCATCAGCTGCAGCATCTTATCCGCCTGCCgg                                                                                                                                                                  <  2:590968/232‑1 (MQ=2)                                                                                   cccccggctggctgGCATCGTGGGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTTTACCCG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGgcag                                                                                                                                                             <  1:787705/233‑1 (MQ=255)                                                                                    ccctcgtggggtggcCTCGTGGGGGCTGGCGGTGGGCAAGTGG‑CAAACCCAGATGGACCCGGTGTACCCG‑CCGGCGGGAAATATCCACCTGTACCACTGCGACCATGGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGgcagc                                                                                                                                                            <  2:278543/228‑1 (MQ=11)                                                                                     cggcggcgggtggcATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTCAACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGgcag                                                                                                                                                             <  1:645456/227‑1 (MQ=255)                                                                                     cggcgcTGGTGGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑ACAATCCAGATGGCACCGGTGTCAACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGgcagc                                                                                                                                                            <  1:780512/235‑1 (MQ=11)                                                                                     cgccgcTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCGGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGGGCGCGGAATACGATGAAGGGGGCAACCTGCTGAATGAATAGAACCCGAATCAGGTGAAACCGCTTAACCCCCGGCCCGGgcagca                                                                                                                                                           >  2:453627/1‑239 (MQ=11)                                                                                     cgccgcGGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGgcagca                                                                                                                                                           <  1:287799/239‑1 (MQ=37)                                                                                       ccgcTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGCCCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGt                                                                                                                                                         <  1:747947/239‑1 (MQ=37)                                                                                             gctggcATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGCATGGTGCGCAGAATACGGTGAATGTGGCAACCTGCTGAATGAAGAGAAACCGCATCAGAGGCACAAGATTATCCCCCGTCCGGGgcaccggaagaaga                                                                                                                                                   >  2:1043375/1‑228 (MQ=11)                                                                                               tggcATCGTGCGGCCGGACGGTGGAGCAGATG‑CAAAACAAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCCCCTGTACCACTGCGACCAGCGCGGCCTGCCTCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCGGGCTGCCGGGGTAGCAGTAtgatga                                                                                                                                                  <  1:314267/238‑1 (MQ=11)                                                                                               tggcATCGGGCCGCCGTACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACCCG‑GCGGCGCGAAAAATCCACCTGTACCACTGCGGCCATCGCGGCCTGCCGCTGGCGCTTGGCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTAtgatga                                                                                                                                                  <  1:612099/238‑1 (MQ=11)                                                                                                       tGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGAATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCg                                                                                                                                          >  1:714156/1‑238 (MQ=37)                                                                                                         gggccggcgggTGGAGCAGATG‑CAAACCCACGTGGCCCCGG‑GTTCACGCCCGGCGCGAAAAACCCCCCTGCCCCACTGGACCCATCGCGGCCTGCCGGTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATTCGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGcc                                                                                                                                       <  2:732065/230‑1 (MQ=11)                                                                                                          ggCCTGACGGGGGGGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGcc                                                                                                                                       <  1:56948/238‑1 (MQ=37)                                                                                                                        gCAGATC‑CAAAACCAGATGGCCCCGGTGAACCCG‑CCGGCGCGAAAAATCCACCTGTACCACGGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGc                                                                                                                         <  1:327232/238‑1 (MQ=255)                                                                                                                            ttG‑CAAAACCAGGGTGACCCGGGTTCCCCG‑CCGGGCCGAAAATCCCCCCTGTCCCAATGCGACCATCGGGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCg                                                                                                                    <  2:183995/238‑1 (MQ=11)                                                                                                                            aTG‑CAAAACCAGGGGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCCCCTGTACCACGGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGc                                                                                                                         <  1:372284/234‑1 (MQ=255)                                                                                                                              g‑CAAAACCAGATGGCCCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCg                                                                                                                                          <  1:997473/215‑1 (MQ=25)                                                                                                                               gaaaaaCCAGATGGACCCGGTGTACCCG‑CCGGCGCGAAAAATCCCCCTGGACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCg                                                                                                                    <  1:931747/235‑1 (MQ=25)                                                                                                                                       aGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCGGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGAATGATGAGGAGTCGGGCCTGTATTAAAACCGCCCCCCCcattactg                                                                                                           >  2:372284/1‑235 (MQ=11)                                                                                                                                          tGGCCCCGGTGTCCACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGGGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGAcc                                                                                                         <  1:452431/237‑1 (MQ=25)                                                                                                                                                  gtgtACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAggg                                                                                               >  1:718378/1‑239 (MQ=35)                                                                                                                                                   tggcaacG‑CCGGCGCGAAAAACCCACCGTTACCACGGCGACCAACGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAggg                                                                                               <  1:1043375/233‑1 (MQ=255)                                                                                                                                                    gtccacG‑CCGGGGCAAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgggg                                                                                              <  1:977685/235‑1 (MQ=25)                                                                                                                                                      acccG‑CCGGCGCGAAAAATCCCCTGGTACCCCTGCGCCCATCGCGGCCTGCCGCTGGCGCTTGTCAGCCCGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAg                                                                                            <  1:813905/235‑1 (MQ=255)                                                                                                                                                         accccGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGta                                                                                         <  1:363263/236‑1 (MQ=37)                                                                                                                                                            ccggccgaaaaaaCCCCCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTatat                                                                                      <  1:624457/231‑1 (MQ=255)                                                                                                                                                             cGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCCGGGGAGGAATATc                                                                                     >  1:14266/1‑239 (MQ=37)                                                                                                                                                              ggCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATc                                                                                     >  1:500911/1‑238 (MQ=37)                                                                                                                                                                cgcgAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCAc                                                                                   >  1:607410/1‑238 (MQ=37)                                                                                                                                                                cgcgAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTAAAACCCCACCCGCCATTATGACCCGCTGCAGGGGAGGTatataaa                                                                                   >  2:217478/1‑235 (MQ=255)                                                                                                                                                                cccgAAAAACCCCCCCGCCCCCCGGCCACCATGGCGCCCGGCCCGTGGCGTTTGCCAGCACGGAAGGGTCAACAGAATGGTGCGCAGAATCCGATGAATGGGGCAACCTGCGGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACt                                                                                  <  2:500911/237‑1 (MQ=11)                                                                                                                                                                                 tACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCg                                                                                                                                                                   <  1:801184/141‑1 (MQ=39)                                                                                                                                                                                 tACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCAACCGCTATTATGACCCGCTGCAGGGGGGGGAAATAACTCAGGAACCGATTGGGcg                                                                 >  2:1000716/1‑238 (MQ=255)                                                                                                                                                                                 tACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGGGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCg                                                                                                                                                                   >  2:801184/1‑141 (MQ=255)                                                                                                                                                                                        gCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACACCCGCCACCGCTATTATGACCCGCTGCCGGGGAGGTATATTACTCAGGATCCGCTTGGGCGGAAgggg                                                          >  2:901115/1‑239 (MQ=255)                                                                                                                                                                                        gCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCGGGATCCGCTTGGGCTGAAgggg                                                          >  1:507295/1‑239 (MQ=37)                                                                                                                                                                                        gCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAgggg                                                          >  1:449020/1‑239 (MQ=37)                                                                                                                                                                                           aCCATCGCGGCCTGCCGCTGGCGCTTGTCAGCAAGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGAAGTATGATGAGGAGTCCGGCCTGTATTTCAACCGCCACCGCTATTATGACCCGTGGCCGGAGAGGTATATAACTCCGAGTCCGATTGGGCTGACGGGGGGa                                                       >  1:513186/1‑239 (MQ=11)                                                                                                                                                                                                         ccGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGCATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGGGGGATATCACTCAGGAt                                                                            >  1:227389/1‑204 (MQ=255)                                                                                                                                                                                                                      gccgcccGGAAGGGAGAAAAGAATGGTCCGAAAAAATCGATGAAGGGGGCAACCTGCTGAATAAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCCCCGGCCGGGGCAGCGGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAAt                            <  2:607410/233‑1 (MQ=11)                                                                                                                                                                                                                      gccgCAGGGAAGGGGCAAAAAAATGTGGCGAAGAATACGAGGAAGGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCGGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAAt                            <  2:718378/236‑1 (MQ=11)                                                                                                                                                                                                                        acgcaggAAGGGCCAACAGAATGGTGCGAAGAATGCGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGTTTATCCCCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATc                           <  1:453627/233‑1 (MQ=255)                                                                                                                                                                                                                          cACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGcc                                                                                                                                       <  2:310819/128‑1 (MQ=35)                                                                                                                                                                                                                          cACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGcc                                                                                                                                       >  1:310819/1‑128 (MQ=35)                                                                                                                                                                                                                          cACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAAc                                                                                                                            >  1:737883/1‑139 (MQ=34)                                                                                                                                                                                                                          cACGGAAGGGGCAACAGAATGGTGCGAAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGAAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAAc                                                                                                                            <  2:737883/139‑1 (MQ=21)                                                                                                                                                                                                                            cggggggggCAAAAGAAGGGGGCCCACAATACGAGGAGTGGGGCACCCTGCGAAAGGAAGAGACCCCGCATCAGCGGCAGGCGCTTACCCGCCGGCCGGGGAAGCAGTATGATGAGGGTGCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGAGCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGGtt                      <  2:637605/234‑1 (MQ=11)                                                                                                                                                                                                                            cGGAAGGGGCAACAGAATGGGGCGCAAAATTCGGTGAAAGGGGCAACCTGCTGAAGGAAAGAAACCCGCACCGGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTAGGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGGtt                      <  2:507295/239‑1 (MQ=11)                                                                                                                                                                                                                            cGGAAGGGGCAAAAGAATGGGGCGTAGAATACGAGGAGTGGGACAACCGGCGGAATGAAGAGAACCCGCATCAGCTGCATCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGGtt                      <  2:14266/239‑1 (MQ=11)                                                                                                                                                                                                                            cGGAAGGGGAAAAAGAAGGGGGCGCAGATGGGGTTAAAGGGGGCCACATGCGGAAGGGAGAGAACCCGCATCGGGTGCAGTAGCTTGTCAGCCTGCCGGGGCAGTAGTATGATGAGGAGGCCGGGCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGGtt                      <  2:513186/239‑1 (MQ=11)                                                                                                                                                                                                                             ggAAGGGGCAACAGAAGGGTGCGCAGAATTCGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCAGCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGGtt                      <  1:901115/238‑1 (MQ=25)                                                                                                                                                                                                                               aaGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATGACTCCGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGTTTccg                   >  1:307125/1‑237 (MQ=21)                                                                                                                                                                                                                               aaGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGGTTCAg                   >  1:971307/1‑239 (MQ=37)                                                                                                                                                                                                                                   ggCAACAGAATGGGGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAA‑GGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGGTTCAGtata               <  1:451404/238‑1 (MQ=255)                                                                                                                                                                                                                                       aCAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATTTCACTCAGGAACCGATTGGGCTG‑AGGGGGGATGGAATTTTTTTAAGTAACCGCGGAAACCGGTTCAGtataaaaat          >  2:451404/1‑234 (MQ=11)                                                                                                                                                                                                                                       aCAGAATGGGGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGGTTCAGTATATAGa           <  1:292720/239‑1 (MQ=39)                                                                                                                                                                                                                                             tGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCCATTGGGCTGAAGGGGGGATGGACTTTTTATCATTATCCGCTGAATCCGCTTCAGAAAAAAGAATCaag     >  2:499335/1‑238 (MQ=11)                                                                                                                                                                                                                                                tGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGGTTCAGTATATAGATTCAATggg  >  1:179341/1‑239 (MQ=255)                                                                                                                                                                                                                                                tGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACACGCATCAGCTGCAGCAGCTTATCCGCCTGACGGGGCAGCAGTATGATGAGGAGCACGGCCTGTATTACAACCGCCACCGCCAATATGCCCCGCGGCAGGGGAGGTATATCACCACGGGTCCGATTGGGCTGCAGGGGGGATGGAAATTTTATCAGGTTCCGCTGAATCGGGTTCAGTATTAAGATTAAATagg  >  2:308266/1‑239 (MQ=11)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               AGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATG‑CAAAACCAGATGGACCCGGTGTACACG‑CCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGCCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGCTGAATCCGGTTCAGTATATAGATTCAATGGG  >  NC_000913/734715‑735190

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