Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	403,108	G→A	L454L (CTG→CTA)	phoA →	alkaline phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	403,108	0	G	A	97.7%	152.5 / NA	44	L454L (CTG→CTA)	phoA	alkaline phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/25);  minor base C (1/0);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AACGGGCGCTGGAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGA  >  NC_000913/402861‑403352                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     aaCGGGCGCTGGAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTaacc                                                                                                                                                                                                                                                   <  2:65562/251‑1 (MQ=255)                 cGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATctct                                                                                                                                                                                                                                   <  2:8029/251‑1 (MQ=255)                     aaaaaGGAGGGTAACCCGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGCGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTCTTGCGGCGTATGGCCCGCATGCCGCCCCTGTTGTTGGACTAACCGGCCCGACCGATCtcttctc                                                                                                                                                                                                                               >  2:260873/1‑250 (MQ=255)                     aaaaaGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTa                                                                                                                                                                                                                               <  1:292298/251‑1 (MQ=255)                       aaaGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTac                                                                                                                                                                                                                              >  1:207115/1‑250 (MQ=255)                                  acacGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGcc                                                                                                                                                                                                                  <  1:294984/251‑1 (MQ=255)                                      gCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGctc                                                                                                                                                                                                              <  2:144074/251‑1 (MQ=255)                                                  cACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAAt                                                                                                                                                                                                  <  1:232068/251‑1 (MQ=255)                                                     cGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATaaa                                                                                                                                                                                               <  1:108698/251‑1 (MQ=255)                                                      gCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATaaaa                                                                                                                                                                                              <  2:85301/251‑1 (MQ=255)                                                         gATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACTAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACcg                                                                                                                                                                                           >  2:282560/1‑251 (MQ=255)                                                                                tgttgCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCt                                                                                                                                                                     <  2:113641/250‑1 (MQ=255)                                                                                 gttgCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTg                                                                                                                                                                    <  2:185117/250‑1 (MQ=255)                                                                                 gttgCCCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGTAACTCCGAATAGGATTCACAAGAACATACCGGCAGTCAGTTGCTTATTGCGGCGTATGGCCCGCATGCCGCCCATGTTGCTGGACTAACCGACCAGACCGATCTCTTCTACTCCATGAAATCCGCTCTTGGTCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGc                                                                                                                                                                   <  1:338786/251‑1 (MQ=255)                                                                                       ccGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGggtg                                                                                                                                                             >  1:378314/1‑251 (MQ=255)                                                                                         ggATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGggtggt                                                                                                                                                           >  2:332993/1‑251 (MQ=255)                                                                                          gagaCCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACAGACCGGCAGTCAGGTGCGTATTGCGGCGTAGGGCCCGCATGCCGCCAATGGTGTTGGACTAACCGACCAGACCGATCTCTTCTACCCCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGtt                                                                                                                                                          <  1:260873/248‑1 (MQ=255)                                                                                               cAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGttttttt                                                                                                                                                     >  1:10647/1‑251 (MQ=255)                                                                                                      ccGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGGATTTTCGCTGCCGGGTGGTTTTTTTGCTGtt                                                                                                                                               <  1:334866/250‑1 (MQ=255)                                                                                                              cACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCa                                                                                                                                       <  2:281057/250‑1 (MQ=255)                                                                                                                                               gatggtgatgAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCAt                                                                                                     >  1:124252/1‑251 (MQ=255)                                                                                                                                                     gatgaGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATGTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTaa                                                                                               <  2:206250/251‑1 (MQ=255)                                                                                                                                                        gaGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTaaaa                                                                                             <  2:110942/250‑1 (MQ=255)                                                                                                                                                            tACGGGAACTCCGAAGAGGATTCACAAGAAAAGACAGGAAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGGTGTTGGACGAACCGACCAGACCGAGCTCGTCTCCACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGTCGTGAAGGTTCGCTGCCGGGTGGGTTTGTTGCTGTTAGTATCGAGACTGAATGGCAGATCACGGgcttatccacccctgggtaaaacatg                                                                                         <  2:54753/250‑25 (MQ=255)                                                                                                                                                               gggAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGa                                                                                        >  1:109618/1‑248 (MQ=255)                                                                                                                                                                      ccGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCCTGCCGCCAATGTTGTTGGACTAACCGAACAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGTGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATgg                                                                              >  2:261305/1‑251 (MQ=255)                                                                                                                                                                      ccGAAGAGGATTCACAAGAACATACCGGCAGTCAGTGGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTTGGGCTGAAATAAAACCGCGCCCGGCAGTGAAGTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATTGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATg                                                                               <  1:407411/250‑1 (MQ=255)                                                                                                                                                                        gAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGGTTTTTTGATGTTAGCAACAAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTg                                                                            >  2:328994/1‑251 (MQ=255)                                                                                                                                                                          cgagGATTCACAAGAACATACCGGCAGTCAGAGGCGTATTGCGGCGTATTGCCCGCATGCCGCCAATTTTGTTGGACTAACCGACAAGACCGATATCTTCTACACCATGAAAGCCGCTCTGGGGCTGAACTAAAACcccgcccg                                                                                                                                                                                     >  2:381892/2‑144 (MQ=255)                                                                                                                                                                          agagGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGACAATGTTGTTGGACTAACCGCCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGccct                                                                                                                                                                                     <  1:381892/144‑2 (MQ=255)                                                                                                                                                                             ggATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATg                                                                       >  2:400445/1‑251 (MQ=255)                                                                                                                                                                              gATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGa                                                                      <  2:124252/251‑1 (MQ=255)                                                                                                                                                                                   acaAGAACATACCGGCAGTCAGTTGCGTATTGCGGAGTATGCCCCGCATGCCGCCCATGTTGTTGGACTCACCGACCAGACCGATCTCTTCTACACCATGAAAGACGATCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTCCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCCTCCGCCCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAAt                                                                  >  1:27398/1‑250 (MQ=255)                                                                                                                                                                                       gAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCGGGGAATTTTCGCTGCGGGGGGGTTTTTTTGCTGt                                                                                                                                                >  1:259958/1‑168 (MQ=255)                                                                                                                                                                                       aaacatgccggcaagacgtagCGTATTGCGGCATATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCAGGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTGTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGTCGCATACGCTCATGGTTAAAACATGAAGAGGGGTGTGGCTATGAAAATAACAt                                                             <  2:144670/231‑1 (MQ=255)                                                                                                                                                                                                ggCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACGCCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAAAGGGATGGTGCTATGAAAATAACATTACTGGTTa                                                    <  1:166244/251‑1 (MQ=255)                                                                                                                                                                                                    gtcagtTGCGTGTTTTGGCGTAGAGCCCGCATGCCGCTTATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACATCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGTCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTCTGTTAGCAATCAGACTTAATAGCTGATCACGGGCGCATACGCTCATGGTGTAAACATGAAGAGGGATGGTGGTATGATAATAACAGTACTGGTTACCt                                                 <  2:42867/250‑1 (MQ=255)                                                                                                                                                                                                             gTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCgg                                       >  1:113772/1‑251 (MQ=255)                                                                                                                                                                                                                      cgTATGGCCCGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGtttt                              >  2:146585/1‑251 (MQ=255)                                                                                                                                                                                                                              ccGCATGCCGCCAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAacca                      <  1:328994/251‑1 (MQ=255)                                                                                                                                                                                                                                    gccgccAATGTTGTTGGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTc                 <  1:43840/250‑1 (MQ=255)                                                                                                                                                                                                                                               tgttgGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCg       >  1:57189/1‑249 (MQ=255)                                                                                                                                                                                                                                               tgtcggACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCg       >  1:57204/5‑249 (MQ=255)                                                                                                                                                                                                                                                gttgGACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTAAACCACCGTCGGCGCTGCCgag     <  2:207115/250‑1 (MQ=255)                                                                                                                                                                                                                                                   ggACTAACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCgagaga  >  2:115840/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     AACGGGCGCTGGAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGA  >  NC_000913/402861‑403352

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