Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,516	G→A	P170P (CCG→CCA)	creB →	DNA‑binding transcriptional regulator CreB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,516	0	G	A	97.6%	122.8 / NA	43	P170P (CCG→CCA)	creB	DNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/15);  minor base C (1/0);  total (28/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

CGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTT  >  NC_000913/4636270‑4636757                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  cGATCGCCTGCTTGGGCTGGAAATGGTTGCTGACGACTATGTGGCTAAACGTTTTTCACCCCGGGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGcc                                                                                                                                                                                                                                               <  2:443227/251‑1 (MQ=255)       ccggcttggggtgagaatgggTGCTGACGGCTGTGTGGCTAAACCGTTTTCCCCCCGCGAAGGGTGCGCCGGGGTGCGCACCTTACGGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCt                                                                                                                                                                                                                                         <  2:78662/232‑1 (MQ=255)                    caaaTTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGc                                                                                                                                                                                                                             >  1:263835/2‑250 (MQ=255)                         gggTGCGGACGACGATGGGGCAAAACCTTTTTCCCCCCGCGAAGGGTGCGCCAGGGTGCGCACTTTCCTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCCCCCGTCACCCGTATTGGTCATTTTGAATGGAAGGACCCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACtaa                                                                                                                                                                                                                       <  1:90796/250‑3 (MQ=255)                             gCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATgg                                                                                                                                                                                                                    >  1:262699/1‑250 (MQ=255)                                   gcctatgTGGCTAAACCGTTTCCCCCCCGCGAAGGGGGCGCCGGGGTGCGCACCTTACTGCGGCGGGTGAAGAAGTTCCCGTCGCCGTCCCCCGTCACCCGTATTGGTCATTTTGAATTGAATGAACCCGGGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAg                                                                                                                                                                                                                              <  2:8859/228‑1 (MQ=255)                                   gacTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAg                                                                                                                                                                                                                              >  1:8859/1‑234 (MQ=255)                                      aatgtgacTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCGTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCCGCTGGTTTGACACGCCATTAACGCTGACGCGGAATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCTAGCCGCGTCTGGTCCCGCCAGCACCTAATGGCTAGCGtgtg                                                                                                                                                                                                          >  2:352122/8‑248 (MQ=255)                                                            cccGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGCAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGCGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGACGCGTCTGGTCCCGCCAGCCAATAATGGTTAGCGGGTGGGTAGAGGCGCAGGACCCCTAc                                                                                                                                                                                    >  2:259117/1‑251 (MQ=255)                                                                      gtgCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCg                                                                                                                                                                          >  1:272823/1‑251 (MQ=255)                                                                       ggcgcCGGGGTGCGCCCCTTAGTGGGGCGGGTGAAGAAGTTCTGGCCGCGGGCCCCCGTCATCGGTTTGGGTCATTTTGAATTGAATGACCCCGCGGCGAAGATAGGCGGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGt                                                                                                                                                                         <  1:482565/250‑1 (MQ=255)                                                                                           aCTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGa                                                                                                                                                                                                    >  2:337048/1‑204 (MQ=255)                                                                                                           aagTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCg                                                                                                                                                                               >  2:482331/1‑209 (MQ=255)                                                                                                           aagTTCTCGCCGCCGTCCCCCGTCACCCGTATTGGTCATTTTGAATTGAAGGACCCCGCGGCGCAGATCAGCTGGTTGGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCg                                                                                                                                                                               <  1:482331/209‑1 (MQ=255)                                                                                                             gTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCAAAGCTgcgcg                                                                                                                                   >  1:79054/1‑251 (MQ=255)                                                                                                               tctcGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGCTCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTATGGTCCTGCCAGCAACTAATGGATAGCGTGTGGGAGGATGCGCAGGCCACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGcc                                                                                                                                 >  2:75592/1‑251 (MQ=255)                                                                                                                         gTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTACTGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCACGCTGCGCGCCATCACccccc                                                                                                                       >  2:180867/1‑250 (MQ=255)                                                                                                                         gTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGCTACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCg                                                                                                                       >  2:308400/1‑251 (MQ=255)                                                                                                                            tcCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGAGCAGATCAGCTGGTTTGACACGCAATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGAACCGTTGATACCAACATTAAAAAGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATc                                                                                                                    >  1:155630/1‑251 (MQ=255)                                                                                                                                gTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTc                                                                                                                >  2:116036/1‑251 (MQ=255)                                                                                                                                  cATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGCTCTTTCAc                                                                                                              >  2:366035/1‑251 (MQ=255)                                                                                                                                                 ttttGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATc                                                                                               >  1:70199/1‑251 (MQ=255)                                                                                                                                                                ccccggggtgaagatctgttGGTTTGACACGCCATAAACGCTGACGCGGTATGATTTTTTATTGCTGAAGACGTTACCCAAGTCGCCAGGCCGCGTCTGGGCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGAtata                                                                                <  2:70199/232‑1 (MQ=255)                                                                                                                                                                     cggcgCAGGACAGCTGGTTTGACACGCCATTAAAGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGACCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGAGCACTCAACCCCGCTCTTTAACCGATTAATACTCCTCGCGGCATGGGAGATAGccat                                                                           >  2:474038/1‑249 (MQ=255)                                                                                                                                                                     cggcgCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGtt                                                                                                                                                                        <  1:411866/158‑2 (MQ=255)                                                                                                                                                                     cggcgCAGATCAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGtt                                                                                                                                                                        >  2:411866/1‑157 (MQ=255)                                                                                                                                                                           ataTCAGTTGGTTTGACACGCCATTAACGTTGACGCGGTATGAGTTTTTAGTGCTGAAGACGTTACTCAAGTCGCCAGGCTGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGc                                                                     <  1:363699/249‑1 (MQ=255)                                                                                                                                                                               cggttgttttGACCCGCCATTACCGTTGACGGGGTAGGATTTTTTATTGCTGAAGCCTTTACCCAAGCCGCCAGGCCGCGTTTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGAGCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGt                                                                 <  2:452968/244‑1 (MQ=255)                                                                                                                                                                               cAGCTGGTTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCCTCTCGGCCTGGGATATAGCCTGAGGGGCCTGt                                                                 >  2:228925/1‑251 (MQ=255)                                                                                                                                                                                       ttGACACGCCATTAACGTTGTCGCGGTATGAGTTTTAATTGCGGAAGACGTTCCTCAGGTCGCCAGTCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGt                                                          <  1:36137/250‑1 (MQ=255)                                                                                                                                                                                         gcaacGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCAGAAGCTGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATc                                                       <  1:176256/248‑1 (MQ=255)                                                                                                                                                                                         aacactccattaacgcggatgcGGTATGAGTTTTTTGTGCGGAAGCCGTTACTCAAGGCGCCAGGCCGCGTCTGTTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACTCCTACGATCGCACTGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTAt                                                        <  1:301728/229‑1 (MQ=255)                                                                                                                                                                                                         tGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCACGATCTTTCACCGATTAATACTCATAGCGGCATGGGATATAGCCTGCGGGCACGGTAATGCGTATCGGACTGCGGTtgctgc                                       >  2:267636/1‑251 (MQ=255)                                                                                                                                                                                                           acgcGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAg                                                                                                                                          >  1:221834/4‑150 (MQ=255)                                                                                                                                                                                                           acgcGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCGGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAg                                                                                                                                          <  2:221834/147‑1 (MQ=255)                                                                                                                                                                                                                  aTGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGCTATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAtt                              >  2:480911/1‑251 (MQ=255)                                                                                                                                                                                                                    gAGTTTTTATTGCTGAAGACGGTACCCAAGTCGCCAGTCGGGGTCGTGCCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATTCCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAtttt                            <  1:30569/251‑1 (MQ=255)                                                                                                                                                                                                                      gTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGCAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGAATCGGCATGCGGTTGCTGCTGGGCTAtttttt                          >  1:452621/1‑251 (MQ=255)                                                                                                                                                                                                                      gTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCgg                                                     >  1:40646/1‑224 (MQ=255)                                                                                                                                                                                                                      gTTTTTATTGCTGAAGACGTTACTAAAGTAGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAtttttt                          >  1:394760/1‑251 (MQ=255)                                                                                                                                                                                                                       tttttATTGCTGAAGGCGTTCCTCAAGGAGCCCTGCCGCGACTGGTCCAGCCAGCAACTAATGGCTAGCGAGTGGGAAGATGCGCAGGACACCTACGATCGCAACGTTGATAACCACATTAAAACGGTACGTGCCAAGCTGAGCGGCATCAACCCCCATCTTTCACCGATTAATACTCATCGCGGCATGGGATAAAGCCTGGAGGGACTGAAATTCGGACCGGAATGCGGTGGCTGCATGGCTAtttttgg                         >  2:325083/1‑249 (MQ=255)                                                                                                                                                                                                                                       cGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTggtggcggtggcgg         >  1:385251/1‑249 (MQ=255)                                                                                                                                                                                                                                              cAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATACGGTTGCTGCTGGGCTATTTTTTATTggtggaggtggcggaccgggt  >  2:381647/1‑242 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  CGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGGTTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTT  >  NC_000913/4636270‑4636757

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