Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,336,901	0	T	A	26.1%	79.7 / 14.1	46	Y44F (TAT→TTT)	adiC	arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (28/6);  new base A (0/12);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 4.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACC‑ACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAAT  >  NC_000913/4336651‑4337143                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    gATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACC‑ACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTGTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGGAAGCGTAAGAAACACCAGGACTTGGGCCGAGGGACG‑ACATTTTAGCGTATACCATCGAGAGCCCCCCTGCACCGATAATCGTCACAAAccacccat                                                                                                                                                                                                                                                     >  2:283702/1‑251 (MQ=255) gATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACC‑ACCATGGCGATATTGCCGATCCAGCAGGCCAGACAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAAccatccat                                                                                                                                                                                                                                                     <  1:110669/251‑1 (MQ=255)     gTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACC‑ACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑t                                                                                                                                                                                                                                                <  1:35724/251‑1 (MQ=255)                  gtacttccagagccggaaagaagTACCTAAAATA‑CCCACCCCCATGACCCACCAGGGCGAGATGGCCGTCCCAGCAGCCCAGCCAGTGGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCTTAAGAACCACCAGGACTTGGGTCGAGGAACG‑TCATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGAAAAGCTTCACCTTCCATCCATAAA‑TGGCAATTCCGCCa                                                                                                                                                                                                                                   <  2:140989/234‑1 (MQ=255)                       ttCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACC‑ACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGGATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCCACCATCCATAAA‑AGGCCATCCCGCCAgc                                                                                                                                                                                                                                 >  2:186468/1‑247 (MQ=255)                           gAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACC‑ACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGcc                                                                                                                                                                                                                          >  1:83085/1‑251 (MQ=255)                              tCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACC‑ACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGGAAGCGTAAGAACCACCAGGAATTGGGTCGAGGAACG‑GCATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCTTCACCCACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAgg                                                                                                                                                                                                                       >  2:216351/1‑251 (MQ=255)                                  gaagaaGTAACTTAAATATCCTACGCCAATGACC‑ACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCCCCAGGACTTGGGTCGAGGAACG‑ACCTTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGCTAATCATTACCAACCAACCATAAA‑TTGCAATCCCACCAggaaaggcaaagtttt                                                                                                                                                                                                                   >  2:239348/1‑236 (MQ=255)                                                 atatCCTACGCCAATGACC‑ACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCCGGCCTTGGGTCGAGGCACG‑ACATTTTGGCGTATACCATCGAGAGCCCCCGTTCACCGGTAATCGTCCCCCACCCTCCATAAA‑CTGCAATCCCGCCAGCAGAGGCCAGGCTTGCCGGTAACAGaaaa                                                                                                                                                                                                     >  2:168086/1‑250 (MQ=255)                                                     ccTACGCCAATGACC‑ACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCAACGAGCGCCCCCGTGCCCCGATAATCGGCACCAACCAACCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGCAAAacaac                                                                                                                                                                                                >  2:441783/1‑249 (MQ=255)                                                       tACGCCAATGACC‑ACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                                                                                                                                                                                              >  1:126073/1‑251 (MQ=255)                                                                        agggCGATAGGGCCGACCCAGCAGGCCAGCCAGTAGGGGCCGTTGTTTTAAAAACC‑GAGAAACGGCCCAAAG‑AGGCGGCGGGCGTAGGCGTAAGAACCCCCAGGACTTGGGTCGAGGAACT‑ACATTTTGGGGTATCCCACGGAGGGCCCAAGTGCCCCGATAACCGACCCCAACCATCCAAAAA‑GGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTccc                                                                                                                                                                              <  1:255554/249‑1 (MQ=255)                                                                        aTGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGGACG‑ACCTTTTGGCGGATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCCTAAA‑TGGCAAACCCGCCAGAAAAGGCCAGGGTTGCCGGTAACAGAAAAACAACTTACCCCATAAAATTccc                                                                                                                                                                              >  2:334349/1‑251 (MQ=255)                                                                            cGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGac                                                                                                                                                                          >  1:368671/1‑251 (MQ=255)                                                                                      gATCCAGCAGGCCAGCCAGTAGAGGACGTTTTTTTGATACCC‑AAGAAACGGGCAAAAG‑CAGCGGGGGGCGTAAGCGAAAGACCCCCCAGGACTTGGGTCGAGGAACG‑CCATTTTGGCGTATCCCACCGAGGGCCCCAGTGCCCCGATAATCTCCACCACCCACCCAAAAA‑GGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAggg                                                                                                                                                                <  1:399074/251‑1 (MQ=255)                                                                                              agGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGCGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCCTAATATTCCCCGACACCCTCAGGGGGGCgggg                                                                                                                                                        >  2:376937/1‑251 (MQ=255)                                                                                                 ccagccagTAGAGGACGTTGTTTTGATACCC‑AAAAAACGGGCAAAAG‑CAGCGGGGGGCTTAAGCGAAAGAACCACCAGGACTTGGGTCGAGGACCG‑ACTTTTTGGCTTATCCCATCGAGACCCCCAGGCCCCCAATAATCGCCCCCACCCACCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAtt                                                                                                                                                     <  1:186468/251‑1 (MQ=255)                                                                                                  caagccagtggggaagTTGGTTTGAAAACCCAAGAAACGG‑CCAAAGGCAGCGGCGGTCGTAGGCGTACAACCCACCGGGCTTTGGTTCGGGAAACG‑ACTTTTTGGCGTAAACCAGCGAGAGCCCCAGTACACCGAGAACCTACACAAACCA‑CCAAAAAGTGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAtt                                                                                                                                                     <  2:300856/236‑1 (MQ=255)                                                                                                            agGACGTTGGTTTTATAACC‑GAGAAGCGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGACCCACCAGGACTGGGGTCGAGGAACG‑ACATTTGGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGCCACCAACCATCAATAAA‑TGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACttt                                                                                                                                          <  1:371661/251‑1 (MQ=255)                                                                                                            agGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACttt                                                                                                                                          >  1:69829/1‑251 (MQ=255)                                                                                                             ggACGTTGGTTTGAAAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtg                                                                                                                                         >  2:339295/1‑251 (MQ=255)                                                                                                               aCGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGACCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgtg                                                                                                                                       >  1:153156/1‑251 (MQ=255)                                                                                                                        ggATAACC‑GGGAACCGGGCCAAAG‑AGGGGGGGGGGGTAGGCGAAAGACCCCCCAGGCCTGGGGACGAGAAACT‑CCATTTGGGCGTACCCCATCAAGAGCCCCGGCGCCCCAAAACTCGCCACAAACCAACCAAAAA‑TGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcagc                                                                                                                              <  2:52407/250‑1 (MQ=255)                                                                                                                                     aaCGGGCCAAAG‑CAGGGGGGGGCGTAAGCGTAAGACCCACCAGGACTGGGGGCGAGGAACG‑ACTTTTGGGGGTATACCCCCGAGAGCCCCGGTGCACCAAAAATCGTCACCAACCATCAAAAAA‑GGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACAt                                                                                                                  <  2:132287/251‑1 (MQ=255)                                                                                                                                                    gcggggGGCGAAAGCGTAAGACCCCCAAGGTCTGGGGCCGAGGAACG‑ACTTTTGGGCGTAGCCCACCGAGGGCCCTAGTGCACCGAAAATCGCCCCCAACCACCCAAAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGc                                                                                                    <  1:83459/251‑1 (MQ=255)                                                                                                                                                           gcgtaagcgtaaGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATCCCATCGAGACCCCCAGTGCACCGATAATCGTCACCACCCATCCAAAAA‑TGGCATTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCg                                                                                             <  1:193889/251‑1 (MQ=255)                                                                                                                                                                agcgtaagAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCtt                                                                                        >  1:113034/1‑251 (MQ=255)                                                                                                                                                                 gcgaaagAACCACCGGGACTGGGGTCGAGAAACA‑ACATTTTGGCGTATCCCATCGAGAGCCCCAGTGCCCCAATAATCGCCCCCACCCATCAATAAA‑GGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTc                                                                                       <  1:376937/247‑1 (MQ=255)                                                                                                                                                                         cccaccaGTACTTGGGTCGAG‑AACGGACTTTTGGGCGTATCCCAGCGAGAGCCCCAGTGCACGGATAACCGTACCCAACCATCCATAAA‑GGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAAt                                                                               <  1:368464/250‑1 (MQ=255)                                                                                                                                                                                 gACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACg                                                                       >  1:390635/1‑251 (MQ=255)                                                                                                                                                                                  aCTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGAAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTACGAAGGCGAGCTTCGTGACAATAAATTACGt                                                                      >  2:48312/1‑251 (MQ=255)                                                                                                                                                                                   cTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                                                     >  1:75653/1‑251 (MQ=255)                                                                                                                                                                                      gggTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                                                   >  2:31936/1‑250 (MQ=255)                                                                                                                                                                                            agaaaCG‑ACTTTTTGGCGAATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCACCCACCCAAAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATAc                                                            <  1:334349/248‑1 (MQ=255)                                                                                                                                                                                                  g‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGc                                                      >  1:194228/1‑251 (MQ=255)                                                                                                                                                                                                        tttGGCGTATCCCCTCGAGGGCCCCAGTTCCCCAATAACCCTCCCCACCCCTCCTAAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGTCAATAAATTACGTCATATCATACGCCTGCATTc                                                  <  1:283702/250‑1 (MQ=255)                                                                                                                                                                                                                  cccATCGAGAGCCCCAGTGCACCGAAAACCCACACCACCCATCAAAAAA‑TGCCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAat                                       <  1:21826/250‑1 (MQ=255)                                                                                                                                                                                                                  aCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGacac                                                                                                                                                                        >  2:384542/1‑122 (MQ=255)                                                                                                                                                                                                                     aTCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat                                    >  2:330400/1‑251 (MQ=255)                                                                                                                                                                                                                       cGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattatt                                   >  1:32025/1‑250 (MQ=255)                                                                                                                                                                                                                              cccAGTGAACCGAAAATCCACACCAACCACACAAAAA‑GGGCACTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg                           <  2:69829/251‑1 (MQ=255)                                                                                                                                                                                                                                      ccaaaaagtgctCACCAACCACCAAAAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCg                   <  2:126073/240‑1 (MQ=255)                                                                                                                                                                                                                                           tAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGAt              >  2:142288/1‑251 (MQ=255)                                                                                                                                                                                                                                             aTCGTCACCAACCATCCATAAA‑TGGCAATACCGCCAGTAGAGGCCAGGATTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGc            >  2:296908/1‑251 (MQ=255)                                                                                                                                                                                                                                                 tCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGtt        >  2:299070/1‑251 (MQ=255)                                                                                                                                                                                                                                                   accaaccaTCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGttat      >  2:357836/1‑251 (MQ=255)                                                                                                                                                                                                                                                    ccaaccaTCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGttatt     >  1:184932/1‑251 (MQ=255)                                                                                                                                                                                                                                                       accaTCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAat  >  1:434979/1‑251 (MQ=255)                                                                                                                                                                                                                                                        ccatccatAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAat  >  2:253559/1‑250 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    GATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACC‑ACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCCAAAG‑CAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACG‑ACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAA‑TGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAAT  >  NC_000913/4336651‑4337143

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