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	97.4%	120.9 / ‑6.5	38	T29K (ACG→AAG)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (22/15);  minor base G (0/1);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCC  >  NC_000913/3485959‑3486451                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         aCAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGGGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCtt                                                                                                                                                                                                                                                      >  2:42874/1‑251 (MQ=255) aCAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGGGCAGTACAGTTGATAGCCCCTCCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTTCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCtt                                                                                                                                                                                                                                                      <  1:372873/251‑1 (MQ=255)   aGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTAt                                                                                                                                                                                                                                                    >  1:317343/1‑251 (MQ=255)    gtactACATTGATGTACTGCATGTATGCAAATGACGTCACATTACCGTGCAGTCCAGTTGATAGCCCCTTCCCCGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACATCGTCGTTATCTGGCTCTGGAGAAAGCTTATCACAGAGGATAACCGCGCATGGTGCGTGGCAAACCGCCAACAGAACCGACTCTCGAATGGTTCTTGTCTCCTTGCCACATTCATACGTACCCATCCACGAGCActcatatt                                                                                                                                                                                                                                                   >  2:239739/5‑244 (MQ=255)                   actgcagttATGCAAAGGTCGTCACATTACCGTGAATTACAGTGGAGAGCCCCTTCCGAGGTAGCGGGAAGCAGATTTCGGCAATCCAGG‑AACATCGGCGTTATCTGGCTCTGGAGAAAGCATATACCGGAGGATAACCGCGCATGGTGCTTGGCAAACCCCAAACAGACCCGACTCTCGAATGGTTCTTGTAACATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGaaaa                                                                                                                                                                                                                                     <  2:460188/242‑1 (MQ=255)                     ggCATGTTTGCAAAGGACGTCACAGTCCCGTNCAGTAATGTTGATAGCCCCTACCCAGTTAGCGGGAGGCATATGTGGGCAATCCAGA‑TACAGCGGCGGTATCTGGCTCTGGAGAAAGCTTATAACGGAGAATAACCGCGCATGGTGCTTGGAAAACCGCAAACAGCCCCGACTCTCGAATGGTTCTTGTCTCATTGCCATATTCATAAGTACCCATCCAAGAGCCAGCTTATGCACCAGGGTGAAAAAGc                                                                                                                                                                                                                                  <  2:218896/250‑1 (MQ=255)                            aTGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAGAAGCGGAAAC‑g                                                                                                                                                                                                                          >  1:218964/1‑251 (MQ=255)                                            aTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGtt                                                                                                                                                                                                           >  2:36473/1‑250 (MQ=255)                                                    gCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGctct                                                                                                                                                                                                  <  1:369918/251‑1 (MQ=255)                                                      aGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGt                                                                                                                                                                                                                      >  2:406314/1‑229 (MQ=255)                                                      aGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGt                                                                                                                                                                                                                      <  1:406314/229‑1 (MQ=255)                                                       gTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAAACGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCAATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCtgtg                                                                                                                                                                                               >  1:438062/1‑251 (MQ=255)                                                         aCAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGc                                                                                                                                                                                             <  1:26515/251‑1 (MQ=255)                                                            gTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAACAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGt                                                                                                                                                                                          >  2:120436/1‑251 (MQ=255)                                                                aTAGCCCCTTCCAAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGATAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACACGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTg                                                                                                                                                                                      >  2:267014/1‑251 (MQ=255)                                                                           ccAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGTCGTTATCGGGCTCTGGAGAAATCTTATAACAGATGATAACCGCGCATGGTGCTTGGCAAACCGCAAACATACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAATGCTCTGTGGCATTGCTGATCAAAGACGa                                                                                                                                                                           <  1:136243/251‑1 (MQ=255)                                                                             aGGTAGCGGGAAGCATATTTCGGTAACCCAGAAGA‑AGCGCCGTTAATTGGCTCTGGAGAAAGCTTATAACAGAGAATAACCGCGCATGGTGCTTGGCAACCCCCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACAGCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAag                                                                                                                                                                         <  1:143872/251‑1 (MQ=255)                                                                                       aaGCATATGTCGGCAACCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑Ga                                                                                                                                                              <  1:184425/251‑1 (MQ=255)                                                                                          caaatTTCGGCAATCCAAG‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCAGGTTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTTCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAAt                                                                                                                                                           <  2:370399/248‑1 (MQ=255)                                                                                          cATATTTCGGCAATCCAGA‑GACAGCGCCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAAt                                                                                                                                                           >  1:362002/1‑251 (MQ=255)                                                                                                 cGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCctc                                                                                                                                                    <  1:458621/251‑1 (MQ=255)                                                                                                  ggCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAAACGCGCATGGTGCTTGGCAAACCGCAAACAGACCAGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCCAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGACGAGGGTAAA‑GAAATGATCctct                                                                                                                                                   >  2:380067/1‑251 (MQ=255)                                                                                                     aaTCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGAAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCt                                                                                                                                                >  2:93534/1‑251 (MQ=255)                                                                                                            a‑gaCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGa                                                                                                                                          >  2:52238/1‑250 (MQ=255)                                                                                                                          aTCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGAttt                                                                                                                             <  1:41660/250‑1 (MQ=255)                                                                                                                               gCTCTGGAGAAAGATTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCCTAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCAGTACTACATCGTTAATGGCTCTGTGGAAGTGGTGATCAACGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATATGAATCAGGGTGATTTTATTgg                                                                                                                       >  2:21799/1‑251 (MQ=255)                                                                                                                                                 aaCAGAGGATAACCGCGCCTGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTg                                                                                                      >  2:393218/1‑250 (MQ=255)                                                                                                                                                           aaCCGCGCATGGTGCTTGACAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGGAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                             >  2:393685/1‑249 (MQ=255)                                                                                                                                                              cgcgCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGCGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAA‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATAAAAGACGAAGAGGGTAAA‑GAAATGAACCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGGGGGCCAGGa                                                                                         >  2:398672/1‑250 (MQ=255)                                                                                                                                                                                aaCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGt                                                                      <  1:59337/251‑1 (MQ=255)                                                                                                                                                                                     caaaccgccccGACTCTCGAATGGTTCTTGTCTCATTGCCACATGCATAAGTACCCATCCAAGAGCAAGCTTAATCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAATGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTTAA‑GAAATGATCCGCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTCTTTGAAGAGGGCCAGGAACGTAGCGCATGAGTACGTGCGaa                                                                 <  2:306238/243‑1 (MQ=255)                                                                                                                                                                                                                  gTCTCATTGACACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGAAGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGGGCGCAAACCGCCTGTGAAGAGGCTGAAATTTCgg                                    >  2:268137/1‑250 (MQ=255)                                                                                                                                                                                                                       aTTGCCACATCAAAAAGGACCCACCAAAGAGCAGGTTTATTCACCAGGGTGAAAAAGCGGAAAC‑GTGGTACTAAATCGTAAAGGGCTCTGGGGAAGTGCTGATAAAAGACGAAGGGGGAAAA‑GAAATTATCCTCTCTTATCTAAATAAGGGTGATTTGATGGGCGAGCGGGGCCTGTTTGAAGAGGGCCAGAAACGTAGCGAATGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaa                               <  1:198465/251‑1 (MQ=255)                                                                                                                                                                                                                         tGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACaaaa                              >  1:145887/1‑250 (MQ=255)                                                                                                                                                                                                                           ccACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCCGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTAGTACAAAAAAt                           >  2:462440/1‑251 (MQ=255)                                                                                                                                                                                                                                       gTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTg                >  1:126396/1‑250 (MQ=255)                                                                                                                                                                                                                                            caaccAAGAGCAAGCGTACTCACCAGGTTGAAAA‑GGAGAAACACCTGTATTCCATCGTTAAAGGCCTTGTGCCGGTGCTGACCAAAGACGAAG‑GGGTAAAAGAAAGGATCCTCTCCAATCTGAAGCAGGGTGATTTTATTGGCGAACTGGGCCGGTTGAAAGTGGGCCAGGAACGTAGCGATTGGGTACGGGCGAAAGCCGCCTGTGAAGTGGCTGAAAGTTCGTACAAAAAATTTCGCCAATTGATTCAg          <  2:36737/248‑1 (MQ=255)                                                                                                                                                                                                                                                 aaGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGAGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTAAGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGAACAAAAAATTTCGCCAATTGATTAAGGTaa      >  1:300287/1‑250 (MQ=255)                                                                                                                                                                                                                                                    agCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGAATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTTGTACAAAAAATTTCGCCAATTGATTCAGGTAAAccc  >  2:84687/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCAGA‑GACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAAC‑GCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAA‑GAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCC  >  NC_000913/3485959‑3486451

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