Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,504,890	A→C	22.2%	V219G (GTG→GGG)	ypdE ←	aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,504,890	0	A	C	22.2%	92.5 / 2.1	58	V219G (GTG→GGG)	ypdE	aminopeptidase
Reads supporting (aligned to +/- strand):  ref base A (21/22);  new base C (6/9);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATCACTCGCTGGTGAGGGAGAACCTGAAAAGTGGTATCAAACGTGA  >  NC_000913/2504657‑2505109                                                                                                                                                                                                                                          |                                                                                                                                                                                                                             tgGGCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGATTTTCTCCCAGCAGGCGGTATCAAGCACAAta                                                                                                                                                                                                                         >  2:547352/1‑238 (MQ=255)   ggCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTGACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCACCAGGCGGGATAAAAAACAAAGGc                                                                                                                                                                                                                       >  1:256606/1‑239 (MQ=255)      aCGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCCACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAAACTGCCGCTGGTTAGCCCCGCCCAAACCAAAGTTTTTCCCCCCTCCGGCGGGGTCAACCCAAAAGGCGaa                                                                                                                                                                                                                    >  1:841346/1‑238 (MQ=255)         ccggtgccggtTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTCACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCCGGCGGTTTCAAGCCCAATGGCGAcct                                                                                                                                                                                                                  >  2:754857/1‑236 (MQ=255)         ccggtgccggtTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCCCCATAATCCAAGTTTTTCCCCCACCAGGCGGGACCAAGACCAATGGCGACga                                                                                                                                                                                                                  >  2:115139/1‑237 (MQ=255)         ccggtgccggtTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGAAGGTTAGCCCCGCCCTAATCAAAGTTTTTTGCCCAGCAGGCGGGATCAAACCAAAAGGCGaacc                                                                                                                                                                                                                  >  2:372122/1‑235 (MQ=255)             tgccggtTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCAAGCAGGCGGTATCAAGCACAATGCGGAcggccggc                                                                                                                                                                                                             >  2:146488/1‑239 (MQ=255)                      aGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGCTGGTTAGCCCCGCCATAATCAAAGGTTTTCGCCCAGCCGGCGGTGTCAAGCACAATGCGAACGTCCGGCGCAACcgc                                                                                                                                                                                                     >  2:686240/1‑238 (MQ=255)                            ccgccccgccGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTGTCAAGCACAATGGCGACGTCCGGCGAAACCGCCCggggg                                                                                                                                                                                              >  2:313902/1‑237 (MQ=255)                            cccacccgccGTCCGGGCCGCCGGTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCGGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTg                                                                                                                                                                                              <  2:155671/235‑1 (MQ=255)                                        cctgtccgccgTTGGTGAACATATCTGCCTGCACCGGAACGCGAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCccg                                                                                                                                                                                  <  2:528656/234‑1 (MQ=255)                                              ccccgTTGCTTAACATTTCCGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGt                                                                                                                                                                            <  2:226272/237‑1 (MQ=255)                                                 cgTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCAGCGGCCCGTTACCAATCTGGCGCTGGTTAGCCGCGCCATAATCAAAGGTTTTCGCCCAGCAGGCGGTATCAAGCACAAAGGCGACCCCCCGCCAACCCGCCCGGGGGGCCGGTTTgccgcccggcaaa                                                                                                                                                                         >  1:855817/1‑232 (MQ=255)                                                     gCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGGGACCCCCCCCGGGAGGCAAGTTGCCCGCCGCGTAATccc                                                                                                                                                                      >  2:547961/1‑238 (MQ=255)                                                     gCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCCCGCGGGTGGCAGTTTGCCCGCCCCGGAATccc                                                                                                                                                                      >  2:504482/1‑238 (MQ=255)                                                              atCTGCCTGCAAGGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGc                                                                                                                                                            <  2:655648/239‑1 (MQ=255)                                                              atCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCCCAATGCCGACGTCCGGCGACACCGCCCGGGGGGCAGTTTGCCCGCCCCGTAATCCCCCctctcccc                                                                                                                                                            >  2:423463/1‑235 (MQ=255)                                                              atCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGGATCAAGCACAATGGCg                                                                                                                                                                                                                      >  2:583296/1‑181 (MQ=255)                                                                ctgcctgcAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGTCGGTAAGTTTTGGCGGCGCAATCCGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCCAAGTTTTTCGCCCAGCCGGCGGGATCAAGCACAACTGGGACGTCCCGCTCCACCCCCCGGGGGGCAGTTTTGCccccgccgaaacacccaccatcccgg                                                                                                                                                          >  1:538012/1‑219 (MQ=255)                                                                ctgcctgcAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGGGGCAGTTTGCCCGCCGCGTAATCCCACCTCTCCGCTg                                                                                                                                                          >  2:178057/1‑239 (MQ=255)                                                                   cctgcAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCCCAATGGCGACGTCCGGCGACACCCCGCGGGTGGCCGTTTTCCCGCCCCCGAATACCACCTCTTCGCTGGa                                                                                                                                                        >  2:656906/1‑238 (MQ=255)                                                                        aaCGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCATTTTGCCCGCCGCGTAATCCCCCCCTCTCGCTGGAAcctt                                                                                                                                                   >  2:795152/1‑235 (MQ=255)                                                                             cccGCCAATTTCTCCCCCCGCGGTTTGCACCCGGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGccacc                                                                                                                                              <  1:809308/236‑1 (MQ=255)                                                                               cccctttttCTCCCCCCTCGTTGTCGACCCCGGTGGTAAGTTTTGGCGTGCCCACCAGGCTCTTGTCGCTTAACCCCAGCCACGCCCCGTTACAAACCGGGGGTTGGTAACCCGGCCCATAATCAAAGGTTTCCGCCCAGCAGGCGGTACCAAGCCCAATGGCGACGTCCGGCGACCCCGCGCGGGGGGCAGTTTGCCCGCCGCGTAAGCCCACCTCTTCGCTGGAACTTGCCACCAGc                                                                                                                                           <  1:13189/233‑1 (MQ=255)                                                                                     tttcgccccctccGGTTTCCACCCGGGGGGTAGTTTTGGGCGGGGCAATACGCGATTTGTCCTTTAACCCCACCACCGGCCCGTTCCCAACCGGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGGGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACt                                                                                                                                     <  1:512272/227‑1 (MQ=255)                                                                                            cgacACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTCCCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTAGCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCACTTCGCTGGAACTTGCCACCAGCCACACTTCCGCa                                                                                                                               <  2:723982/235‑1 (MQ=255)                                                                                              ccACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGGGGCAGTTTGCCCGCCCCGTAATCCCACCTCTTCGCTGGAACTTGCCACCACCCCCACTTCCCCAgga                                                                                                                            >  2:401753/1‑238 (MQ=255)                                                                                                          cccgggcgtgatgTTTTGGGGGCGCAATCCGGGACTTGTCCGTTAACCCCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGGGCCAAAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACCCCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGcgtcg                                                                                                                <  1:432178/227‑1 (MQ=255)                                                                                                                    aaGTTTTGGCGGCGCAATCAGGGACTGGTCGCTTAACCCCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTcg                                                                                                       <  2:354174/238‑1 (MQ=255)                                                                                                                     aGTTTTGGCGGGGCAATCAGCGACTTGTCGGTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTcgcg                                                                                                     <  2:529061/239‑1 (MQ=255)                                                                                                                         ttGGGGGCGAAATCAGTGATTTGTCGTTTAACCCCAGCCTCGGCCCGTTCCCAATGCGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTCCGCCCACCAGGCGGTATCAAGCCCAATGGCGACGTCCGGCGACCCCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGt                                                                                                 <  1:154184/239‑1 (MQ=255)                                                                                                                                    aTCAGCGACTTGTCGCTTAACCCCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACcag                                                                                      <  2:841346/239‑1 (MQ=255)                                                                                                                                               gTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCTGCGTCGTGCAGTTCGCGCAGTAACGTCACCACCAGATAgaca                                                                            >  2:128239/1‑235 (MQ=255)                                                                                                                                                 cgctcaccccctgcctCGGCCCTTTACCAATCGTGCGATGGTTGGCGCCGCCATAATCAAAGTTTTTCGCCCAGCTGGCGGTATCAAGCACAATGGCGACGCCCGGGGACCCCGCGCGGGTGGCAGTTGGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCCCCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCg                                                                         <  1:558548/224‑1 (MQ=255)                                                                                                                                                     tAACCCCAGCCTCGGCCCGTTACCAATCGGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGc                                                                     <  2:299772/239‑1 (MQ=255)                                                                                                                                                      accccccccaccGGCCCGTTTCAAATTTGGCGATGGTTCGCCGCGGCATAATAAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCCCAATGGCGACGCCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGGCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCg                                                                    <  1:138850/228‑1 (MQ=255)                                                                                                                                                          ccAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGACCTTGCACCCAGCCACCCTTCCGCAGGAACGTCGGCGTCGTGCATTTCGCGCCTGAACGTCCACACCAGAAAGCACCCGGAGCggccc                                                                >  1:194474/1‑236 (MQ=255)                                                                                                                                                             gCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGtcatca                                                             <  2:492735/239‑1 (MQ=255)                                                                                                                                                             gCATCGGCCCGTTACCAATCTGGAGGTGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGtcatca                                                             <  2:820225/239‑1 (MQ=255)                                                                                                                                                                   gccgttttccattttggtcgtGGTTACCCGCGCCATAACCAAAGTTTTCCGCCCGGCGGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCt                                                      <  1:372122/219‑1 (MQ=255)                                                                                                                                                                            cactttggcatttgtttccccccccAAAATCAAAGTTTTTCGCCCCGCAGGCGGTTTCAAGCCCAATGGCGACGTCCGGCGACACCGGGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGTGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATc                                             <  1:332929/216‑1 (MQ=255)                                                                                                                                                                            cAATCTGGCGTTGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATc                                             <  2:641082/239‑1 (MQ=255)                                                                                                                                                                            cAATCGGGCGATGGGTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGGAGGCGGTATCAAGCCCAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGGCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATc                                             <  2:680677/239‑1 (MQ=255)                                                                                                                                                                                cTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACCCCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATCACTc                                         <  2:199807/239‑1 (MQ=255)                                                                                                                                                                                 gggggggggtttcccGCCCCAAAATCAAAGTTTGCCGCCCAGCAGGGGGTAGCAAGCCAAGGGGGGACGGCCGGCTACCCCGGGCGGGGGGCAGTTTGCCCGCCGCGTAATCCCACCTCTGCGCTGGAACTTGCCCCCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATCACTCg                                        <  1:754857/226‑1 (MQ=255)                                                                                                                                                                                  ggggttggttcgcccgcccATAATCAAAGTTTTTCCCCCAGGGGGCGGTATCAAGCAACATGGCGCCGTCCGGCGACCCCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATCACTCGc                                       <  1:523629/222‑1 (MQ=255)                                                                                                                                                                                   gCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATCACTc                                         >  2:207076/1‑236 (MQ=255)                                                                                                                                                                                   gCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGGCATCAAAG‑GCTTTCCCCATccccc                                         >  2:332929/1‑232 (MQ=255)                                                                                                                                                                                     gATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTGGTTTCGCGTCCGTCCGTTCGAGCACTAACGCATCCACCCGATCGCAACCGAGCCCGGCACCAAAGAGCTTCCCCCacaaccccgag                                     >  1:2437/1‑229 (MQ=255)                                                                                                                                                                                       tGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGGCCGGCGACACCGCGCGGGTGGCAGTTGGCCCGCCGCGTAATCCCACCTCTGCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATCACTCGCTGGt                                   <  2:256606/238‑1 (MQ=255)                                                                                                                                                                                                   ccATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGCCATCAAAG‑GCTTTCCCCATCCCTCGCTGGGGAGGGAGAAccgg                      >  1:179125/1‑237 (MQ=255)                                                                                                                                                                                                          cAAAGTTTTTCGCCCGGCAGGGGGTATCAAGCCCAAGGGCGCCGTCCGGCCACACCGCGCGGGTGGCAGTTGGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTccc                                                 <  1:92984/205‑1 (MQ=255)                                                                                                                                                                                                             aGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGGCATCAAAG‑GCTTTCCCCATCACTCGCTGGGGAGGGGGAACCTGAAAAggggg             >  2:712724/1‑237 (MQ=255)                                                                                                                                                                                                                ttttCGCCCAGCAGGCGGTATCAAGCCCAATGGCGACGTCCGGCGACCCCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATCACTCGCTGGTGAGGGAGAACCTGAAAAGTGGTAt           <  2:820921/237‑1 (MQ=255)                                                                                                                                                                                                                tttcCGCCCAGCGGGCGGTTTCAAGCACAATGGCGCCTGCCGGCGCCCCCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTCCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATCACTCGCTGGTGAGGGAGAACCTGAAAAGTGGTAt           <  1:178057/237‑1 (MQ=255)                                                                                                                                                                                                                    cGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCACAG‑GCTTTCCCCATCACTCGc                                       >  2:92984/1‑205 (MQ=255)                                                                                                                                                                                                                       ccgccgggCGGTATCAAGCCAAAGGGCCCCGTCCGGCGACCCCGGGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATCACTCGCTGGTGAGGGAGAACCTGAAAAGTGGTATCAAACGTGa  <  1:547352/233‑1 (MQ=255)                                                                                                                                                                                                                               cGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCg                                                                    <  2:584406/166‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                             TGGGCACGCCGGTGCCGGTTAAGTGCACCGCCCCGCCGTCCGTGCCGCCGTTGCTGAACATATCTGCCTGCAACGGCACGCCAATTTCTGCCGCCACGGTTTCGACCCAGGCGGTAAGTTTTGGCGGCGCAATCAGCGACTTGTCGCTTAACACCAGCATCGGCCCGTTACCAATCTGGCGATGGTTAGCCGCGCCATAATCAAAGTTTTTCGCCCAGCAGGCGGTATCAAGCACAATGGCGACGTCCGGCGACACCGCGCGGGTGGCAGTTTGCCCGCCGCGTAATCCCACCTCTTCGCTGGAACTTGCCACCAGCCACACTTCCGCAGGTAGTTCGGCGTCGTGCAGTTCGCGCAGTAACGTCACCAGCAGATAGCAACCGAGGCGGTCATCAAAG‑GCTTTCCCCATCACTCGCTGGTGAGGGAGAACCTGAAAAGTGGTATCAAACGTGA  >  NC_000913/2504657‑2505109

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