Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,336,901	0	T	A	31.7%	61.5 / 12.6	42	Y44F (TAT→TTT)	adiC	arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (20/8);  new base A (0/13);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.59e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCA‑GCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATT  >  NC_000913/4336661‑4337140                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                aCCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCA‑GCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGATAATCGTCACCAGCCATCCATAAATGGCAAt                                                                                                                                                                                                                                       <  1:135361/251‑1 (MQ=255)  ccAATGGATCTTCCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGAACCAGCAGGCCA‑GCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGTTAATCGTCACCAACCATCCATAAATGGCAATc                                                                                                                                                                                                                                      <  2:400375/251‑1 (MQ=255)          tCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCA‑GCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAAAGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCCCCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCAACGATAATCGTCCCCAACCATCCATAAATGGCAATCCCGcccgt                                                                                                                                                                                                                              >  1:83177/1‑248 (MQ=255)            tttCAGAATCGGGAAGAAGTAACTTAAATCTCCTCCGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCA‑GCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGt                                                                                                                                                                                                                              >  1:169507/1‑249 (MQ=255)                                 aCTTAAATATCCTACCCCAATGCCCACCATGGGGATATTGCCGCCCCAGCAGGCCA‑GCCGGAAGAGGACGTGGGTTTGAAAACC‑GAAAAAGGGGCA‑AAA‑GGGGCGGGGGGGGTAAGCATAAGAACCCCCAGTATTGGGGTCGAGGAACGATATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCCCCGATAATCGTCCCCACCCATCAAAAAAGGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACa                                                                                                                                                                                                       <  1:310952/251‑1 (MQ=255)                                  cTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCA‑GCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGAATTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAAC‑CCCAGTGCACCGATAATCGTCACCAAACATCCATAAATGGCAATCCCGCCAGTAAAGGCCAGGTTTGCAGGTAACa                                                                                                                                                                                                       >  2:429932/1‑250 (MQ=255)                                              aCGCCAATGACCCCCATGGCAATATGGCCGACCCAGCAGGCCA‑GCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGACCCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCt                                                                                                                                                                                            <  1:494193/249‑1 (MQ=255)                                                    aTGACCACCATGGCGATATTGCCGATCCAGCAGGCCA‑GCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCC‑AAA‑GCAGCGGCGGGCGTACGCGTAAGAACCACCAGGACTTGGGACGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGCTAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGGAGAGGCCAGGTTTGCAGGGACAAGAAAAAAAACTGAccccc                                                                                                                                                                                     >  2:109774/1‑249 (MQ=255)                                                        ccacctgggCGATATTGCCAACCCAGC‑GGCCAAGCCAGTAGAGGAGGTTGGTTTGATAACC‑GAGAAACGGCCC‑AAA‑GCAGCGGCGGGCATAAGAGAAAGAACCCCCGGGTCTTGGGTCGAGGAACGACATTTTGGCTTATACC‑ACCGAGAGC‑CCCAGTCCCCCGATAACCGCCCCCACCCATCCAAAAATGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataata                                                                                                                                                                                <  1:18864/244‑1 (MQ=255)                                                        ccaccaTGGCGAAATTGCCAATCCAGCAGCCCA‑GCCAGTAGAGGGCGTTGGTTTGAAACCC‑GAAAAACGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGACCCCCCAGGCTTTGGGTCGAGGAACGACATTTTGGCTTATACC‑ATCAAGAGC‑CCCAGTGCCCCGATAATCGTCCCCACCCTTCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataata                                                                                                                                                                                <  1:105790/251‑1 (MQ=255)                                                                           gccccAGCGGGCCA‑GCCAGTAGAGGTCGTTGTTTTGATACCC‑GAGAA‑CGGGCCCAAA‑GCGGCGGGGGGCGTAAGCGAAAGAACCACCAGGACTGGGGTCGAGGAACGTTTTTTTGGCTTATCCC‑ACCGAGAGC‑CCCAGTCCACCGAAATTCGCCACCACCCACCCAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAggg                                                                                                                                                             <  1:217517/248‑1 (MQ=255)                                                                                 gcgggCCA‑GCAAGTAGGGCACGTGGGTTTGA‑AACCAGAAAAACGGCCA‑AAA‑GAGGGGGCGGGCGAAAGCGTAAGAACCACAAGTTCTTGGGGCGAGTAACGATATTTGGGCGTATCCC‑ATCGAGAGC‑CCCAGTGCACCGAAAATCGCCTCCAACCACCCAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgg                                                                                                                                                       <  1:49609/248‑1 (MQ=255)                                                                                       ca‑gccagTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAtt                                                                                                                                                  >  1:169201/1‑250 (MQ=255)                                                                                        a‑gccagTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAAAGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGATAATCGTCACCAACCATCCCTAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCACCGACACCATCCGGGTGACGGGGATTa                                                                                                                                                 >  1:118495/1‑250 (MQ=255)                                                                                                  gggACGTGGGTTTGAAAACC‑GAGAAGCGGGCA‑AAA‑GCGGCGGCGGGCGAAAGCGTAAGACCCCCCAGGACTTGGGGCGTGGAACGACATTTTGGCGTATACC‑ATCGAGCGC‑CCCAGTGCCCCTATAACCGCCACCACCCACCCAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATAGTCCCCGACACCAGCAGGGTGACGGGGATTAAGCCCACttt                                                                                                                                       <  2:340934/250‑1 (MQ=255)                                                                                                                 tAACC‑GAGAAACGGGCC‑AAA‑GCAGCGGCGGGCGTACGCGGAAGAAACACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGA‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGGAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGAGAAGGGGATTAAGCCCACTTTGTGagcatcagcatc                                                                                                                        >  2:399858/1‑251 (MQ=255)                                                                                                                  aaCC‑AAGAAACGGGCA‑AAA‑GCAGGGGCGGGCGAAAGCGTAAGAACCACCAGGACTGGGGTCGAGGAACGAAATTTGGGCGAATACC‑ATCGAGACC‑CCAAGTGCACCGAAAAGCGCCACCAACCATCAAAAAATGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATcc                                                                                                                       <  2:202065/251‑1 (MQ=255)                                                                                                                   aCC‑GAGAAACGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCg                                                                                                                      >  1:477496/1‑251 (MQ=255)                                                                                                                          aaaagggCC‑AAAAGCAGGGGCGGGCGAAAGCGTAAGAACCCCCAGTACTGGGGTCGAGAAACGCCTTTTTGGCTTATCCC‑ATCAAGACC‑CCAAGTGACCCGAAAATCGTCACCACCCACCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGAc                                                                                                                 <  1:462092/247‑1 (MQ=255)                                                                                                                           aaCGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGAGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACAt                                                                                                               >  1:515581/1‑251 (MQ=255)                                                                                                                                            cggggggcgTAATCGGAAGATCCACCAGCACGTGGGGCGACAATCGACATTTTGGGGTAACCC‑ATCGAG‑GCTCCCAGTGCCCCGTAAACCGACCCCACCCACCCAGAAATAGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCTCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCg                                                                                                <  1:416171/247‑1 (MQ=255)                                                                                                                                                       agcttaacaccccccaGGACGTGGGCCGATGAACCAGATTTTGGCGTATCCC‑ATCGGGAGC‑CCCAGGGCCCCAATAATCTGCACCACCCATCCATAAAGCGCACCCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAGTAAGCCCACTTTGGGAGCATCAGCAGCCGAAGACATAATTAAACTCCTGCGAAGTCGAGCtt                                                                                     <  2:205040/238‑1 (MQ=255)                                                                                                                                                        gcgtaagAACCCCCCGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataat                                                                                                                                                                                 >  2:64059/1‑158 (MQ=255)                                                                                                                                                        gcgtaagAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataat                                                                                                                                                                                 >  2:144319/1‑158 (MQ=255)                                                                                                                                                        gcggaagaccccccaGGACTTGGGGCGAGGAACGACATTTGGGCGTATACC‑ATCGAGAGC‑CCCAGTGCCCCGATAATCGTCCCCACCCATCAATAAAGGGCACTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTTAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTc                                                                                    <  2:379829/239‑1 (MQ=255)                                                                                                                                                                       ggACTTGGGTCGAGGAACGAAATTTTGGCTTATACC‑ATCGAGAGC‑CCCAGGGCCCCGATATGCGTCACCAACCACCAATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTAc                                                                     <  1:540361/251‑1 (MQ=255)                                                                                                                                                                         cctgggttCGAGAAACGACATTTTGGCG‑ATACCCATAGAGCGC‑CCCAGTGCCCCGAAAACCGCCCCCAACCACCCATAAATTGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCGGCGAAGGCGAGCTTCGTGACAATAAATTACGt                                                                   <  1:462430/244‑1 (MQ=255)                                                                                                                                                                           tggggTCGAGAAACTACTTTTTGGCGTATCCC‑ATCGAGAGC‑CCAAGTGCCCCGAAAATCGTCACCAACCATCCAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGt                                                                   <  1:427642/247‑1 (MQ=255)                                                                                                                                                                             ggttCGAGGACCTACAGTTTGGCGAATCCC‑ATCGAGAGC‑CCCGGTGCACCGAAAACCGTCACCACCCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcata                                                               <  1:413930/248‑1 (MQ=255)                                                                                                                                                                                 gggggAACGTCATTTTGGCGTATACC‑ACCGAGACC‑CCAAGTTCACCGATAACCGTCACCACCCATCCAAAAATGCCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcat                                                           <  2:169507/248‑1 (MQ=255)                                                                                                                                                                                   ggaaaggACATTTTGGCTTATACA‑ATCGAGAGC‑CCCAGTGCCCCAAAAATCGTCACCAACCAACCAAAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATAc                                                         <  2:150885/245‑1 (MQ=255)                                                                                                                                                                                        cGACACTTTGGCGTATACC‑ATCGAGAGA‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACACCTg                                                    >  2:362691/1‑251 (MQ=255)                                                                                                                                                                                           acatttgggcgaaacccaacGAGACC‑CCCATGGCACCGAAAACCGTCACCACCCATCAAAAAATGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTg                                                    <  2:188203/230‑1 (MQ=255)                                                                                                                                                                                                              cGAGAGA‑CCAAGTGCACCGATAATCGTCACCAACCATCCATACCTGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGAGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattatt                                >  2:239870/1‑250 (MQ=255)                                                                                                                                                                                                               gagaGC‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTaa                                                                                                                                                >  1:182009/1‑137 (MQ=255)                                                                                                                                                                                                                agagC‑CCCCGTGCACAGATCATCGTCACCAAGCATCCATAAATGGCAATCCCGCCAGGAGAGGACAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATCATATTCCCCGACACCATCAGGGTGACGGGGCTTAAGCCCACTTTGTGAGCCTCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAACTTACGTCATATCATACGCCTGCATTCGCTAGTTGAATTATTATc                             >  2:534044/1‑251 (MQ=255)                                                                                                                                                                                                                agagC‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATc                             >  1:450546/1‑251 (MQ=255)                                                                                                                                                                                                                         aGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGt                     >  2:126914/1‑251 (MQ=255)                                                                                                                                                                                                                                cGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCt              >  2:441595/1‑251 (MQ=255)                                                                                                                                                                                                                                 gATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTg             >  2:151806/1‑251 (MQ=255)                                                                                                                                                                                                                                        gTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGCGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGt      >  1:190961/1‑251 (MQ=255)                                                                                                                                                                                                                                            ccaaccaCCCATAAAAGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGttatt  <  2:477496/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                ACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCA‑GCCAGTAGAGGACGTTGGTTTGATAACC‑GAGAAACGGGCC‑AAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACC‑ATCGAGAGC‑CCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATT  >  NC_000913/4336661‑4337140

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