Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,809,142	G→A	V201V (GTG→GTA)	waaA →	KDO transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,809,142	0	G	A	98.4%	203.6 / ‑7.3	62	V201V (GTG→GTA)	waaA	KDO transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/28);  minor base C (0/1);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

ATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCT  >  NC_000913/3808895‑3809378                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             aTAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGATCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAAcc                                                                                                                                                                                                                                           <  2:230792/251‑1 (MQ=255)    aaGTCGACCCTAAACGGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGtt                                                                                                                                                                                                                                        <  2:373452/251‑1 (MQ=255)               aaaCTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCt                                                                                                                                                                                                                             >  1:371564/1‑251 (MQ=255)               aaaCTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGGCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCTCGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGcc                                                                                                                                                                                                                              <  2:32600/250‑1 (MQ=255)                aactgggtttGATTGTGGAAACCGAACTAGGGCGTAACCTGATATCGGCGAGACAGAACCGTCACATGCCTCTGGTGACCGCTATCGCGCGAATCTCTGCCCGCTCGGCCGCAGGTTATTCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTTATTGCCGCCCCA‑AATGAAGAAGATGGTGCACGGTTTGTGGCGCTGGGCTCAAAAAAGAACAAGGTCACCTTTTCCGGTAGCCTg                                                                                                                                                                                                                            <  2:269578/243‑1 (MQ=255)                aaCTGGTGTTGATTATGGAAACCGAAATATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTg                                                                                                                                                                                                                            >  1:7568/1‑251 (MQ=255)                          gATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGat                                                                                                                                                                                                                   >  2:305821/1‑250 (MQ=255)                              aTGGAAACCGAACTATGGCCTAACCTGATTGCGGAGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATCTTTc                                                                                                                                                                                                              >  2:169485/1‑251 (MQ=255)                               tGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCt                                                                                                                                                                                                             >  2:130821/1‑251 (MQ=255)                                     ccGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATAGTGCTGTAAcg                                                                                                                                                                                                       <  1:141057/251‑1 (MQ=255)                                       gAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAAcgcc                                                                                                                                                                                                     <  1:311668/251‑1 (MQ=255)                                                 cTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGgct                                                                                                                                                                                           >  1:51035/1‑251 (MQ=255)                                                   aaCCTGATTGCGGCGCTAAATAAACGCAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGACAAACTGGGCAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCACAAAATAACCAGGTAACCGTTCCCGGTAGCCTGAATTTCGCTATTTCTGTAACGCCGCAGTTGgctgc                                                                                                                                                                                         >  2:289548/1‑251 (MQ=255)                                                    aCCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGgctgct                                                                                                                                                                                        >  1:163872/1‑251 (MQ=255)                                                         aTTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGc                                                                                                                                                                                   <  2:279318/251‑1 (MQ=255)                                                              ggcgCTACCTAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGCCTCTCTGCCCGCTCGGCAGCAGGTTATGACAAACTGGGTAAATTCGTCCGTAGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGATACCGGTAGCCTGAAATTCGCTATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTg                                                                                                                                                                               >  1:182848/1‑250 (MQ=255)                                                                    aCATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCt                                                                                                                                                                          <  2:7568/249‑1 (MQ=255)                                                                               aaaTTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTaa                                                                                                                                                                                      >  2:153065/1‑226 (MQ=255)                                                                               aaaTTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTaa                                                                                                                                                                                      <  1:153065/226‑1 (MQ=255)                                                                               aaaTTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTgg                                                                                                                                                             >  1:190323/1‑251 (MQ=255)                                                                                 aTTCCGCTGGAGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑ACTGAAGCAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTggg                                                                                                                                                            >  2:161967/1‑250 (MQ=255)                                                                                               cGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCgg                                                                                                                                             >  1:40555/1‑251 (MQ=255)                                                                                                gCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGt                                                                                                                                            >  2:122541/1‑251 (MQ=255)                                                                                                   aaCGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATg                                                                                                                                         >  2:122616/1‑251 (MQ=255)                                                                                                     cgcgcgACTCTCTGCCCGCTCGGCCGCCGGTTCTGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGAGACGGATTACGCTGATTGCCGCCCAA‑AATGAAGAAGCTGGTGCACGTTTTGTGGCGCTGGGCGCAAAAACTACCCAGGTAACCGTTACCGGTATCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCAGCCAGTGGGCCCCAACCCGCCCGGTATgg                                                                                                                                        >  2:430922/1‑250 (MQ=255)                                                                                                          gACTCTCTGCACGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGC‑CAACAATGAAGAAGCTGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGcc                                                                                                                                  >  2:167191/1‑251 (MQ=255)                                                                                                                          ggCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGGATGGATTGCCACCAGCACTCACGatg                                                                                                                  >  1:152495/1‑249 (MQ=255)                                                                                                                            ccGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGc                                                                                                                >  1:148201/1‑251 (MQ=255)                                                                                                                            ccGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCGGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGc                                                                                                                <  2:371564/251‑1 (MQ=255)                                                                                                                                     aTGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgt                                                                                                       <  1:122541/251‑1 (MQ=255)                                                                                                                                        ccAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCCGCACTCCCGAAGGCGAagggggggtg                                                                                                    >  2:416733/1‑243 (MQ=255)                                                                                                                                               gggTAAAGTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTTGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTTACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATc                                                                                              <  1:394745/250‑1 (MQ=255)                                                                                                                                                    aaTTCGTCCGTCGCTTGCTGCGTCGTATTACGCTTATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCACTGGGCGCAAAAAATAACCAGGTACCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTTGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAATGCGAAGAGAGTGTTGTGATCGCCGca                                                                                        <  2:273994/251‑1 (MQ=255)                                                                                                                                                    aaTTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgccgc                                                                                         >  1:430504/1‑250 (MQ=255)                                                                                                                                                               cGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGACGAAGATGGTGCACGTTTTGTGGCGCTAGGCGCAAGAAATACCCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCCGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGCAGAGAGTGTGGTGATCGCCACACATCAAGCAtt                                                                             >  1:404114/1‑251 (MQ=255)                                                                                                                                                                    gctgcGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTa                                                                         <  1:223410/250‑1 (MQ=255)                                                                                                                                                                       gcGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGCAGGCGAAGAGAGTGTGGTGATCGCCGCCCATCAGGCATTGTTACAGc                                                                     >  2:352488/1‑251 (MQ=255)                                                                                                                                                                       gcGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACa                                                                       <  1:182821/249‑1 (MQ=255)                                                                                                                                                                       gcGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTGTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTCCCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAGGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGc                                                                     <  2:173414/251‑1 (MQ=255)                                                                                                                                                                         gtcgtATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCCCATCAGGCATTGTTACAGCa                                                                    >  2:17185/1‑250 (MQ=255)                                                                                                                                                                          tcgtATTCCGCTGCTTGCCGCCCAA‑AATGACGCAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCATGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAAt                                                                  >  2:420499/1‑251 (MQ=255)                                                                                                                                                                           cgtATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAAtt                                                                 <  2:124556/251‑1 (MQ=255)                                                                                                                                                                           cgtATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGGAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAAtt                                                                 <  2:121099/251‑1 (MQ=255)                                                                                                                                                                              aTTACGCGGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTccc                                                              <  1:17185/251‑1 (MQ=255)                                                                                                                                                                                   gCTGATTTCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCTCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAAtt                                                         <  1:262609/251‑1 (MQ=255)                                                                                                                                                                                                  a‑aaTGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGc                                                                                >  1:320449/1‑213 (MQ=255)                                                                                                                                                                                                  a‑aaTGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTCCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGc                                                                                <  2:320449/213‑1 (MQ=255)                                                                                                                                                                                                     aTGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTa                                        <  2:148201/251‑1 (MQ=255)                                                                                                                                                                                                       gaagaaGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTAcc                                      >  1:122361/1‑251 (MQ=255)                                                                                                                                                                                                        aagaagATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTAccc                                     >  1:70360/1‑251 (MQ=255)                                                                                                                                                                                                          gaagaGGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGGGGTGATCGCCGCACAGCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCg                                   <  1:284517/251‑1 (MQ=255)                                                                                                                                                                                                                  tGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAAtt                                                                 >  1:253627/1‑213 (MQ=255)                                                                                                                                                                                                                  tGCACGTTTTGGGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAAtt                                                                 <  2:253627/213‑1 (MQ=255)                                                                                                                                                                                                                         tttGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCGGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATGGCTCATCCTGGTACCCCGTCATCCAGAACGCt                     <  2:51035/250‑1 (MQ=255)                                                                                                                                                                                                                              ggCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCgg               >  1:348174/1‑251 (MQ=255)                                                                                                                                                                                                                              ggCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCg                <  1:130821/250‑1 (MQ=255)                                                                                                                                                                                                                              ggCGCTGGGCCCAAAAAATAACCAGGTAACCGTTACCTGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGATTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAAGTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCgt               <  2:240387/251‑2 (MQ=255)                                                                                                                                                                                                                              ggCGCTGGGAGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCgg               <  2:207185/251‑1 (MQ=255)                                                                                                                                                                                                                                  cTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGACCGGTATGGATTGCCACCAGCACTCACGAAGACGAAGAGTGTGTGGTGATCGCAGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTACCGGATGc           >  1:214017/1‑251 (MQ=255)                                                                                                                                                                                                                                          aaaaaaTAACCAGGTATCCGTTCCCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCa                                                                    <  1:239767/186‑1 (MQ=255)                                                                                                                                                                                                                                          aaaaaaTAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGGGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTAACGAAGGCGAAGAGAGTGTGGAGATCGCCGCACATCAGGCATTGTTACAGCa                                                                    >  2:239767/1‑186 (MQ=255)                                                                                                                                                                                                                                           aaaaaTAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTaatct  >  2:424813/1‑248 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             ATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCT  >  NC_000913/3808895‑3809378

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