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	97.6%	119.3 / ‑4.5	41	L463L (CTC→CTT)	proP	osmolyte:H(+) symporter ProP
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (9/31);  total (9/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

ATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGA  >  NC_000913/4331649‑4332135                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       atatCCGTTACAGCGCGCGGGCGGCGGCATTTAATAGTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCAt                                                                                                                                                                                                                                               <  2:457563/251‑1 (MQ=255)         tACAGGGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTTACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGGAGTGGCGGTGGTTGGTTTAAGCACCGGCGTAACCATGAAAGAGACGGCAAAGCGGCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGata                                                                                                                                                                                                                                        <  2:469794/250‑1 (MQ=255)                   tGGCGTCGGCAGTTAAGAGTGCGGTGCTGGTTTCCGGTCTGACGCCTACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGATTCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGGTTAAACACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGc                                                                                                                                                                                                                               <  2:1552/249‑1 (MQ=255)                    ggcggcggcATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGa                                                                                                                                                                                                                            <  1:430764/251‑1 (MQ=255)                               ttAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTAc                                                                                                                                                                                                                                            <  2:238062/224‑1 (MQ=255)                                      ttCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTTGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGAc                                                                                                                                                                                                           <  1:469617/250‑1 (MQ=255)                                               ggTTGCCGGTCTGACGCCAACGCTGGCGTCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGAt                                                                                                                                                                                                   <  1:113526/249‑1 (MQ=255)                                                gTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGc                                                                                                                                                                                                <  1:155765/251‑1 (MQ=255)                                                   gCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCg                                                                                                                                                                                              <  2:286862/250‑1 (MQ=255)                                                      ggTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGAttt                                                                                                                                                                                          <  1:366029/251‑1 (MQ=255)                                                           gACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCCCGCAGAAACTGATGATGCCTGCCTATTACCTGATGTTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCa                                                                                                                                                                                       <  2:63976/249‑1 (MQ=255)                                                                   cGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTTGTTGGGTTAATCGCCGGCGTATCCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGGCATATATACAGGAATCGAAGGAAATTCTTGGCGAGCATTGCGATAATATCGAGCAGAAAATCTATGATATTGGCCACGAGATTGCCGATTTGCAGGCGAAACGt                                                                                                                                                                             <  1:21595/251‑1 (MQ=255)                                                                   cGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGt                                                                                                                                                                             <  2:79135/251‑1 (MQ=255)                                                                            ccTGGCTGTCCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTCCCTGATGGTAGTGGCGGTGGTTGGTTTAATCCCCGGTGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTg                                                                                                                                                                    <  1:260853/251‑1 (MQ=255)                                                                                   ggTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTgcagc                                                                                                                                                             <  1:273832/251‑1 (MQ=255)                                                                                        aaagtccGCAGAATCTGATGATGCCTGCTTTTTCCCTGAGGGTATTGGCGTGGGTGGGTTTAATCCCCGGGGTAACCATGAAAGAGACGGCAAATCGTCTGTTGAAAGGGGCGACACCGGGGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCGGGCGAAACGTACCCGCCTGGTGCAGCAACAt                                                                                                                                                        <  2:479860/245‑1 (MQ=255)                                                                                          aGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGaaa                                                                                                                                                                                                                          <  1:155610/183‑1 (MQ=255)                                                                                                                  gcctATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGAAAATCGGCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATcc                                                                                                                                                      <  2:429162/227‑1 (MQ=255)                                                                                                                                        gggtggtTGGTTTAATCACCGGCGTAACCATGAAAGAGACGTCAAATCGTCCGTTGAAAGGGGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGt                                                                                                        <  1:328786/250‑1 (MQ=255)                                                                                                                                                ggttTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATc                                                                                                 >  1:310226/1‑249 (MQ=255)                                                                                                                                                   ttAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGcc                                                                                            >  1:146786/1‑251 (MQ=255)                                                                                                                                                    tAATCACCGGCTAAACCATGAAAGAGACGACAAAGCGTCCGTGGAAAGGTGCGACACGGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTTTTGTCGAGCATTACGATAATATCGAGCATAAAATCGATGATATTGACCACGATATTGCCGATTTGCGGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGTTGTGCCGCTGT‑CTTATCAGGct                                                                                            <  2:69955/250‑2 (MQ=255)                                                                                                                                                     aaTCACCGGCGTAACCATGAAAGAGACGGCAAATGGTCCGTTGAAAGGGGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTTCCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCa                                                                                           <  1:277399/250‑1 (MQ=255)                                                                                                                                                       ccACCGGCGTAACCATGAAAGAGAGGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAAtt                                                                                        <  2:211016/250‑1 (MQ=255)                                                                                                                                                         ccctgtgtacccaggaaagaacgggCAAATCGTCCGTTGAAGTGTGCGGCACCGGCGGCGCCAGTTATACAGGAAGGGAAGAAAATTTTTGGCCGGCATTACGATAATACCGAGCAGAAAATCGATGATATTGACGACGATATTTCGGTTTTGCGGGAGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGa                                                                                      <  2:415933/228‑1 (MQ=255)                                                                                                                                                                tAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGAGGGCCTGATGGGCCGCTGT‑CTTATCAGGCCAATTGAACTCtt                                                                                >  1:273126/1‑250 (MQ=255)                                                                                                                                                                 aaCCATAAAAGATACGACAAACCGTCCGTTGAAAGGTGCGTCACCGGCGGCGTCACATATAAAGGAAGCGAAGAAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGTCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGAACTCTTa                                                                               <  1:264746/250‑1 (MQ=255)                                                                                                                                                                    cATGAAAGAGACGGCAAACCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGAACTCTTAAGGt                                                                           <  1:286335/251‑1 (MQ=255)                                                                                                                                                                      tGAAAGAGACGGCAAATCGTCCGTTGAAAGGAGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGATAATTACGATAATATCGCGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGACGAAACGTACCAGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGC‑GTGCTTATCAGGCCAATTGAACTCTTAAGGtt                                                                          >  2:45193/1‑250 (MQ=255)                                                                                                                                                                                  cAAATCGTCCGTTGAAAGGTGCGCCACCGGCGGCGTCAGATATACAGGCAGCGAAGGAAATTCTTGTCGCGCATTCCGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGGACCCGCCTGGTGCAGCAACATCCGCGCATTGATGACTAAGCTGAAACGGATGGCCTGCAGTGACGCTGT‑CTTATCAGGCCAATTGATCTCTTATGGTTCACTTaaactttg                                                             >  2:244632/1‑245 (MQ=255)                                                                                                                                                                                          ccGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCAt                                                     >  2:11886/1‑251 (MQ=255)                                                                                                                                                                                                 aaGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTctcc                                              >  1:231365/1‑251 (MQ=255)                                                                                                                                                                                                  aGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGAt                                                                                                                  >  2:51019/1‑183 (MQ=255)                                                                                                                                                                                                  aGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGAt                                                                                                                  <  1:51019/183‑1 (MQ=255)                                                                                                                                                                                                        gACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTTAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACCCTCCTCCAGGt                                       >  2:26072/1‑251 (MQ=255)                                                                                                                                                                                                          cacCGGCGGCGTCAAATATACAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTa                                     <  1:411211/251‑1 (MQ=255)                                                                                                                                                                                                                             aCAGGAAGCGAAGGAAATTCTTGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATg                   <  1:58654/250‑1 (MQ=255)                                                                                                                                                                                                                                      gaagaaAATTCTTGTGGGGCATTACGAAAATATCGGCCAGAAAAGCGAGGATATTGACCACGAGATGGCCGTTTTGCAGGCGAACCGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCgtg         <  2:195064/251‑1 (MQ=255)                                                                                                                                                                                                                                      aaaGGAAATTCTTGTCGAGAATTACGATAATACCGAGCAAAAAATCGATGATATGGACCACGAGATTGCCGATTTGCAGGCGAAACTTCCCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCgtg         <  2:341743/250‑1 (MQ=255)                                                                                                                                                                                                                                       aagaaaaTTCTTGTCGAGCATTACGATAATTTCGAGAAGAAAACCGATGATATTGCCCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCgtgt        <  1:101031/247‑1 (MQ=255)                                                                                                                                                                                                                                             aTTCTTGTCGAGCATTACGATAAGATCGAGCAAAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGa  <  2:310226/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       ATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTCGAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGT‑CTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGA  >  NC_000913/4331649‑4332135

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