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

GGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGC  >  NC_000913/4336654‑4337148                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        ggTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCTTCACCAACC‑ATCCATaaa                                                                                                                                                                                                                                                      <  2:533580/251‑1 (MQ=255)    tAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACTACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATg                                                                                                                                                                                                                                                    <  1:586040/250‑1 (MQ=255)          cAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATcc                                                                                                                                                                                                                                             <  1:549604/251‑1 (MQ=255)          cAATGGATCATTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATcc                                                                                                                                                                                                                                             <  1:259061/251‑1 (MQ=255)                            cGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ACCCATAAATGGCAATCCCGCCAGTAGAGGCCAggg                                                                                                                                                                                                                           >  2:400644/1‑250 (MQ=255)                             gggAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCAAGCAAGAAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGAGTCGAGGAACGACATTTTGGCGTCTACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATAGCAATCCCGCCAGTAGAggacacgtt                                                                                                                                                                                                                          >  1:278145/1‑244 (MQ=255)                                              aaaccctcccacaaagACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGAAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCGGCGGCGGGCGTA‑AGCGTAAGACCCACCAGGACTTGGGTCGAGAAACGACATTTTGGCGTACACCATCGAGAGCCCCAGGGCACCGATAACCGTCCCCAACC‑ATCCATAAATGGAAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGaaaaa                                                                                                                                                                                                         <  1:554368/236‑1 (MQ=255)                                                    tACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGAAACAGAAAAACACCTg                                                                                                                                                                                                   >  2:419411/1‑251 (MQ=255)                                                             gACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAACGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGCCTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCata                                                                                                                                                                                           >  1:239299/1‑250 (MQ=255)                                                             aacccccaTGGCGATATTGCCGACCCAGCAGGCCAGCCAGAAGAGGACGTTGGTTTGA‑AACCCGAGAAACGGGCCAAAGCAGCGGGGGGCGTA‑AGCGTAAGAACCACCAGGACTGGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCCCCGAAAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataa                                                                                                                                                                                          <  1:178072/246‑1 (MQ=255)                                                              accaccATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCg                                                                                                                                                                                                                                           >  1:144776/1‑201 (MQ=255)                                                              accaccATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCg                                                                                                                                                                                                                                           <  2:144776/201‑1 (MQ=255)                                                                     ggggattgtggcctactcgacagGCCACCCGGTAGAGGAGGTTGGTTGGA‑TAACAGGAAAACGGCACAGAGCAGCGGCGGGGGTA‑AGAGTAACACCCCCGATGCCTGGGGTCGAGAAACGACTTTTTGGCGTACCCAATCGAGAGCCCCAGTGCACGGAAAACCGTCACCACCC‑ACCCAAAAATGGAAATCCCGCCAGTAGAGGCCGGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTcccc                                                                                                                                                                                  <  2:490051/231‑1 (MQ=255)                                                                        cGATATTGCCGATCCAGCAGGCCCGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCCATCCCGCCAGTAGAGGCCAGGTTTGCCGGGTACAGAAAAAAAACTGACCCCATAATATTtccccac                                                                                                                                                                               >  2:1026/1‑248 (MQ=255)                                                                        atatactGCCGAGCAAGCGGGCCAGCAAGTAGAGGACGTGGGTTGAA‑AAACCGAGAAACGGGCAAAAGCAGCGGCGGGCGAA‑AGCGCAAGACCCACTAGGACTGGGGTCGAGAAACCACATTTTGGGTTTTCCCATCGAGAGCCCCAGTCCCCCGATAATCGCCACCACCA‑ACCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGac                                                                                                                                                                               <  1:46341/245‑1 (MQ=255)                                                                             gtGCCGATCCAGCAGGCCGGCCAGGAGGGGACGTGGGTGTGA‑AAACAGAAAACCGGGCAAAAGCAGCGGCGGGCGTAGAG‑GAAAGAACCGCCAGGACGTGGGCCGGGGAACGACATTTGGGCTTATCCCACCGAGACCCCCTGTGCCCCAATACTCGGCACCAAAC‑TTCCATAAATGTCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCa                                                                                                                                                                           <  1:485417/249‑1 (MQ=255)                                                                                       agcagGCCAGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCCGGGTGACg                                                                                                                                                                >  1:63279/1‑251 (MQ=255)                                                                                                  cagTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGCGAGCCCCCGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCCGGTAACAGGAAAACACCTGACCCCCTAATATTCCCCGACACCATCAGGGGGACGGGGATTAAGcc                                                                                                                                                     >  2:289211/1‑251 (MQ=255)                                                                                                       gagGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTAGGTCGAGGAACGCCATTTTGGCGTATACCATCGAGAGCACCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTCGAGGCCCGGTTTGCAGGTAACAGAAAAACACCTGCACCCATAATATTCCCCGACACCATCAGGGTGCCGGGGATTAAGCCCACtt                                                                                                                                                >  1:554912/1‑251 (MQ=255)                                                                                                       agaggagggggTTTGAATAACCGAGAACCGGACAAAAGCAGCGGCGGTGGTA‑GGCGAAAGAACCGCCAGGACTGGGGTGGAGGAACGACATTTTGCCGTATACCATCGAGAGCCACAGTGCACCGATAACCGCCCCAAACC‑ATCAAAAAATGGCATTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCAcg                                                                                                                                                 <  2:86701/242‑2 (MQ=255)                                                                                                              ttGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑CGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCCGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagc                                                                                                                                         >  2:569191/1‑251 (MQ=255)                                                                                                                   tttgaaacccGAGAAAGGGGCAAACGCAGCGGGGGGCATA‑AGCGAAAAAACCACCAGGACTGGGGGCGAAGAACAACATTTGGGCGTATCCCAACAAGCCCCCCAGGCCCCCGAAACTCTTCCCCAACC‑ACCCAAAAATGGAAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATGCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatc                                                                                                                                      <  2:109101/242‑1 (MQ=255)                                                                                                                     gA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGATTTGGGTCGAGGAACGACTTTTGGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCCCCAACC‑ACCAAAAAATGCCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcagca                                                                                                                                  <  1:343184/251‑1 (MQ=255)                                                                                                                     gA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTGTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCACGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcagc                                                                                                                                   >  1:338103/1‑250 (MQ=255)                                                                                                                     gA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGGTTGCAGGTAACAGAAAACCACCTTACCCCATAATATTTCCCGACACCATCAGGGTGACGGGGGTTAAGCCCACTTTGTGcgcatcagc                                                                                                                                   >  2:140442/1‑250 (MQ=255)                                                                                                                          aCCGAGAAACGGGCCAAAGCGGGGGCGGGGGAA‑AGCGTAAGACCCACCAGGACTTGGGTCGAGGAACGACATTTTGGGTTATACCATCGAGAGCCCCAGTGCACCGATAAGCGTCACCACCC‑ATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCg                                                                                                                              <  1:108013/251‑1 (MQ=255)                                                                                                                                  ccGGGCCAAGGCAGGGGCGGGCGAA‑CGCGAAAGAGCCACCAGGACGTGGGTAGAGAATCGACATTTGGGGGAATACCACCGAGACCCCCATCGCCCCAATAATCGTCACCAACC‑ATCGATAAAGGGGAAGGCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACGCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                                                                                                                      <  1:549261/250‑1 (MQ=255)                                                                                                                                          aaGCGGCGGCGGGCGTA‑AGCGAAAAACCCCCCCGGACTTGGGCCGAGAAACGATGTTTTGGCGTAGCCCATCGAGAGCCCCAGTGCACCGATAAGCGTCACCAACC‑AGCCAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACt                                                                                                              <  1:480694/251‑1 (MQ=255)                                                                                                                                               ggtgcgggcgTA‑AGCTGAAGAACCACCAGGACTTTGGGCGAGGAACTACATTTTTGCGTATCCCATAGAGAGCCCCAGTGCCCCGATAAGCGTCACAAACC‑ATCCAAAAATGGAAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACCCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGc                                                                                                         <  2:420396/248‑1 (MQ=255)                                                                                                                                                                taagAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                                                         >  2:191217/1‑251 (MQ=255)                                                                                                                                                                    aaCCACCAGGACTTGGGTCGAGGAACGCCATTTAGGCGTATACCATCGAGAGCCCCAGTCCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCAATTTGTGAGCATCAGCATCCGACGACATAATTAAACTCCTGCGAAGGCGAGCTTCGCGACa                                                                                      >  2:43469/1‑250 (MQ=255)                                                                                                                                                                          caGGACTGGGTTCGAGGAACGACATTTTGGCGTATCCCAGCGAGCGCCCCAGTGCACCGATAAGCGTCACCAACC‑ACCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAAtt                                                                               <  1:535105/251‑1 (MQ=255)                                                                                                                                                                               cTTGGGTCGAGGAACGACATTTTGGCGTATCCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ACCCAAAAAGGACAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAg                                                                                                <  1:150981/229‑1 (MQ=255)                                                                                                                                                                               cTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAg                                                                                                >  2:150981/1‑229 (MQ=255)                                                                                                                                                                                       gAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGAAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcat                                                                   >  1:309709/1‑250 (MQ=255)                                                                                                                                                                                         ggAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAAAGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGAGACGGGGCTTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTACTGCGCAGGCGAGATTCGTGACAATAAATTACGTCATATCATACg                                                                >  1:81562/1‑251 (MQ=255)                                                                                                                                                                                          gAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGc                                                               >  1:223799/1‑251 (MQ=255)                                                                                                                                                                                                      ggcgaatccaaacgggagCCCAAGGGCACCGAAAATCGTC‑CCAACCCACCCAAAAATGTCAACCCCGCCAGTAGAGGCCAGGTTGGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATa                                                   <  2:216767/236‑1 (MQ=255)                                                                                                                                                                                                      ggCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAatat                                                                                                                                                                                       >  1:321287/1‑119 (MQ=255)                                                                                                                                                                                                          tatcccATGGAGAGCCCCAGTGCACCAATAATCGTCCCCAACC‑AACCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGtt                                                <  2:550985/246‑1 (MQ=255)                                                                                                                                                                                                                   gagaGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATtattat                                      >  2:482203/1‑251 (MQ=255)                                                                                                                                                                                                                         cccAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg                                >  2:505332/1‑251 (MQ=255)                                                                                                                                                                                                                               gccccGATAATCGCCCCCACCC‑ACCCATAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTa                           <  1:160255/247‑1 (MQ=255)                                                                                                                                                                                                                                 aCCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTCCGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCg                        >  1:306612/1‑251 (MQ=255)                                                                                                                                                                                                                                     aTAAGCGTCACCAACC‑ACCCATAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGa                    <  1:569191/251‑1 (MQ=255)                                                                                                                                                                                                                                            tCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGtt             >  1:445909/1‑251 (MQ=255)                                                                                                                                                                                                                                                        ccatAAATGGGAATCCCGCCAGTAGAGGCCAGGAATGCAGGTAACAGAAAAACAACTGCCCACATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGACTGCATTCGATAGTTGAATTATTATCGATCTGGTTATCGCTGATGCCGTTATTAATATAGc  >  2:563164/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        GGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGA‑TAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTA‑AGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACC‑ATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGC  >  NC_000913/4336654‑4337148

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