Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,946	T→C	Q464Q (CAA→CAG)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,946	0	T	C	94.3%	88.1 / NA	35	Q464Q (CAA→CAG)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/14);  minor base A (2/0);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGT‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGT  >  NC_000913/66705‑67182                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              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)                                                                                                                                                                                                    tatCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGTAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGatcaatc                                  >  1:121731/3‑251 (MQ=255)                                                                                                                                                                                                             c‑aatcaCCGTGATTTCAATGTCGTGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCCGGGAGTGCGGTGTTTTCAACGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGCGATCAATCAGGCTGGCg                         >  2:10024/1‑251 (MQ=255)                                                                                                                                                                                                                                     tGCA‑TCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGC‑GCATGGCTGACGAAGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGCAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGATTGGGt  >  1:320392/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              GACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATC‑AATCACCGTGATTTCAATGTCGTGCA‑TCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGT‑GCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGT  >  NC_000913/66705‑67182

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: ‑