Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,010,878	A→G	T379A (ACC→GCC)	unknown →	Type III restriction‑modification system restriction subunit (EC 3.1.21.5)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,010,878	0	A	G	100.0%	140.3 / NA	49	T379A (ACC→GCC)	unknown	Type III restriction‑modification system restriction subunit (EC 3.1.21.5)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/25);  total (24/25)

AGAAATTGAACATCAAAACATCGCCAGCGAAGAAATTGAAGGTGTAACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTACCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAACCGTTTAAAAAACGCATCAAAAATCAAC  >  1/1010631‑1011089                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                    aaaggcaaaatcactgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGccc                                                                                                                                                                                                                  <  1:733705‑M2/205‑1 (MQ=255)             actgatgcgaaaaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAag                                                                                                                                                                                                      <  1:162789‑M2/217‑1 (MQ=255)                        aaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGt                                                                                                                                                                                            >  1:75313‑M2/24‑250 (MQ=255)                        aaccaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGt                                                                                                                                                                                            >  1:571232‑M2/24‑250 (MQ=255)                         accaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCg                                                                                                                                                                                          >  1:94995‑M2/23‑251 (MQ=255)                         accaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCg                                                                                                                                                                                          >  1:605010‑M2/23‑251 (MQ=255)                         accaaaggcgatgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCg                                                                                                                                                                                          >  1:230783‑M2/23‑251 (MQ=255)                                   atgttattcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCACAAGAGGTGAAAGTCGTGTTTaaaa                                                                                                                                                                                 >  1:858713‑M2/13‑250 (MQ=255)                                        attcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGt                                                                                                                                                                           >  1:865805‑M2/8‑251 (MQ=255)                                          tcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGttt                                                                                                                                                                         >  1:492655‑M2/6‑251 (MQ=255)                                          tcctaACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGttt                                                                                                                                                                         <  1:485454‑M2/246‑1 (MQ=255)                                               aCCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAggg                                                                                                                                                                    <  1:253462/251‑1 (MQ=255)                                               aCCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAggg                                                                                                                                                                    <  1:554768/251‑1 (MQ=255)                                               aCCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAggg                                                                                                                                                                    <  1:482110/251‑1 (MQ=255)                                                       tACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGGGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCTAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAggg                                                                                                                                                            >  1:881450/1‑251 (MQ=255)                                                                          tCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGaa                                                                                                                                         <  1:748938/251‑1 (MQ=255)                                                                          tCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGaa                                                                                                                                         <  1:391/251‑1 (MQ=255)                                                                             aTCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTg                                                                                                                                                                        >  1:155846/1‑217 (MQ=255)                                                                                 aaaTCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATa                                                                                                                                  >  1:695091/1‑251 (MQ=255)                                                                                                   gCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGt                                                                                                                  <  1:738679/249‑1 (MQ=255)                                                                                                       gACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGAc                                                                                                                                                                >  1:888838/1‑199 (MQ=255)                                                                                                        aCTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAggg                                                                                                                                                                    <  1:160252/194‑1 (MQ=255)                                                                                                           tATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAgg                                                                                                        <  1:116365/251‑1 (MQ=255)                                                                                                               ccccTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGAtt                                                                                                    >  1:557460/1‑251 (MQ=255)                                                                                                                   tGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTAcc                                                                                                <  1:776771/251‑1 (MQ=255)                                                                                                                    gATTTTGCCTATATCGGCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGAc                                                                                                                                                                <  1:439843/186‑1 (MQ=255)                                                                                                                                atCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCa                                                                                   <  1:7127/251‑1 (MQ=255)                                                                                                                                       aaaccaaaacAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAAc                                                                             >  1:278563/1‑250 (MQ=255)                                                                                                                                        aaccaaaacAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGt                                                                                                                                                           <  1:128099/171‑1 (MQ=255)                                                                                                                                        aaccaaaacAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAAccg                                                                           >  1:599202/1‑251 (MQ=255)                                                                                                                                             aaacAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGAc                                                                                                                                                                >  1:170098/1‑161 (MQ=255)                                                                                                                                                       gaTTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGcaa                                                             >  1:474134/1‑250 (MQ=255)                                                                                                                                                        aTTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTAttt                                                                                         <  1:317629/221‑1 (MQ=255)                                                                                                                                                                      gATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACaaa                                             >  1:764292/1‑251 (MQ=255)                                                                                                                                                                                   ggCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACaaa                                                         <  1:826957/226‑1 (MQ=255)                                                                                                                                                                                       tGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAAc                             <  1:561528/250‑1 (MQ=255)                                                                                                                                                                                       tGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAAc                             <  1:306041/250‑1 (MQ=255)                                                                                                                                                                                             gTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAAtg                                                       <  1:134857/218‑1 (MQ=255)                                                                                                                                                                                             gTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAACCGTTTaa                      >  1:666200/1‑251 (MQ=255)                                                                                                                                                                                             gTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAACCGTTTaa                      >  1:550755/1‑251 (MQ=255)                                                                                                                                                                                             gTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAACCGTTTaa                      >  1:133911/1‑251 (MQ=255)                                                                                                                                                                                                 ggATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAACCGTTTaaaaaa                  >  1:740409/1‑251 (MQ=255)                                                                                                                                                                                                  gATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAACCGTTTaaaaaa                  >  1:614610/1‑250 (MQ=255)                                                                                                                                                                                                      aCTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAACCGTTTAAAAAACGCAt             >  1:180023/1‑251 (MQ=255)                                                                                                                                                                                                          aCGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAACCGTTTAAAAAACGCATCaaa         <  1:246662/251‑1 (MQ=255)                                                                                                                                                                                                             aaaaGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAACCGTTTAAAAAACGCATCaaaaa       <  1:375951/250‑1 (MQ=255)                                                                                                                                                                                                             aaaaGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAACCGTTTAAAAAACGCATCaaaaa       <  1:171740/250‑1 (MQ=255)                                                                                                                                                                                                                  gCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAACCGTTTAAAAAACGCATCAAAAATCAAc  <  1:25620/250‑1 (MQ=255)                                                                                                                                                                                                                                aTCAAACACGCCGAACGTTTATTTGCCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGc                                                    <  1:41740/186‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                    AGAAATTGAACATCAAAACATCGCCAGCGAAGAAATTGAAGGTGTAACCGTGTTTACTAAAATCCCGAAAAATTCTATCAAAATCCCCGTGGCAGGAGGCGGGACTTATTCCCCTGATTTTGCCTATATCGTCAAAACCAAAACAGGCGAGATTTTGAATTTTGTGATTGAAGCAAAAGGCGTGGATGGTTCGGATAACTTACGAAAAGGCGAAGAACGCAAAATCAAACACGCCGAACGTTTATTTACCCAAATCTCAAAAGAGGTGAAAGTCGTGTTTAAAACTCAGTTTGAGGGGGACAGGGTGGCTGAGTTGATTATGAAAACGATAACCTCTTTATCCGGGTAGTTTGGAGGGATTTACCACTATTTTTTCCATCAAACCGCCGCCGACATGCAACAAATGTGCCGACAAATTCGCGGCGGCAGAACCGTTTAAAAAACGCATCAAAAATCAAC  >  1/1010631‑1011089

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