Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,047,415	G→A	T150T (ACG→ACA)	polA →	DNA polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,047,415	0	G	A	91.8%	189.9 / NA	61	T150T (ACG→ACA)	polA	DNA polymerase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/32);  minor base C (4/1);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACGCTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGTCTTGGCGGACTGGATACGCTGTATGCCGAGCCAGAAAAAATTGCTGGGTTGAGCTTCCGTGGCGCGAAAACAATG  >  NC_000913/4047173‑4047661                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        cGTGATGAACTGTTTGAACATTACAAATCCCATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑a                                                                                                                                                                                                                                                 <  2:101438/250‑1 (MQ=255)  gTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AAt                                                                                                                                                                                                                                               <  1:344486/251‑1 (MQ=255)   tgatgaACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACCCTTTTC‑AAt                                                                                                                                                                                                                                               >  2:61406/1‑250 (MQ=255)              ttGAACATTACAAATCCCATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGGCCGGTGCTGCTCAGCACTGGCGATAAAGATATGGCGCAGCTGGGGACGCCAAATATTTCCCTTTTC‑AAAACCATGACgca                                                                                                                                                                                                                                    >  2:219776/1‑248 (MQ=255)                 aaCATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATAcc                                                                                                                                                                                                                                >  1:319334/1‑251 (MQ=255)                    ataaCAAATCACACCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATc                                                                                                                                                                                                                             <  2:111218/248‑1 (MQ=255)                       aCAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGGAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AAAACCATGACGAATACCATCCTc                                                                                                                                                                                                                          >  1:330528/1‑251 (MQ=255)                        cAAATCACAGCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCg                                                                                                                                                                                                                         <  1:307078/251‑1 (MQ=255)                         aaaTCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGGGACGCCAAATATTACACTTATC‑AATAACATGACGAATACCATCCTCgg                                                                                                                                                                                                                        >  1:136817/1‑251 (MQ=255)                           atcccaTCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGAc                                                                                                                                                                                                                      <  2:98907/247‑1 (MQ=255)                            ccgcaTCGCCCGCCTATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGGTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCGCTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTTTC‑AATACCATGACGAATACCATCCTCGGAcc                                                                                                                                                                                                                     <  1:223145/248‑1 (MQ=255)                            ccacaCCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGAcc                                                                                                                                                                                                                     <  1:219776/250‑1 (MQ=255)                                                     aTCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGGCGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCCCTGGCGATAAAGAAATGGCGCAGCCGGTGACGCCAAATATTACACTTATC‑AATAACATGACGAAGACCCTCCCCGGCCCGGAAGAGGTGGTGAATAAGTACGGc                                                                                                                                                                                            >  2:263162/1‑251 (MQ=255)                                                         gCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGcct                                                                                                                                                                                          <  1:86366/249‑3 (MQ=255)                                                         gCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGt                                                                                                                                                                                          <  1:191358/249‑1 (MQ=255)                                                         gCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTAGCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGTGATAAAGATATGGTGCAGCTGGTGACGCCAAATATTACACTTAGC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACg                                                                                                                                                                                              <  1:335329/245‑1 (MQ=255)                                                         gCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAACGGACGACGGTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGGCCGGTGCTGATCAGCACTGGCGGTAAAGATATGGCGCAGCTGGTGACGCCAAATATTACCCTTATC‑AATACCATGACGAATACCCTCCTCGGCACGGAAGAGGGGGTGAATAAGTACGGCgggg                                                                                                                                                                                        >  2:155031/1‑248 (MQ=255)                                                          cGTGCACCAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCAGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGGGACGCCAAATATTACACTTTTCTAAAACC‑TGACGAATACCATCCTCGGACCGGGAGGGGTGGTGAATAAGTACGGCGTgcc                                                                                                                                                                                       >  2:220934/1‑251 (MQ=255)                                                             gCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTgccgcc                                                                                                                                                                                    <  1:32227/251‑1 (MQ=255)                                                               acaAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCAATGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTgccgcc                                                                                                                                                                                    >  1:35943/1‑249 (MQ=255)                                                                       aaCCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGatc                                                                                                                                                                          <  2:310040/251‑1 (MQ=255)                                                                       aaCCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGGGGTGAATAAGTACGGCGTGCCGCCCGAACTGatc                                                                                                                                                                          >  1:242833/1‑251 (MQ=255)                                                                        aCCCTTGCACGCGATGGGTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGttca                                                                                                                                                                         <  2:138448/251‑1 (MQ=255)                                                                              gCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGGGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGGGAATAAGGACGGCGTGCCGCCAGAACTGatccacgatt                                                                                                                                                                   >  2:32300/1‑244 (MQ=255)                                                                                   cgATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTg                                                                                                                                                              <  2:64739/251‑1 (MQ=255)                                                                                    gATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTgg                                                                                                                                                             <  2:68251/251‑1 (MQ=255)                                                                                              gCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGGGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATggg                                                                                                                                                   <  2:331087/251‑1 (MQ=255)                                                                                                     gACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCAACGATTTCCTGGCGCTGATGGGTGActcc                                                                                                                                            >  1:163990/1‑251 (MQ=255)                                                                                                         gCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCctct                                                                                                                                         <  1:356346/250‑1 (MQ=255)                                                                                                                  ggCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTc                                                                                                                               >  1:96396/1‑251 (MQ=255)                                                                                                                  ggCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGGGACTCCTCTGATAACATTc                                                                                                                               >  2:238110/1‑251 (MQ=255)                                                                                                                    cGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCt                                                                                                                             <  1:231664/251‑1 (MQ=255)                                                                                                                     ggTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTg                                                                                                                            <  1:38317/251‑1 (MQ=255)                                                                                                                         tCTGGAGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATGCCATCCTCGGACCGGAAGAGGTGGGGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTTCTGGCGCTGATTGGTGAGTTCTTTGATGACATTCCTGGCg                                                                                                                         <  2:168235/250‑1 (MQ=255)                                                                                                                                tAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCAGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGc                                                                                                                 >  1:208685/1‑251 (MQ=255)                                                                                                                                      cGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGt                                                                                                           >  1:344540/1‑251 (MQ=255)                                                                                                                                                 aTCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCTGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGTCGCTGATGTGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACcgc                                                                                                <  1:213621/251‑1 (MQ=255)                                                                                                                                                   cGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCgcgc                                                                                              >  1:55770/1‑251 (MQ=255)                                                                                                                                                    ggTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCgcgc                                                                                              >  1:210322/1‑250 (MQ=255)                                                                                                                                                    ggTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCa                                                                                             <  2:277322/251‑1 (MQ=255)                                                                                                                                                       aCTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCCGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGc                                                                                          <  2:272240/251‑1 (MQ=255)                                                                                                                                                                    aaGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGt                                                                             <  1:84705/251‑1 (MQ=255)                                                                                                                                                                                 cGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGt                                                                                                               >  2:276260/1‑204 (MQ=255)                                                                                                                                                                                 cGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGt                                                                                                               <  1:276260/204‑1 (MQ=255)                                                                                                                                                                                 cGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGTCTTGGCGGACTgg                                                                >  1:30091/1‑251 (MQ=255)                                                                                                                                                                                      gTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCCAATATTACCCTTATC‑AATACCCTGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGActcc                                                                                                                                            >  1:171555/1‑170 (MQ=255)                                                                                                                                                                                      gTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGTTGAATAAGTACTGCGTGCCGCCAGAACTGATCATCGATTGCCTGGCGCTGATGGGTGActcc                                                                                                                                            <  2:171555/170‑1 (MQ=255)                                                                                                                                                                                                     gCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGTCTTGGCGGACTGGATACGCTGTATGCCGAGcc                                             >  1:37179/1‑250 (MQ=255)                                                                                                                                                                                                      cACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGTCTTGGCGGACGGGATACGCTGTATGCCGAGCCAg                                           <  2:37179/251‑1 (MQ=255)                                                                                                                                                                                                          ggctaTAAAGATATGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGGGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGGCTTGGCGGAGTGGATACGCTGTATGCCGAGCCAGaaaa                                       <  2:316436/247‑1 (MQ=255)                                                                                                                                                                                                                      atGGCGCAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGTCTTGGCGGACTGGATACGCTGTATGCCGAGCCAGAAAAAATTGCTGGGt                            <  2:302259/250‑1 (MQ=255)                                                                                                                                                                                                                          ggcAGCTGGTGACACCAAATATTACCCTTATC‑AATATGCTGACGAATACCATCCTCGGACTGGAAGCGGTGGTGAATGAGTACGGTGTGTCGATAGACCTGATCATCGATTTCCTGGCGTTGATGGGTGACTGCTCTGATAACAGTCCTGGCGTAACTGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAATGTCTTGGTGGACTGGAAGCGCCGTCGGCCGAGTCCGAAAAAATTGTTTGGTTGAga                       <  2:343079/250‑2 (MQ=255)                                                                                                                                                                                                                           gcAGCTGGTGACGCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGTCTTGGCGGACTGGATACGCTGTATGCCGAGCCAGAAAAAATTGCTGGGTTGAGCt                      <  1:263162/251‑1 (MQ=255)                                                                                                                                                                                                                                       gCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGTCTTGGCGGACTGGATACGCTGTATGCCGAGCCAGAAAAAATTGCTGGGTTGAGCTTCCGTGGcgcg           >  1:142656/1‑250 (MQ=255)                                                                                                                                                                                                                                       gCCAAATATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGTCTTGGCGGACTGGATACGCTGTATGCCGAGCCAGAAAAAATTGCTGGGTTGAGCTTCCGTGGcgcg           >  1:222344/1‑250 (MQ=255)                                                                                                                                                                                                                                          aaaTATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGTCTTGGCGGACTGGATACGCTGTATGCCGAGCCAGAAAAAATTGCTGGGTTGAGCTTCCGTGGCGCGaaa        >  1:250219/1‑250 (MQ=255)                                                                                                                                                                                                                                          aaaTATTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGAGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGAGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAAACGCGCAGGCATTGCTGCAAGGGCTTGGCGGACTGGATACGCTGTATGCCGCGCCAGAAAAAATTGCTGGGTTGAGCTTCCGTGGCGCGaaa        >  1:247148/1‑250 (MQ=255)                                                                                                                                                                                                                                            atatTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATtgctgc                                                                                  >  2:362360/1‑174 (MQ=255)                                                                                                                                                                                                                                            atatTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATtgctgc                                                                                  <  1:362360/174‑1 (MQ=255)                                                                                                                                                                                                                                              atTACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGTGCGTCGGTGAAAAAACCGCGCAGGCATTTCTTCAAGGTCTTGGCGGACTGGATACGCTGTATGCCGAGCCAGAAAAAATTGCTGGGTTGAGCTTCCGTGGCGCGAAAACAAt   <  1:211540/251‑1 (MQ=255)                                                                                                                                                                                                                                               ttACACTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGACCTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGTCTTGGCGGACTGGATACGCTGTATGCCGAGCCAGAAAAAATTGCTGGGTTGAGCTTCCGTGGCGCGAAAACAATg  >  2:182429/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        CGTGATGAACTGTTTGAACATTACAAATCACATCGCCCGCCAATGCCGGACGATCTGCGTGCACAAATCGAACCCTTGCACGCGATGGTTAAAGCGATGGGACTGCCGCTGCTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACGCTTATC‑AATACCATGACGAATACCATCCTCGGACCGGAAGAGGTGGTGAATAAGTACGGCGTGCCGCCAGAACTGATCATCGATTTCCTGGCGCTGATGGGTGACTCCTCTGATAACATTCCTGGCGTACCGGGCGTCGGTGAAAAAACCGCGCAGGCATTGCTGCAAGGTCTTGGCGGACTGGATACGCTGTATGCCGAGCCAGAAAAAATTGCTGGGTTGAGCTTCCGTGGCGCGAAAACAATG  >  NC_000913/4047173‑4047661

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