Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,331,890	C→T	L463L (CTC→CTT)	proP →	osmolyte:H(+) symporter ProP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,331,890	0	C	T	95.2%	119.2 / NA	43	L463L (CTC→CTT)	proP	osmolyte:H(+) symporter ProP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (13/27);  minor base G (0/2);  total (13/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 = 1.00e+00

CGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTA  >  NC_000913/4331643‑4332132                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    cGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAAGTCttgt‑c                                                                                                                                                                                                                                                 <  2:319792/251‑1 (MQ=255)           ccGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACg                                                                                                                                                                                                                                       <  1:542655/251‑1 (MQ=255)             ggtccAGCTCGCTTGCGGCGGCATTTAATATTGCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTTGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCTGCGGCGTCAGATAGACAGTAAGCGAATGATATTCTTGT‑CTATCATTACGat                                                                                                                                                                                                                                     <  2:325667/247‑1 (MQ=255)                                                  gCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACgaga                                                                                                                                                                                                <  2:168945/251‑1 (MQ=255)                                                      gTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTTATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGTCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGc                                                                                                                                                                                            <  2:431584/251‑1 (MQ=255)                                                                 gACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGTGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAAGTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAg                                                                                                                                                                                  <  1:488547/250‑1 (MQ=255)                                                                      cAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCTGAATCTGATGATGCCTGCCTATTTCCTGATGGTAGTGGCGGTGGTTGGGGTAAGCACCGGGGTAACCATGGAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATGACTATAATATCGAGCAGAAAATCGATGATAATGACCACGAGATTGCCGATTTGCAGTCGaag                                                                                                                                                                            <  2:30412/251‑2 (MQ=255)                                                                                                    cGCAGAATCTGATGATGCTTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCcgc                                                                                                                                                <  1:126786/249‑1 (MQ=255)                                                                                                    cGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGAtt                                                                                                                                                                                       <  1:283918/210‑1 (MQ=255)                                                                                                                    gcctgccatttcccTGATGGTAGGGGGGGGGGTTGGTTTAATCACAGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGGGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCCCGATATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCt                                                                                                                               <  1:205584/238‑1 (MQ=255)                                                                                                                        gcctATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAACCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAAc                                                                                                                          <  1:204955/251‑1 (MQ=255)                                                                                                                            aTTACCTGATGGTAGGGGCGTTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAAACGTCCGTTGAAAGTTGCGGCTCCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTTACCACGAGATGGCCGATTTGCAGTCGAAACTTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAGACGGa                                                                                                                       <  2:486189/250‑1 (MQ=255)                                                                                                                                ccTGATGGGAGTTGCGGAGGTTGGTTTAATCACCGGTGGAACCATGAACGAGACGGCAAAGCGTCCGTGTAAAGGTGCGACTCCGGCGGCGTCAGATATACAGGATGCGAAGGAAATTCTTGT‑TTAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGAGTGCCGATTTGCGGTCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGcc                                                                                                                  <  1:354131/251‑1 (MQ=255)                                                                                                                                  ggaggggaGTGGCGGTGGTGGGTTTAATCACCGGGGTAACCATGAAAGAGACGGAAAATCGTCCGTTGAAAGGGGCGACACGGGCGGCGGCAGATATAAAGGAAGGGAAGGAAATGCTGGT‑GGAGCATTACGATAATATCGAGAAGAAAATCGATGAGATTGACCACGATAGGTCCGATTTGCAGGCGAAACGTACCCGCCTGTTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTg                                                                                                                <  1:289039/244‑1 (MQ=255)                                                                                                                                   gATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGG‑CGAGCATTACGATAATATCGAGCAAAAAACCGATGATATTGCCCACGAGAGTGCCGATGTTCAGGCAAAACGTCCCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAATCTGAAACGGAGGGCCt                                                                                                                 <  2:510773/249‑1 (MQ=255)                                                                                                                                   gATGGTAGGGGCGGTGGTTGGTTTAATCACCGGCGTAACCAGGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGCCACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGAGGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCt                                                                                                                 <  1:502870/249‑1 (MQ=255)                                                                                                                                       gTAGTGGGGGGGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGGCCGTTGAAAGGGGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGAtgtg                                                                                                           <  2:94380/251‑1 (MQ=255)                                                                                                                                         aGTGGCTGTGGTTGGTTTAAGCACCGGCGTACCCATGAAAGAGCCGGCAAATCGACCGTTGAAAGTTGCGACACCGGGGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATAACGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCAAACCGTCCCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGAc                                                                                                         <  2:85116/251‑1 (MQ=255)                                                                                                                                               ggtggtTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTc                                                                                                   <  1:203793/251‑1 (MQ=255)                                                                                                                                                   gttggttTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTAt                                                                                               >  1:427260/1‑251 (MQ=255)                                                                                                                                                       gttTAACCTCCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTGGAAAGGTGCGACACCGGCGGCGTCAGAAAAACAGGAAGCGAAGGAAATTCT‑GTGCTAGCATTACGAAAATATCGAGCAGAAAACCGATTATATTGACCACGAGATGGCCAATTGGCAGGCGAAACGTACCCGCCTGGGGCAGCAACATCCGCGAATGGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAgg                                                                                           <  1:324543/251‑1 (MQ=255)                                                                                                                                                                ccGGCGTACCCATGAAAGAGACGGCAAATCTGCCGTTGAAAGGTGCGACCCCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCCCGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGaa                                                                                                                            <  1:482371/209‑1 (MQ=255)                                                                                                                                                                ccGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGaa                                                                                                                            >  2:482371/1‑209 (MQ=255)                                                                                                                                                                           aTGAAAGAGAAGGCAAATAGTCCGTTGAAAGGTGAGACACCGGCGGCGTCAGATATACCGGAAGCGAAGGACATTCTTGC‑CGAGCATTAAGATAATATCGAGCGGAAAATCGATGCTATTGCCCACGAGCTTGCCGATTTACAGGCGAAACATCACCGCCTGGGGCAGCACCATCCGCGAATTGATGGATAAGCTGAAACGGATGGCCTGATGTGACGCTGTATTATCAGGCCAATTGCACTCTTAAGGtt                                                                       >  2:477984/1‑251 (MQ=255)                                                                                                                                                                             aaaaGAGGCGGCAAATCGTCCGTTGAAAGGTTCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCa                                                                     <  1:80321/250‑1 (MQ=255)                                                                                                                                                                              aaaGAGACGGCAAATCATCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCAc                                                                    <  1:337983/251‑1 (MQ=255)                                                                                                                                                                               aaGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGCAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACt                                                                   >  2:266132/1‑251 (MQ=255)                                                                                                                                                                                     cGGCAAATCGTCCGTTGAAAGGGGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATACGCGAATTGATGACTAAGCTGAAAAGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGACCTCTTAAGGTTCACTTAAtct                                                             >  2:406397/1‑251 (MQ=255)                                                                                                                                                                                               tCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGCGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCACTTGAACTCTTACGGTTCACTTAATCTCTGACgcgcc                                                   >  2:221760/1‑250 (MQ=255)                                                                                                                                                                                                 cGTTGAAAGGGGAGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCTTTACGATAATATCGAGCAGAAAATAGATGATATTAACCACGAGATTGCCGAGTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATa                                                 <  1:549019/251‑1 (MQ=255)                                                                                                                                                                                                                    ggcggaGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACGCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACgg                              <  2:290392/251‑1 (MQ=255)                                                                                                                                                                                                                         cgTCAGATATACAGGAAGAGAAGGAAATTCTTGT‑CGAGCATTCCGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGCAACGTAACCGCCTGGTGCAGCAACATCAGCAAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCCTACTCTCCTCCAGGTTAACggagg                           >  2:524338/1‑249 (MQ=255)                                                                                                                                                                                                                                 ataCAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGGGGAGAGt                      >  2:339222/1‑246 (MQ=255)                                                                                                                                                                                                                                 agaCAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAACCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGt                      <  1:339222/244‑1 (MQ=255)                                                                                                                                                                                                                                    cAGGACGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTc                                                                      >  2:529567/1‑195 (MQ=255)                                                                                                                                                                                                                                    cAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCTATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTc                                                                      <  1:529567/195‑1 (MQ=255)                                                                                                                                                                                                                                    cAGGAAGCGAAAGAAATTCCTGT‑CGAGCCTTCCGAACCTATCGAGCAGAAAATCGATGATATTGCCCACGAAATTGCCGCTTTGCAGGCGAAACGTACCCGCCTGAAGCAGCAACCTCCGCGAATTGCTGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCCGGCCAATTAAACTCTTATGGATCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGagagcccaatgc               >  2:321001/1‑242 (MQ=255)                                                                                                                                                                                                                                     aGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGaa              >  1:61599/1‑250 (MQ=255)                                                                                                                                                                                                                                     aGGAAGCGAAGGAAATTCTTGT‑CAAGCATTACAATATTACCGAGCAGAAAAGCGATGATATTGACCACGAGATGGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCGTATCAGGCCAATTGAACGCTTAAGGTTCACTTAATCTCTGACGCGCATACGCTCCTCCAGGTTAACGGAGTGGAGTGCAATGaaa             <  2:83866/251‑1 (MQ=255)                                                                                                                                                                                                                                       gAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTGCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGaaaaa           <  2:205539/251‑1 (MQ=255)                                                                                                                                                                                                                                          gCGAAGGAAAGTCTGGC‑CGAGCATTACTATAATATCGAGCAGAAAAGCGATGATATTGACCCCGATATTGCCGATTTGAAGGCGAAACTTACCCGCCTGGGGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGTCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCGCCGCCAGGTTAACTGAGGAGAGTTCAATGAAAAACCg        <  2:362377/251‑1 (MQ=255)                                                                                                                                                                                                                                           cGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCgg       >  1:281982/1‑250 (MQ=255)                                                                                                                                                                                                                                            gAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCgtg      >  2:7142/1‑251 (MQ=255)                                                                                                                                                                                                                                                gAAATTCTTGG‑CGAGCATTACGATAATATCGAGTAGAAAACCGTTGATATTGACCATGAGCTGCCCGATTTGCAGCCGAAAGGTACCCGCCTGGTGCAGCATCATCCTCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCGCCTCCAGGTTAACGGAGGAGAGTGCAATGAATAACCGTGTTTa  <  1:406397/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    CGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTA  >  NC_000913/4331643‑4332132

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