Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,630,553	A→G	intergenic (‑31/‑180)	ettA ← / → slt	energy‑dependent translational throttle protein EttA/soluble lytic murein transglycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,630,553	0	A	G	95.7%	143.8 / ‑2.4	48	intergenic (‑31/‑180)	ettA/slt	energy‑dependent translational throttle protein EttA/soluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (22/23);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

GCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCT  >  NC_000913/4630335‑4630788                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            gcTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCttttt                                                                                                                                                                                                             >  1:499609/1‑251 (MQ=255)   ttCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTcc                                                                                                                                                                                                           >  1:442337/1‑251 (MQ=255)     ccccTTCGATGTTTTTATCAATTCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTccc                                                                                                                                                                                                          <  1:270314/248‑1 (MQ=255)          tCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGACGTGTAGCCTTTTTCCCCGCCt                                                                                                                                                                                                     >  2:174816/1‑250 (MQ=255)               gTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGCCCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCCTGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTTATGATGAGTTTTCCAAGGCGAAATGAAGCCTTTTTCCCCGCCTAATccc                                                                                                                                                                                               >  2:187957/1‑251 (MQ=255)               gTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATccc                                                                                                                                                                                               >  1:183395/1‑251 (MQ=255)                   tcacccagccccGCCAAAATTCTCATCAGGGTAGACTTCCCCGCGCCGTTCAGACCCAGGACCCCAATTTTGGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTTCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTTTGTTTATCGTTATGATGAGTTTTCAAAGGCGAAGTGTCGCCTTTTTCCCCGCCTAATCCCAGc                                                                                                                                                                                            <  2:67634/241‑1 (MQ=255)                    tatcaagcccctccATAATGCGCAGCAGTGTAGATTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTCCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAgtt                                                                                                                                                                                           <  1:447576/238‑3 (MQ=255)                       cAATGCCCGCCATAATGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATccccgttggtc                                                                                                                                                                                       >  2:156202/1‑243 (MQ=255)                          cccccGCCATAATGCTCAGCAGGGTAGACATACTCGCGCCGTTCAGACCCAGGACACAAATTTTTGCCCCAGGGAAGAAAACCAGAGAGATGTTTTT‑CAAAATATGACGTGTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGAt                                                                                                                                                                                    <  2:356292/249‑1 (MQ=255)                           ccccGCCATAATGCGCAGCAGGGGAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATAGGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTGTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATc                                                                                                                                                                                   <  1:57670/250‑1 (MQ=255)                                 cATAATGCGCAGCAGGGTAGACTTACTCGCGCCGTTCAGACCCAGCACACCAATTTTTGCCCCAGGGCAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGCACCACTTTGCCGACACGATGCATGGTCTAAACGAATTGAGCCACGTTGGACTTCGCCTCTCTGTTTATCGTGATGATGGGTTTTCACAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATccacgttggtcgttcactttc                                                                                                                                                                             >  2:62715/1‑232 (MQ=255)                                      tGCGCAGCAGGGTAGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCagt                                                                                                                                                                        >  1:184363/1‑251 (MQ=255)                                                                            aggaaacaaaTTTTGGCCCCAGGGAAAAAACTCAGGGAGA‑GTTTTTGTAAAAGTTGCAGTTTCGGCGGAACATCTTTGCCTACACGATGCATGGTATAAACGAATTGAGCCTCGTGGGACTCCGCCTCTTGGTTTATCGTGATGATGAGTTTGCAAAGGCGAAGTGTAGCATTTTACCCCGCCTAATCCCAGTTGGTCGAGCACCTTCGCAGTAAAAGTAAAAAAGGGGCCGTAACGTGGCGTAATCGGCa                                                                                                                                  <  2:44867/243‑1 (MQ=255)                                                                               acacCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAAAGAATTGAGCCAAGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAACCGGCAATg                                                                                                                               >  2:412175/1‑251 (MQ=255)                                                                                           gCCCCAGGGAAGAAACTCAGAGATATGTTTTT‑AAAAATATGACGTTTTGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGCCTGCGCCTCTATGTTTATCGTGAGGATGAGGTTTCAAAGGCGAAGTGTAGCATTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCa                                                                                                                    <  2:342631/250‑1 (MQ=255)                                                                                           gCCCCAGGGAAGAAACTCAAAGAGAGGTTTTT‑CAAAAGATGCCGTTTCGGCGGAACAACTTTGCCGACACGATGCATTGTATAAACGTATTGATCCACGTTGGACTTCGTCTCTATGTTTATCGTGATGATGAGTTTGCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGGCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGcc                                                                                                                    <  2:408631/250‑2 (MQ=255)                                                                                             cccAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATaa                                                                                                                 <  1:162622/251‑1 (MQ=255)                                                                                              ccAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACACCTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGGTGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCAGCCTAATCCCAGTTGGTCGATCACCTTCGCTGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATaa                                                                                                                 >  2:262683/1‑250 (MQ=255)                                                                                                aGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGCCGTTTCGGCGGCACACCTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGCTGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAACGTGTCAGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATCAATc                                                                                                              >  2:32720/1‑251 (MQ=255)                                                                                                   gaagaaACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTAt                                                                                                           >  1:133864/1‑251 (MQ=255)                                                                                                   gaagaaACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTAt                                                                                                           <  1:463027/251‑1 (MQ=255)                                                                                                     agaaACTCAGAGAGATGTTTTT‑AAAAATATGACGTTTCGGCGGAACAACTTTGCAGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCACTATGTTTATCGGGATGATGAGATTTCAAAGGCGAAGCGTAGCCTTTTTCCCCGCCTAATCACAGTTGGACGATCACCTTCGAAGTAAAAGTAAAAAAGGGGCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATaaagatatt                                                                                                          >  2:229383/1‑244 (MQ=255)                                                                                                         acccagaggtagtttttt‑CAAAATATCGCGTTCCGGCGGAACAATTTTGCCGACACGATTCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGGGATGATGAGTTTTCAAAGGCGAAGTGTAGCTTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTAcgcg                                                                                                     <  2:409157/238‑1 (MQ=255)                                                                                                                 agctGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGATGGACTTCGCCTCTATGTTTATCGAGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCAACTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCAATTACGCGGCCTGACg                                                                                             >  1:455632/4‑251 (MQ=255)                                                                                                                 agaTGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTAtt                                                                                                          <  1:34384/238‑1 (MQ=255)                                                                                                                 agaTGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTAtt                                                                                                          >  2:34384/1‑238 (MQ=255)                                                                                                                 agaTGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACg                                                                                             >  1:455641/1‑251 (MQ=255)                                                                                                                               aaTATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGt                                                                                 >  2:420969/1‑250 (MQ=255)                                                                                                                                  atGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTAttt                                                                             >  1:267482/1‑251 (MQ=255)                                                                                                                                  atGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTAtt                                                                              <  1:46262/250‑1 (MQ=255)                                                                                                                                  atGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTAtt                                                                              <  1:174816/250‑1 (MQ=255)                                                                                                                                    gacgtgttgggggaaaaatgtggtcGACACGATGCATGGTATAAACGAATAGACCCTCGTTGGAGTGCGCCTCTATTTTTTTCGTGATGATGAGTTTTCAAAGGCGAATTGTAGGATTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTTTCTGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCATCATGACTCTGCATTGATGCATTTac                                                                           <  2:455378/227‑1 (MQ=255)                                                                                                                                    gACGTTTCGGCGGAACAACTTTGCCGACACGATGCTTGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTa                                                                            <  1:438606/250‑1 (MQ=255)                                                                                                                                                 aacaacTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGAATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATgc                                                              >  2:270817/1‑251 (MQ=255)                                                                                                                                                           ccGACACGATGAGTGGTATAAACGTGTTGAGCCACGTTGGACTTCGCCTCTATGTTTATAGAGATGATGAGTTTTCAAAGGCGAAGTGTAGCATTTTTCCCCGCCTAATCCCAGTTGGTCGATCACATTCGCAGTAAAAGAACAAAAGTGTCCGTAACGTATCGTAAACGGCAATGACTGGTGAGCATAAATATATTACGCGGCCTGACGCCGCATTGATGTATTTACCCTTAGAGGATgcgcgtgcgga                                                     >  2:91667/1‑243 (MQ=255)                                                                                                                                                               cacGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTc                                                                                                                                                                             >  1:126371/1‑126 (MQ=255)                                                                                                                                                               cacGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAAGGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTc                                                                                                                                                                             <  2:126371/126‑1 (MQ=255)                                                                                                                                                                 cGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATTTTTAGCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGcc                                              <  1:183563/251‑1 (MQ=255)                                                                                                                                                                    ggcttGGTATAAACGAATGTAGCAACTTTGGACTTCGCCGATAGGTTTATCGTGAGTAGGAGTTTTAAAAGGCGAAGTGTAGCCTTTTTCAACGCCGAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATGGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCaaa                                           <  2:206908/247‑1 (MQ=255)                                                                                                                                                                             taAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAAAAAGTTAcc                                  >  2:166143/1‑251 (MQ=255)                                                                                                                                                                              aaaCGAATTGGGCCACGTTGGACGTCGCCGCGATGTTTATCGGGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCt                                 <  2:183395/251‑1 (MQ=255)                                                                                                                                                                               aaCGAATTGAGCCACGTTGGACTGCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTg                                <  1:49148/251‑1 (MQ=255)                                                                                                                                                                                  gAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTggcg                             <  1:440184/251‑1 (MQ=255)                                                                                                                                                                                   aaTTGAGCCACGTTGGACTTCGCTTCTATTTTTATCGTGTTGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTGACCTggc                              <  1:318526/249‑1 (MQ=255)                                                                                                                                                                                           cACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGGGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAACCGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTGCACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTAGGc                    >  2:220361/1‑251 (MQ=255)                                                                                                                                                                                                        gCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGtgtctgt       <  1:410233/251‑1 (MQ=255)                                                                                                                                                                                                              aTGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTctgct  <  1:156202/250‑1 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            GCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGTGGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTT‑CAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCT  >  NC_000913/4630335‑4630788

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: ‑