Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,581,139	G→A	6.3%	C108C (TGC→TGT)	hsdS ←	specificity determinant for hsdM and hsdR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,581,139	0	G	A	6.3%	191.7 / 3.4	63	C108C (TGC→TGT)	hsdS	specificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base A (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.45e-01

CTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATA  >  NC_000913/4580929‑4581384                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      cTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAtt                                                                                                                                                                                       >  2:22037/1‑275 (MQ=255)      ccAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGa                                                                                                                                                                                                                >  4:451864/1‑245 (MQ=255)          cagcGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGAtgtg                                                                                                                                                                                                            <  2:451869/245‑1 (MQ=255)                                             cGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACAt                                                                                                                                                                       >  2:80737/1‑247 (MQ=255)                                                 aaGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCa                                                                                                                                                                   <  4:80736/247‑1 (MQ=255)                                                              ggTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCaaa                                                                                                                                                                                                                             >  2:78984/1‑176 (MQ=255)                                                               gTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCg                                                                                                                                                                                                                                                       >  8:248803/1‑149 (MQ=255)                                                                  tatTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCGGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCaa                                                                                                                                                                                                                                                    >  5:278264/1‑149 (MQ=255)                                                                  tatTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGt                                                                                                                                                                                                                         <  4:78983/176‑1 (MQ=255)                                                                    tttATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCaaaa                                                                                                                                                                                                                                                  >  7:382157/1‑149 (MQ=255)                                                                     ttATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAAt                                                                                                                                                                                                                                                 <  7:174755/149‑1 (MQ=255)                                                                      tATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATg                                                                                                                                                                                                                                                >  7:311173/1‑149 (MQ=255)                                                                       atcaaatcaaAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCACAAAATgc                                                                                                                                                                                                                                               >  8:198126/1‑149 (MQ=255)                                                                          aaatcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACAAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATTTAAGTTTTTCAGGACGTAATACACCGCAAAATGAGcc                                                                                                                                                                                                                                            >  6:8805/1‑149 (MQ=255)                                                                            atcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGa                                                                                                                                                                                                                >  1:10418/1‑175 (MQ=255)                                                                             tcaaaGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGaa                                                                                                                                                                                                                                         >  5:303712/1‑149 (MQ=255)                                                                                    ttGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAt                                                                                                                                                                                                                                  <  7:131720/149‑1 (MQ=255)                                                                                            ttAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCACAAAATGCGCCGAAACTACATTCAAATgg                                                                                                                                                                                                                          <  7:198126/149‑1 (MQ=255)                                                                                                 attattAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCACAAAATGCGCCGAAACTACATTCAAATGGTAGAt                                                                                                                                                                                                                     >  7:81711/1‑149 (MQ=255)                                                                                                     ttaATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTTCGGATTTACCAACTACGGATTTGCTCCCTGATGACAtt                                                                                                                                                                      >  2:373242/1‑192 (MQ=255)                                                                                                           tttGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGc                                                                                                                                                                                                           >  6:395924/1‑149 (MQ=255)                                                                                                             ggCAACAGCCAAAAGTGAGTGAAA‑TTAGTTTGAGTAAAAAGAAGATTTTGTGAAATGAACAATAAAAACAGAAAATAAAAGTTTTACAGGGCGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGGGCGGATTTACCAACTACTGATTTGCTCCCTGATGAGATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAATATTTTTAGGGACAAAAACCAAGTCCGTAGTATCAAACTTGCCAt                                                                           <  4:22037/274‑1 (MQ=255)                                                                                                               cACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTc                                                                         >  1:114610/1‑275 (MQ=255)                                                                                                                 ccAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGCGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGAc                                                                                                                                                                         <  1:66318/177‑1 (MQ=255)                                                                                                                   agcagAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTAc                                                                                                                                                                                                   <  8:382157/149‑1 (MQ=255)                                                                                                                      agAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCaa                                                                                                                                                                                                >  8:67674/1‑149 (MQ=255)                                                                                                                       gAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAAc                                                                                                                                                                                               <  5:350968/149‑1 (MQ=255)                                                                                                                              a‑tgaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCa                                                          <  2:73660/275‑1 (MQ=255)                                                                                                                                tgaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAt                                                                                                                                                                                        >  2:369991/1‑148 (MQ=255)                                                                                                                                 gaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAttt                                                                                                                                                                                      <  7:248803/149‑1 (MQ=255)                                                                                                                                 gaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAttt                                                                                                                                                                                      <  6:134126/149‑1 (MQ=255)                                                                                                                                 gaAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAGTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAttt                                                                                                                                                                                      <  6:278264/149‑1 (MQ=255)                                                                                                                                    aTTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGc                                                                                                                                                                                    <  4:369986/148‑1 (MQ=255)                                                                                                                                     ttttGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTc                                                                                                                                                                                  <  8:96093/149‑1 (MQ=255)                                                                                                                                        tGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCt                                                                                                                                                                               <  8:311173/149‑1 (MQ=255)                                                                                                                                         gTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCtg                                                                                                                                                                              >  7:81714/1‑149 (MQ=255)                                                                                                                                              gATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGAc                                                                                                                                                                         >  6:294185/1‑149 (MQ=255)                                                                                                                                              gATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGAc                                                                                                                                                                         >  6:546/1‑149 (MQ=255)                                                                                                                                               aTAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGAt                                                                                                                                                                                        >  1:69895/1‑133 (MQ=255)                                                                                                                                               aTAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACa                                                                                                                                                                        >  8:17055/1‑149 (MQ=255)                                                                                                                                                tAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACAt                                                                                                                                                                       <  7:127796/149‑1 (MQ=255)                                                                                                                                                   agagAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGc                                                                                                                                                                                    <  3:69895/133‑1 (MQ=255)                                                                                                                                                     agaagaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCaa                                                                                                                                                                  >  5:379800/1‑149 (MQ=255)                                                                                                                                                        agaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATaa                                                                                                                                                               <  6:326054/149‑1 (MQ=255)                                                                                                                                                         gaTTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCACAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAAc                                                                                                                                                              <  8:81711/149‑1 (MQ=255)                                                                                                                                                           ttttGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACaa                                                                                                                                                            >  6:19943/1‑149 (MQ=255)                                                                                                                                                                  aaaTGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTc                                                                                                                                                     <  6:303712/149‑1 (MQ=255)                                                                                                                                                                    aTGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACt                                                                                 >  4:1411/1‑215 (MQ=255)                                                                                                                                                                        gCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGcc                                                                             <  2:1411/215‑1 (MQ=255)                                                                                                                                                                         cAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGc                 >  2:23055/1‑274 (MQ=255)                                                                                                                                                                            tAAAACCAGAAAATTTAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATAtt                                                                                                                                           <  5:8805/149‑1 (MQ=255)                                                                                                                                                                            tAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATgg                                                                                                                                                                                                                          >  2:514078/1‑70 (MQ=255)                                                                                                                                                                                aCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGa                                                                                                                                                                                                                      <  4:514073/70‑1 (MQ=255)                                                                                                                                                                                       aaTATAAATATTTCAGGAAGTAATAAACCGCAAAATGCGCCGAAACTAAATTCCAAAGGTAGAAGCTAAGTTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGGAATAACAATATCTTCAGGAGATAATTTTTGACTTTATTTAACAAGATTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTtata  <  4:23055/275‑1 (MQ=255)                                                                                                                                                                                            aaGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTaa                                                                                                                           <  7:67674/149‑1 (MQ=255)                                                                                                                                                                                               tttttCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACaa                                                                                                                        >  8:370704/1‑149 (MQ=255)                                                                                                                                                                                               tttttCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACaa                                                                                                                        >  8:473969/1‑149 (MQ=255)                                                                                                                                                                                                      ggACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGAttttt                                                                                                                 <  5:294185/149‑1 (MQ=255)                                                                                                                                                                                                      ggACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGAttttt                                                                                                                 <  7:17055/149‑1 (MQ=255)                                                                                                                                                                                                          gTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGa                                                                                                             >  7:380088/1‑149 (MQ=255)                                                                                                                                                                                                           tAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACGTTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGaa                                                                                                            <  8:81714/149‑1 (MQ=255)                                                                                                                                                                                                            aaTACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAAc                                                                                                           >  6:278098/1‑149 (MQ=255)                                                                                                                                                                                                              tACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACaa                                                                                                         >  6:472287/1‑149 (MQ=255)                                                                                                                                                                                                               acacCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACaaa                                                                                                        >  5:199795/1‑149 (MQ=255)                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      CTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGA‑TGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATA  >  NC_000913/4580929‑4581384

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