Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	664,146	G→A	Q349* (CAA→TAA)	rlpA ←	rare lipoprotein RlpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	664,146	0	G	A	98.2%	193.8 / NA	57	Q349* (CAA→TAA)	rlpA	rare lipoprotein RlpA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/32);  minor base C (1/0);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGT  >  NC_000913/663852‑664441                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                        tCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACNACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTNACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCGACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAAc                                                                                                                                                                                                                                                                                                    >  2:549993/1‑300 (MQ=255)                   ggTCATCGGCATGTGCAGCAGAGATAAAGTCTGTGCAAGGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGGGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAGATGCACTACCGGCAGGCATCTGACTTTCCGCTGGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTgg                                                                                                                                                                                                                                                                                  <  2:119380/300‑1 (MQ=255)                     tCATCGGCATGTGCAGCATAGATAAAGGCTGTGCAAAGAGCCGTGGGGAGCGCCAGGCGCTTCATGAGACGAGCGGAAAAAATGGCATTCATGGTCTGAACTACGACATCCGGNATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGNGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCt                                                                                                                                                                                                                                                                                <  2:540889/300‑1 (MQ=255)                     tCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGATCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCt                                                                                                                                                                                                                                                                                <  1:226418/300‑1 (MQ=255)                     tCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCt                                                                                                                                                                                                                                                                                <  1:312176/300‑1 (MQ=255)                     tCATCGGCATGTGCAGCAGAGATAAAGGCGGTGCAAAGAGCCGTGGGGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCt                                                                                                                                                                                                                                                                                <  1:213161/300‑1 (MQ=255)                     tCATCGGCATGTGCAGCAGAGAGAAAGGCTGTGCAAAGAGCCGTGGTGGGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGNACTACCGGCAGGCATCTGACTTTCCGNGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCt                                                                                                                                                                                                                                                                                <  2:450953/300‑1 (MQ=255)                               gtgCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCNACTATAGCAAATGCACTACCGGCAGGNATCTGACTTTCCGCGTGACTTTGTTAACGTCATTCACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATCGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGATTTGc                                                                                                                                                                                                                                                                      >  2:455552/1‑300 (MQ=255)                               gtgCAGCAGAGATAAAGGCTGTGCAAAGAGCAGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGc                                                                                                                                                                                                                                                                      >  1:33144/1‑300 (MQ=255)                                      agagATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACAAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGaa                                                                                                                                                                                                                                                               <  2:295113/300‑1 (MQ=255)                                       gagaTAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGNAAATGCACTACCGGCAGGCATCTGACNTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTAGGCTTTGCTGGCGAAt                                                                                                                                                                                                                                                              >  2:462361/1‑300 (MQ=255)                                         gaTAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATgg                                                                                                                                                                                                                                                            >  1:157819/1‑300 (MQ=255)                                          aTAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATggg                                                                                                                                                                                                                                                           <  1:239213/300‑1 (MQ=255)                                                gCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGc                                                                                                                                                                                                                                                     >  2:207389/1‑300 (MQ=255)                                                 cTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAGATGGTATTCATGGTCTGAACTACGACGTCTTGGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGNCATCTGACTTTCCGCGTGACTTTGTTAATGTCATTTACAGAAAGTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCt                                                                                                                                                                                                                                                    <  2:164674/300‑1 (MQ=255)                                                   gtgCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATAACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTgg                                                                                                                                                                                                                                                  >  1:298550/1‑300 (MQ=255)                                                   gtgCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCNTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTgg                                                                                                                                                                                                                                                  >  2:609240/1‑300 (MQ=255)                                                     gCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATNCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACANACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGCATGGGCCAAGCTGGAt                                                                                                                                                                                                                                                >  2:634730/1‑300 (MQ=255)                                                            gCCGTGGTGAGCGCAAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACAAATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTAGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAg                                                                                                                                                                                                                                         >  2:128549/1‑300 (MQ=255)                                                              cGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGAc                                                                                                                                                                                                                                       <  1:393508/300‑1 (MQ=255)                                                                      gcgcAAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCTGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTAACTTTCCGCGTGACTTTGTTAACGCCATTTACAGAAATTGACACATCAGATACCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTCCAAACGTTGCTGCAAGGTACTGGTTTCGGCTTTGCCGGCGCATGGGCAAAGCTGAATCCGCCAGACCgcgcaac                                                                                                                                                                                                                               >  2:233094/1‑297 (MQ=255)                                                                            ggCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAg                                                                                                                                                                                                                         <  2:387503/300‑1 (MQ=255)                                                                             gcgcTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGt                                                                                                                                                                                                                        <  1:102422/300‑1 (MQ=255)                                                                                cTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGGAATTGACACATCAGATGCCTGCTTTAAGCTACTGCGCGGTAGTAATAAGTGACTGTAATTGGGCTTCGGTTTACACCCGTTGCTGCACGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCGAGCTGGATCCGACAGACCGAGCCATTTTGAGTTAc                                                                                                                                                                                                                     >  2:91971/1‑300 (MQ=255)                                                                                                                 ggTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCATACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTg                                                                                                                                                                                    >  1:632837/1‑300 (MQ=255)                                                                                                                 ggTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTAACGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGACGAACTTCTGTCCCAGTTg                                                                                                                                                                                    >  2:209958/1‑300 (MQ=255)                                                                                                                      gAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCGTTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTCTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGCACTTCTGTCCCAGTTGCTGtt                                                                                                                                                                               >  2:46017/1‑300 (MQ=255)                                                                                                                         cTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGt                                                                                                                                                                            <  1:444617/300‑1 (MQ=255)                                                                                                                              aCATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATTTGACTTTCCTCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAAGAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGt                                                                                                                                                                       <  2:137929/300‑1 (MQ=255)                                                                                                                                aTCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGNATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCTGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTg                                                                                                                                                                     <  2:511943/300‑1 (MQ=255)                                                                                                                                aTCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTg                                                                                                                                                                     <  1:593198/300‑1 (MQ=255)                                                                                                                                   cGTGATNGAATTAAAAAAGTGCCCTACGATAGCAAATGCACTACCGGCAGGCATNTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTgcgc                                                                                                                                                                  <  2:561944/300‑1 (MQ=255)                                                                                                                                     tGATGGAATTAAAAAAGTGCCCTACGATATCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGGGACTTTGTGAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAACGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCAc                                                                                                                                                                <  2:271240/300‑1 (MQ=255)                                                                                                                                     ggATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCTTGACTTTGTTAACGTCATTTACAGAAATTGACACGTCAGATGCCTGCTTTACGCTACTGCGCGTTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCAc                                                                                                                                                                <  2:237076/299‑1 (MQ=255)                                                                                                                                        tGGAATGAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGTTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAg                                                                                                                                                             <  1:622228/300‑1 (MQ=255)                                                                                                                                          gAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGcc                                                                                                                                                           <  2:347389/300‑1 (MQ=255)                                                                                                                                                 aaaaGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGc                                                                                                                                                    <  1:1486/300‑1 (MQ=255)                                                                                                                                                 aaaaGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTGAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGTTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGc                                                                                                                                                    <  1:188482/300‑1 (MQ=255)                                                                                                                                                            actaTAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGGCTTTGTTAACGTCATTTACAGAAATTAAAACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAACCGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGCTGCGCACGAGCCTGATCGCTTACGGccacg                                                                                                                                         >  1:579341/1‑297 (MQ=255)                                                                                                                                                               ataGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACt                                                                                                                                      >  1:237934/1‑300 (MQ=255)                                                                                                                                                                   cAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCa                                                                                                                                  >  1:400457/1‑300 (MQ=255)                                                                                                                                                                    aaaTGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATTCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCAc                                                                                                                                 <  1:590897/300‑1 (MQ=255)                                                                                                                                                                    aaaTGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCAc                                                                                                                                 <  1:662978/300‑1 (MQ=255)                                                                                                                                                                    aaaTGCACGACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCAc                                                                                                                                 <  1:294459/300‑1 (MQ=255)                                                                                                                                                                               cGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGc                                                                                                                      <  1:553348/300‑1 (MQ=255)                                                                                                                                                                                     gcaTCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTgg                                                                                                                >  2:391510/1‑300 (MQ=255)                                                                                                                                                                                                  cgcgTGACTTTGTTAATGTCATTTACAGAAATTGACACATCAGATGCCGGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGAc                                                                                                   <  2:436702/300‑1 (MQ=255)                                                                                                                                                                                                                   gTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCa                                                                                  <  1:111331/300‑1 (MQ=255)                                                                                                                                                                                                                     cATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCAtt                                                                                >  2:11152/1‑300 (MQ=255)                                                                                                                                                                                                                                                   gCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGa                                                  >  1:604322/1‑300 (MQ=255)                                                                                                                                                                                                                                                      ttACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCt                                               >  1:527201/1‑300 (MQ=255)                                                                                                                                                                                                                                                          gctactacgtgttagtaATAAATGACGTTAATCGGGCTTCGGTTCACAAACGTCACTACAAGGTACTGGCTTTTGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGTGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGGTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCGGTaa                                           <  2:124142/288‑1 (MQ=255)                                                                                                                                                                                                                                                             aCTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTaacaa                                        <  1:686942/300‑1 (MQ=255)                                                                                                                                                                                                                                                                     gtagtaATAAATGACTGTAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCtgt                                >  1:242418/1‑300 (MQ=255)                                                                                                                                                                                                                                                                        gtgataaATTACTGTAATTGGGCTTCGGGTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGgc                             <  2:281804/297‑1 (MQ=255)                                                                                                                                                                                                                                                                                      tAATTGGGCTTCGGTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTc               >  2:310283/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                   gTTTACAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCAAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGt  >  2:239342/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                        TCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCGGTAGTAATAAATGACTGTAATTGGGCTTCGGTTTGCAAACGTTGCTGCAAGGTACTGGCTTCGGCTTTGCTGGCGAATGGGCCAAGCTGGATCCGCCAGACCGCGCCATTTTGAGTTACGCGACCGGGGACGCCGAACTTCTGTCCCAGTTGCTGTTGGTACTGTTGCGCACGAGCCTGATCGCTTACGGCCCCGACTTGCACCATAAAGTTGCCGCTGGCGCTTTGCGAGACGGCTTGCGGTGTCACCATTGCAGGCGAGGTTGCAGGCGTTGTCGACGGAGCTGTAACAACGGGCTGTGGCGCAGGCGTCGGTTCGCTGCCTTCCAGT  >  NC_000913/663852‑664441

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