Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,901	G→A	D479D (GAC→GAT)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,901	0	G	A	97.5%	99.0 / ‑6.3	40	D479D (GAC→GAT)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/20);  minor base T (1/0);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

TATCCGGCCTACAAAATCG‑CAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGT‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCG  >  NC_000913/66667‑67141                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    tATCCGGCATACAAAATCGGCA‑CGCGTAGGCCTGCTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGAGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTCCTTAGCGACGAAATCCGTACTACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGCTAAC‑AATCAcc                                                                                                                                                                                                                                     >  2:116385/1‑250 (MQ=18)  aTCCGGCCTACAAAATCG‑GAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGt                                                                                                                                                                                                                                   <  1:241561/251‑1 (MQ=21)      ggCCTACAAAATCG‑GAGCGCGTAGGCCTTATAAGACGCGTCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTGGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTTATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCTACTACTTAGCGACGAAATCCGTAATACACTTCGGTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGTGGGTATCGTTATC‑AATCACCTTGAtt                                                                                                                                                                                                                               <  2:282992/251‑1 (MQ=18)             caaagtcggactgcgTAGGCCTGATAAGACGTGCCTGTGTACCATCAGGCTGTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACGCTTCGTTCCAGCGCATCGTGTCTTTAAACGCTGGCATGCGGGTATCGTTATC‑AATCACCGTGATTTCAATg                                                                                                                                                                                                                         <  2:190500/239‑1 (MQ=21)                      aacgcgTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑tc                                                                                                                                                                                                               <  2:176297/249‑1 (MQ=25)                           gTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGc                                                                                                                                                                                                          >  1:281273/1‑251 (MQ=25)                                   gATAAGACGCGCCTGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCGACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTAGC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGaactg                                                                                                                                                                                                    <  2:210989/249‑4 (MQ=25)                                   gATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGaactg                                                                                                                                                                                                    <  2:281273/249‑4 (MQ=25)                                        gACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCata                                                                                                                                                                                             <  1:143835/251‑1 (MQ=25)                                            cgcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGGTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCgg                                                                                                                                                                                         >  1:104158/1‑250 (MQ=25)                                               cAGCTTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCCACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTACA‑TCTCGGCGAACTGGCGCATATCGTTGa                                                                                                                                                                                      <  2:85971/251‑1 (MQ=25)                                                               gTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATC‑AATCACCGGGATTTTAATGTCGTGCA‑TCCCGGCGAACCGGCGCCTATCGTTTAGGTTTAGC‑GCATGGcc                                                                                                                                                                     >  2:175021/1‑250 (MQ=25)                                                                      gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACGGGCGCATATCGTTGAGGGTCAGC‑GCATGGCTGAAGAc                                                                                                                                                               <  2:180317/250‑1 (MQ=37)                                                                       ccGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGGTATC‑AATCACCGGGGTTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGCTGAGGTTCAGC‑GCCTGGCCGAAGCCgg                                                                                                                                                             >  2:284497/1‑251 (MQ=25)                                                                             gcgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGC‑GCATGGCTGAAGACGGTATggg                                                                                                                                                       >  1:200112/1‑250 (MQ=38)                                                                             gcgcTTTGCTTAACCGGCCTACACAATAGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAAAGCTCGCAGGCGGGTATCGTTATC‑CATCACCGTGATTTCAATGTCGGGCA‑TCTCGGCGAAACTGCGCATATCGTTGAGGTTCAGC‑GCACGGCTGAAGACGGTCTGGt                                                                                                                                                       >  2:39374/1‑251 (MQ=255)                                                                              cgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCATCGCGTCGTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATTACCGTGATTGCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTg                                                                                                                                                      <  1:299076/251‑1 (MQ=37)                                                                                              ccTACAAAATCTCATCGTGGAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCATGCGGGTATCGTTATC‑AATCGCTGTGATTTCAAGGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTGCAGC‑GCATGGCTGAATACGGTATGGGGCGGGCAACCTGCGAgg                                                                                                                                      <  2:203387/251‑1 (MQ=255)                                                                                               cTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTTTC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGGTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGGGCGCGCCACCAGCggggg                                                                                                                                     >  2:229452/1‑247 (MQ=255)                                                                                                 aCAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGATc                                                                                                                                   >  1:95748/1‑251 (MQ=255)                                                                                                     aaTCGCAGCGTGTAGGCCTGATAAGACGCGTCCGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTTTCGTTATA‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAAATGGCGCATATCGTTGAGGGTCAGC‑GCATGGCTGAAGAAGGGATGGTGCGCGCCACCAGAGAGGATACACg                                                                                                                               >  2:54290/1‑251 (MQ=255)                                                                                                                        gATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCa                                                                                                             <  2:2163/250‑1 (MQ=255)                                                                                                                                       cGTCGAATCAGGCGTTACACACCTGATGCGACTACCTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGGACTGGCGCATATCGTTGAGGTTCAGC‑GCACGGCTGAAGACGGTATGGTGCGCGCCCCCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCttt                                                                                             >  1:272495/1‑251 (MQ=255)                                                                                                                                                   cGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCCAGGATCCCCGCTTCGGAAGCAGTTGGCAGATCCGGGTGCGCTTTCCACAGCGCAtt                                                                                 >  2:219762/1‑251 (MQ=255)                                                                                                                                                            cctgaagcgacTACTTAACGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGAAACGTTATCCACCCACC‑TGGTTTCAATGTCGTGCAATC‑CTGCGAACTGGCGACTATGGTTGCGGTT‑ACTCGCATGGGCGGAGACCGTGTGGTGCGGGCCCCCAGCGGGGGTCCCCACCTCGGCAGCAGTTGGGAGAGCAGGTTGCCCTTCCCACAGCGCCTTCCCCACCgg                                                                        >  2:317909/11‑251 (MQ=255)                                                                                                                                                             ctgATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGc                                                                       >  1:65910/3‑251 (MQ=255)                                                                                                                                                               gATGCGACTACTGAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAg                                                                     <  2:166759/251‑1 (MQ=255)                                                                                                                                                                aTGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGGCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCAACGGCAGt                                                                    <  2:23067/251‑1 (MQ=255)                                                                                                                                                                       tACTTAGCGACGCAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATC‑AACCCCCGGGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCCTATCGATGAGGGTCCGC‑GCACGGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGAACCACGCTTCGGGAGCCGTTGGCAGATCCGGTTGCGCGTTCCCCAGGGCCTTCACCACCGGCAGGTTccgc                                                              >  2:199984/1‑247 (MQ=255)                                                                                                                                                                                  gaAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGAGCCACGCCTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAggtttggccggggggggggg                                                  >  2:72544/1‑232 (MQ=255)                                                                                                                                                                                   aaatccGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCggg                                                                                                     >  2:1366/5‑198 (MQ=255)                                                                                                                                                                                   aaatccGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGt                                                                                                     <  1:1366/195‑1 (MQ=255)                                                                                                                                                                                                       tCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATA‑AATCACCGGGATTTCAATGGCGTGCA‑TCACGGCGAACTGGCGCCTATCGTTGAGGGTCAGA‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGCTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGGGTGCGGGGGTTTCAACgtggtggtac                               >  1:45018/1‑242 (MQ=255)                                                                                                                                                                                                        ccAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCCACGTGCTTTCAATGTCGTGCA‑TCTCGGCGGACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGGGTGCGGTGTTTTCACCGTGTCGATACAGt                            >  1:204025/1‑251 (MQ=255)                                                                                                                                                                                                                   cgTCTTTAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGa                 <  2:104158/251‑1 (MQ=255)                                                                                                                                                                                                                        ttAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACg               <  1:175021/248‑1 (MQ=255)                                                                                                                                                                                                                         tAAACGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTa            <  1:29114/250‑1 (MQ=255)                                                                                                                                                                                                                          aaaCGCTGGCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTaa           <  1:43541/250‑1 (MQ=255)                                                                                                                                                                                                                                 ggCAGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCg     <  1:128417/249‑1 (MQ=255)                                                                                                                                                                                                                                    aGGCGGGTATCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCg  >  1:316291/1‑249 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    TATCCGGCCTACAAAATCG‑CAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGT‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCG  >  NC_000913/66667‑67141

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