Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,486,205	C→A	T29K (ACG→AAG)	crp →	DNA‑binding transcriptional dual regulator CRP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,486,205	0	C	A	95.6%	146.0 / NA	46	T29K (ACG→AAG)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (25/18);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAA  >  NC_000913/3485963‑3486447                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    taataaATTGATGTACTGAATGCATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGCAGCATATTTCGGCAATCCAGAGACAGCCGCGTTATCTGGATCTGGAGAAAGCTTCTAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCCTGTCTCATTGCCACATTCATAAGTACCAATCCAAGAGCACGCTTATTc                                                                                                                                                                                                                                             >  2:61800/1‑251 (MQ=255)     ccATTGATGTACTGCATTTATGCAAAGGACGTCCCATTACCGTGCAGTCCAGTTGAAAGCCCCTTCCCAGGTAGCGGGAAGCATATTCCGGCAATCAAGAGACAGCGTCGTTATCTGGCTCGGGAGAAAGCTTAAAACACAGAATAACCGCGCTTGGTCCTTGCCACACCCCAATCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCCCATTCATAAGGACCCATCCAAGAGCAAGCTTATTCCCCa                                                                                                                                                                                                                                         <  2:298035/250‑1 (MQ=255)     aCATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCAAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAcc                                                                                                                                                                                                                                          >  1:225327/1‑250 (MQ=255)          gATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTg                                                                                                                                                                                                                                    >  2:193181/1‑251 (MQ=255)             gTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaaa                                                                                                                                                                                                                                 <  1:335181/251‑1 (MQ=255)                     ggaaTGCAACGGACGACACATTACCGTTGATTACAGTCGATCACCGATTCCAAGGTAGCGAGAAGCATATTTCGTGAATCAAGAGAAAGCGGAGCTATCTGGCTCTGGAGAAAGTTTATCACAGAGGATAACCGCGCATGGGGCTTGGCAAACACCAAACAGACCCGACTCTCGAAGGGCTCTTGTCTCATTGCCACATTCATAAGTACCTATCCAAGAGCAGACTTATTCACCAGGGTTAACAA‑GCNGa                                                                                                                                                                                                                         >  2:370018/4‑250 (MQ=255)                          gCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCGTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCATAGACAGGGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCAGGGTGCTGGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCt                                                                                                                                                                                                                   <  1:177166/251‑1 (MQ=255)                           cAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTg                                                                                                                                                                                                                  >  2:29550/1‑251 (MQ=255)                               ggACGTCACATTCCCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGtact                                                                                                                                                                                                              <  2:341094/251‑1 (MQ=255)                                  cgccacaTTACCGTGCATTACAGTTGATACCCCCTTCCCAGGTAGCGGGAAGCAAATTTCGGCATTCCAGAGACAGCGGCGTTATCTGGCTCTGGAAAAAGCTTAGAACAGAGGATAACCGGGCATGGGGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACa                                                                                                                                                                                                           <  2:404369/248‑1 (MQ=255)                                        aTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGtt                                                                                                                                                                                                      >  2:209378/1‑250 (MQ=255)                                               ttcagtacagtaGATCGCCCCTGCCACGGTACCGGGAAGCATATTTCGGCAACCCAGAGACTGCGTCGTTATCTGGCTCTGGATAAAGCTTATACCAGATGATAACCGCGCATGGTGCGTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTGCTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAACGGct                                                                                                                                                                                               <  2:339608/248‑1 (MQ=255)                                                          tGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTg                                                                                                                                                                                    <  1:222860/250‑1 (MQ=255)                                                                    ttCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCaa                                                                                                                                                                            >  1:321440/1‑248 (MQ=255)                                                                               cGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGt                                                                                                                                                               >  1:213409/1‑250 (MQ=255)                                                                               cGGGAAGCATATTTCGGAAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCAGAAACGCTGTACTACATCGTTAAAGTCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTa                                                                                                                                                              <  1:136758/251‑1 (MQ=255)                                                                                            tCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATcc                                                                                                                                                  <  1:66782/250‑1 (MQ=255)                                                                                                  aTCCAGAGACAGCGGCTTTATCTGGCTCTGGAGAAAGCTTCTAACAGAGGATAACCGTGCATGCTGCTTGGCAAACCTCAAACCGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCTGAAACGCTGTACTACATCGTTAAAGTCTCTGTGGCAGTGCTGATCAAAGACGAAGAGTGTAAAGAAATGATCCTCTCCTa                                                                                                                                           <  1:221722/251‑1 (MQ=255)                                                                                                                  gTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTTAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAt                                                                                                                           >  2:71636/1‑251 (MQ=255)                                                                                                                  gTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGACAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGGAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAt                                                                                                                           >  1:275856/1‑251 (MQ=255)                                                                                                                    tATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGAttt                                                                                                                         >  1:239561/1‑251 (MQ=255)                                                                                                                            tctGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCg                                                                                                                 >  1:319594/1‑251 (MQ=255)                                                                                                                                          ataACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGtt                                                                                                    <  2:56690/250‑1 (MQ=255)                                                                                                                                           taACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAACTGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTg                                                                                                  >  2:117547/1‑251 (MQ=255)                                                                                                                                            aaCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTg                                                                                                  >  1:12675/1‑250 (MQ=255)                                                                                                                                            aaCAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACAAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTg                                                                                                  >  1:119259/1‑250 (MQ=255)                                                                                                                                                      aaCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                         >  1:173675/1‑249 (MQ=255)                                                                                                                                                          gcgcATGGTGCTTGGCAAACCGCAAACAGTCCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATc                                                                                                                                                                                                         <  1:385207/133‑1 (MQ=255)                                                                                                                                                          gcgcATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATc                                                                                                                                                                                                         >  2:385207/1‑133 (MQ=255)                                                                                                                                                            gcATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGt                                                                                 >  2:310014/1‑251 (MQ=255)                                                                                                                                                                 gTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAgcgc                                                                            >  1:239785/1‑251 (MQ=255)                                                                                                                                                                      tgacAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATGGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATg                                                                         <  2:404889/246‑1 (MQ=255)                                                                                                                                                                         cAAACCGCAAACAGACCAGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCATGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTAc                                                                    >  1:167459/1‑251 (MQ=255)                                                                                                                                                                             ccGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGc                                                                >  1:317697/1‑251 (MQ=255)                                                                                                                                                                                  aaCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTTGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGaaaa                                                           <  2:225327/251‑1 (MQ=255)                                                                                                                                                                                   aCAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAAc                                                          <  2:328886/251‑1 (MQ=255)                                                                                                                                                                                          cGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCtgt                                                   >  2:66763/1‑251 (MQ=255)                                                                                                                                                                                            aCTCTCGAATGGTTCTTGTCTCATAGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCtgtg                                                  >  2:150733/1‑250 (MQ=255)                                                                                                                                                                                                 cGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTgg                                            >  2:336788/1‑251 (MQ=255)                                                                                                                                                                                                  gAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATACACCAGGGTGAAAAA‑GCGGACACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGc                                           >  1:225005/1‑251 (MQ=255)                                                                                                                                                                                                                tcATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTAc                              >  1:373540/1‑250 (MQ=255)                                                                                                                                                                                                                     gCCACATTCATAAGTCCCCATCCAAGAGCAAGCTTATCCACCAGGGTAAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGTTGATTTTATTGGCGACCTGGGCCTGTTTGAAGAGGGCAAGAAACGTAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaaaa                         <  1:66763/250‑1 (MQ=255)                                                                                                                                                                                                                               tAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAAt              >  1:307544/1‑251 (MQ=255)                                                                                                                                                                                                                                 aGTACCCATCCAAGAGCACGCATCTTCAACCGGGCGAATAC‑GCGGAAACGCTGTACTCCATCGATAACGGCTCTGTGGCCGTGCTGATCAAAGACGCAGAGCGTGAAGAAATGATCCTCTCCTATCTGAATCCGGGTGATTTTCTTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGCACGTGCGAACAACGCTTGTGACGTGGCTGAACTTTCGTACAAAAAAGTTCGCCAATTg            >  2:34429/1‑251 (MQ=255)                                                                                                                                                                                                                                          caaaGAGCAAGCTTATTCACCGGGGTGAAAAA‑GCGGAAACGCTGTACTACAGCGTAAAAGGCTCTGGGGCATTGCTGATAAAAGACGAAGAGGGAAAAGAAAGGTCCCCCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTGGAAGAGGGCAAGGAACGTGGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTa   <  1:193181/249‑1 (MQ=255)                                                                                                                                                                                                                                           aaaGAGCAAGTTTATTCACCAGG‑TGAAAAAAGGGGAAACGCTGTACACCATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGACTCAGGGTGATTTTATGGGCAAACTGGGCCTTTTTGAAGAGGCCAAGGAACGTAGCGCATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTaa  <  1:281738/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    TAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAA‑GCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAA  >  NC_000913/3485963‑3486447

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