Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,263,180	G→T	18.1%	D107E (GAC→GAA)	lolB ←	lipoprotein localization factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,263,180	0	G	T	18.1%	149.5 / 5.4	70	D107E (GAC→GAA)	lolB	lipoprotein localization factor
Reads supporting (aligned to +/- strand):  ref base G (4/52);  new base T (12/2);  total (16/54)
Fisher's exact test for biased strand distribution p-value = 1.36e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAGGTTGCGTTTTGGTGTC‑ATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATACTGATTAAGATTGCGCACGTCTTGCTGATGCTGACGCCATTGTGGCGAATCCGGGCTTTTGC  >  NC_000913/1262952‑1263409                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        tgcgggTTGGGTTTTGGTGTCTTTAACCTAATAAAACATCCTTCCAGTTTTCGTCATTTTGGGTGTAGGGAATTTCGCTCAGGCGGTAGTGGTCGTCCAGTTTGTAGTCGGTTGAATCACCAGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGAAATTGGCATACCGGTAAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAAc                                                                                                                                                                                                                               <  2:572468/235‑1 (MQ=255)    gggTTGCGTTTTGGGTTT‑ATAACCCCCCAAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCa                                                                                                                                                                                                                           <  2:1171445/238‑1 (MQ=255)     gggtGCGTTTTGGTGTT‑ATAACCATCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTTATTTCGCTCAGGCGGTACTGGTCGTACAGTTTGTAGTCGGTTGCATCTCCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTACCTGCa                                                                                                                                                                                                                           <  1:785140/235‑1 (MQ=255)                tttggtactAACCCCCATAACCACCTTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCTAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTGTATTGTC‑GCCTAACTGCACG‑TTACCCGGt                                                                                                                                                                                                               <  2:1011458/231‑1 (MQ=255)                       aTAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCAt                                                                                                                                                                                                        <  1:1129800/239‑1 (MQ=255)                         atcccccttAAACACCTTTCCAGTTTTTGCCGTTTGGGATGTAGGTAATTCAGCTCAGGGCGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTc                                                                                                                                                                                                      <  2:1068261/231‑1 (MQ=255)                              cctttatccacCTTTCACTTTTTGCCTTTCTGGTTGTGGGTAATTTCGCTCAGGCGGTACTGGTCGCCCAGTTTGTAGTCGGTTGCATAACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTc                                                                                                                                                                                                 <  2:354386/230‑1 (MQ=255)                              ccaTAAACATCTTTTCAGTTTTTTCCATTCTGGCTGTAGGTAAGTTCCCCCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTGGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTc                                                                                                                                                                                                 <  2:1019755/239‑1 (MQ=255)                              ccaTAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTc                                                                                                                                                                                                 <  1:1212863/239‑1 (MQ=255)                              ccaTAAAATAGCTTTCATTTTTTGCCCTTTTGGCTGTAGGTTATTTCGCCCAGGGCGTACTGGTCGTCCAGTTTGTGGTCGGTTGCATCAACCGGTAAACCTAAAATCCTCTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTc                                                                                                                                                                                                 <  2:638578/239‑1 (MQ=255)                              cactAAACTACCTTTCGGTTTTTGCCATTTTGGCTGTAGGTAATTCCGGCCAGGCGGTACTGGGCGGCCAGGTTGGAGTCGGTTGCATCACCCGGTAACCCTAAAATCCTCTGGCGCAAGCTGTTTAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTACGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTc                                                                                                                                                                                                 <  2:918056/236‑1 (MQ=255)                                 tAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCCGACCCTTTTTTTT‑GGCTAATCTCACG‑TTTCCCGGGTGGGCAttacgcccc                                                                                                                                                                                                >  1:1145919/1‑230 (MQ=255)                                 tAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGGATCAAGCTTACCTTTCTTTTT‑TACTAACTTCCCA‑TTACCCGGTTAAGCAtttcacaccca                                                                                                                                                                                              >  1:8176/1‑230 (MQ=255)                                   accaacCTTTCATTTTTTGCCATTCTGGCTGTAGGTAATTTCGCGCAGGCGGTACTGGTCGTCCAGTTTTGAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGtt                                                                                                                                                                                            <  2:594297/237‑1 (MQ=255)                                        ccTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGGCCCTTATTTTC‑GCCTAAATGCACG‑TTACCCGGGTGAGCCTTCCGCTCCAGTTCCGt                                                                                                                                                                                        >  1:780596/1‑238 (MQ=255)                                         cTTCCAGTTTTTGCCATTCTTGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGtgc                                                                                                                                                                                      <  1:947764/239‑1 (MQ=255)                                         cTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTCACCCTAAAATCCACTGGCTCAAGCTGTTGAGCGGAATTGGGATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTTTAACGCTTCCCCTTCTTTTC‑GCCCAACTCCACC‑CTCACCCGTTGAGCATTCGTCCCCAGCTCCGtgc                                                                                                                                                                                      >  2:792386/1‑239 (MQ=255)                                         cTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGCCCTTTTTTTTC‑GACTAACTGCACC‑TTACCCGGTTTAGCATTCACGTCCCGTTCCGtgc                                                                                                                                                                                      >  1:103235/1‑239 (MQ=255)                                         cTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGtgc                                                                                                                                                                                      <  1:30587/239‑1 (MQ=255)                                         cTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAAAGCCTACCCTTTTTTTC‑GACTAAACGCACGTTTCCCCGGTTG‑GACTTCAGCTCCCGTTCCGt                                                                                                                                                                                        >  1:594297/1‑237 (MQ=255)                                               gTTTTTTCCATTTGGGGTGTGGGTAATTTCGCCCAGGGGGTACCGGTTGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCACCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCa                                                                                                                                                                                <  2:780596/239‑1 (MQ=255)                                               gTTTTTGTCTTTTTGTCTTTTTGTATTTTCGCTCAGGGGGTCCTGGTCGGTCCGTTTGTAGTCGGTTGCATGACCCGGTAAACCTAAAATCCTCTGGCGCAAGCTTTTGAGGGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCGCGTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCa                                                                                                                                                                                <  2:74386/239‑1 (MQ=255)                                               gTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCa                                                                                                                                                                                <  1:472281/239‑1 (MQ=255)                                                 ttttGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGGCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCa                                                                                                                                                                                <  1:867037/237‑1 (MQ=255)                                                            tGGCTGTAGGTAATTTCGCTCAGGCGGTACGGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGag                                                                                                                                                                   <  1:179485/239‑1 (MQ=255)                                                                   tgtttgtttCGCTGAGGCTGTCATGGTCGTCCAGTTTGTAATTGGGTGCTTCACCCGGTAAACCTAAAATCCCCTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCCGCTCTTCGGCGGCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTGGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGgaga                                                                                                                                                            <  2:971502/233‑1 (MQ=255)                                                                   aGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTGTAAACCTGACCCTTATTTTCTGCCTAA‑TTCCCG‑TTTACCGCTTGAGCATTCACCCCCCGTTCCCGGCTGCCCCCTTGGGTAGAGagccgacga                                                                                                                                                            >  1:102670/1‑233 (MQ=255)                                                                   aGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGa                                                                                                                                                            <  1:1237236/239‑1 (MQ=255)                                                                   aGGTAATTTCGCGCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGa                                                                                                                                                            <  1:881593/239‑1 (MQ=255)                                                                    tgTAATTTCGCCCAGGCGGTTCTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGa                                                                                                                                                            <  1:105356/237‑1 (MQ=255)                                                                                             tcgtcCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAAACCTGACCTTTATTTTC‑GACTAACTGCCCG‑TTACCCGGTTGAGGCTTCAGCTCCCGTTCCGTGCTGCCCAATGGGGTAGTGAGCAGCGGACGGGAGCGGTCCTGGCCGGGTtgct                                                                                                                                   >  1:1025415/1‑238 (MQ=255)                                                                                                 ccAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCCTTATTTTT‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCCATGGGTTAGTGAGGAGCAGCCGGGAGCGATCCTGGCCGGGTTGGTGcccg                                                                                                                              >  1:824484/1‑237 (MQ=255)                                                                                                       tGTAGTTGGTTGCTTTCCCCGGTAAACCTAAAACTTACCGGGGCAGGTTTGGGACCGGAATTTGGATTCCGGTCAATTTGCC‑AATCACCGCTTCGGCGGCATCGTCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAg                                                                                                                        <  2:102670/239‑1 (MQ=255)                                                                                                       tGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAg                                                                                                                        <  1:37278/239‑1 (MQ=255)                                                                                                                ttGCCTCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTa                                                                                                                <  1:829429/238‑1 (MQ=255)                                                                                                                        cccTGTAAATCTAAAATCCCCCGGGGCGAGCTTTTGAGCGGAATTGGCATGCCGGTCAATTGGTC‑AATCATCTCTTCGGCGGCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTttg                                                                                                       <  2:8176/239‑1 (MQ=255)                                                                                                                         ccGGTAAACCTAAAACCCCCTGGCGCACGTTGTTTGGCGGAATTGGCTTTCCGGTCAATTTGCC‑AATCATCTCTTCGGGTGCATCGGCGGTATAACGCTGACCTTTAGTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCACTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGTTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTttgt                                                                                                      <  2:669000/239‑1 (MQ=255)                                                                                                                            gttacccTTAAACCCACTGGCCCAAGTTTTTGAGCGGAATTGGCATTCCGGTAAATTTGCC‑AATCGTCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTgg                                                                                                   <  2:723259/234‑1 (MQ=255)                                                                                                                             tAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGGATAACGCTGACCTTTTTTTTC‑GCCTAAATGCCCG‑TTCCCCGGTTGAGCATTCAGCTCCAGTTCCGGGCTTCCCCCTTGGTTAGTGTGGAGCAAACGGTAGCGACACCTGCCGGGTTTTTGCCCGCAAAAAAGGGGGGCCACTTTTTTTTGGt                                                                                                  >  2:962184/1‑239 (MQ=255)                                                                                                                                cccTAAAACCCCCGGGCGCAAGGTTGTGAGCGGAATTGGCATTTCGGTCAATTTGCC‑AACTCTCTTTTCGGCGGCCTCGGGGGTTAAACGCTGCCCTTTATGTTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAg                                                                                               <  2:1145919/238‑1 (MQ=255)                                                                                                                                aCCTAAAATTTCCTGGGGCAAGGTGTGGGGCGGAATGGGCATGCCGGTCATTTTGCC‑AATCATCTCTTCGGCGGCAGCGGCGGTATAACGCTGACCTTTATTGGC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAg                                                                                               <  2:446924/239‑1 (MQ=255)                                                                                                                                aCCTAAAACCCAGTGCCGTGAGCTGTTGGGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCCTCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTGCACTTTTTGTTGGTCAg                                                                                               <  2:1084413/239‑1 (MQ=255)                                                                                                                                 ccTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAAAGCCCACCTTTCTTTTC‑GCCTAACTGCCCC‑CTTCCCCGCTTAGCATTCAACTCCACTTCCTCGCTGCCCCCTGGGGCTGTGTGCAGGCAACGCTACCGACACCCGCCGGCTTTCTGCCCGAAAACACGGGCGGGCACTTTTTCTTGGTCaca                                                                                              >  2:158729/1‑237 (MQ=255)                                                                                                                                  cTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTTTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTTCGCTCCCGTTCCGTGCTGCCCAATGGGGTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGGTTGCTGCCCGAAAAAACGGGCGGACACTTTTTGTTGGTCAGaa                                                                                             >  1:237825/1‑239 (MQ=255)                                                                                                                                  cTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAg                                                                                                                                                      <  1:371481/182‑1 (MQ=255)                                                                                                                                                       tGTTGGGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGcgccgc                                                                          <  1:323365/237‑1 (MQ=255)                                                                                                                                                       tGTTGGGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGcgccgc                                                                          <  1:1290444/237‑1 (MQ=255)                                                                                                                                                       tGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGcgccgc                                                                          <  1:1018062/237‑1 (MQ=255)                                                                                                                                                       ggTTGAGCGGAATTGGCTTTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGGC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGa                                                                        <  1:892067/238‑1 (MQ=255)                                                                                                                                                              cGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGa                                                                  <  1:768623/238‑1 (MQ=255)                                                                                                                                                               ggAATTGGCTTTCCGGTCAATTTTCCAAATCCTC‑CTTGGGGGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATa                                                                <  2:1239630/239‑1 (MQ=255)                                                                                                                                                               ggAATTGGCATTCCGGTCAATTTGCC‑AATCATCGCTTCGGCGTCATCGGCAGTATAACGTTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATa                                                                <  1:182372/239‑1 (MQ=255)                                                                                                                                                                taagtggatttCCGGTTAATTTGGC‑AGTCTGCCTTTCGGAGTCATCGGCGGTTTATCGCGTCCCTTTATTTGC‑GAGTAACTGCACG‑TGACCCGGTTGAGCATTCAGCTCGAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTGTGATa                                                                <  2:292065/229‑1 (MQ=255)                                                                                                                                                                gAATTGGGATTCCGGGCAATTTGCC‑AATCATTCGTTCGGGGTCATCGGCGGTATAACGCTGACCTTTATTGGC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATa                                                                <  1:960827/238‑1 (MQ=255)                                                                                                                                                                gAATTGGCATTCCGGTCAATTTGCC‑AATCATCTTTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGGC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATa                                                                <  1:407674/238‑1 (MQ=255)                                                                                                                                                                gAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATa                                                                <  1:158729/238‑1 (MQ=255)                                                                                                                                                                gAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATa                                                                <  1:1163719/238‑1 (MQ=255)                                                                                                                                                                gAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATa                                                                <  1:608964/238‑1 (MQ=255)                                                                                                                                                                gAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATa                                                                <  1:341844/238‑1 (MQ=255)                                                                                                                                                                    tGGCATTCCGGTCAATTTGCC‑AATCATCTCTGCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATACTGAt                                                           <  1:1222481/239‑1 (MQ=255)                                                                                                                                                                               tCAATTGGCC‑AATCATCTCTTCGGCGTCAGCGGTGGTATACCGCTGACCTTTATTGTC‑GACTAACTGCACG‑TGCCCCGGTTGAGCAGTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATACTGATTAAGATTgcg                                                 <  1:967400/238‑1 (MQ=255)                                                                                                                                                                                  aTTTTCC‑AATCATCTCTTCGGGGTCTGCGGCGGTTTAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCGCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATACTGATTAAGATTGCGCACg                                             <  2:679866/239‑1 (MQ=255)                                                                                                                                                                                           atcatcTCCTCGGCGGCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATACTGATTAAGATTGCGCACGTCTtgctg                                     <  1:971150/239‑1 (MQ=255)                                                                                                                                                                                                ctctTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCg                                                                                                                      <  1:227982/153‑1 (MQ=255)                                                                                                                                                                                                                         cGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATACTGATTAAGATTGCGCACGTCTTGCTGATGCTGACGCCATTGTGGCGAATCCGGGCt       <  1:862586/239‑1 (MQ=255)                                                                                                                                                                                                                         cGCTGACCTTTATGGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATACTGATTAAGATTGCGCACGTCTTGCTGATGCTGACGCCATTGTGGCGAATCCGGGCt       <  1:1108472/239‑1 (MQ=255)                                                                                                                                                                                                                              aCCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCCATGGGTTAGTGAGCAGCAGACGGTACCGATCCTGGCCGGTTTGCTGCCAGAAAAACCGGGCGTACACTTTTTGTTGTTCAGAAATATAAGCGAACACCCCGCCGGCCTGATATTGATTAAGATTGCGAAAGTCTTGCTGATGACGACTCCATTTGGGGCAATCAGGGCttttt   >  2:800059/1‑237 (MQ=255)                                                                                                                                                                                                                              aCCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTTGTCAGAAATATAAGCGAACGAGCCGCGAGTCTGATACTGATTAAGATTGCACCCGCCTTGCTGAGTCTGGCGCCATTTTGGCGAAAACAGGCTTTTGc  >  2:511659/1‑239 (MQ=255)                                                                                                                                                                                                                              aCCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGTTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATACTGATTAAGa                                                      >  1:713695/1‑187 (MQ=255)                                                                                                                                                                                                                              aCCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATACTGATTAAGATTGCGCACGTCTTGCTGATGCTGACGCCATTGTGGCGAATCCGGGCTTTTg   >  2:792075/1‑238 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        CGCAGGTTGCGTTTTGGTGTC‑ATAACCACCATAAACAACCTTCCAGTTTTTGCCATTCTGGCTGTAGGTAATTTCGCTCAGGCGGTACTGGTCGTCCAGTTTGTAGTCGGTTGCATCACCCGGTAAACCTAAAATCCACTGGCGCAAGCTGTTGAGCGGAATTGGCATTCCGGTCAATTTGCC‑AATCATCTCTTCGGCGTCATCGGCGGTATAACGCTGACCTTTATTGTC‑GACTAACTGCACG‑TTACCCGGTTGAGCATTCAGCTCCAGTTCCGTGCTGCCCAATGGGTTAGTGAGCAGCAGACGGTAGCGATCCTGGCCGGTTTGCTGCCAGAAAAAGCGGGCGTACACTTTTTGTTGGTCAGAAATATAAGCGAACGCGCCGCGAGTCTGATACTGATTAAGATTGCGCACGTCTTGCTGATGCTGACGCCATTGTGGCGAATCCGGGCTTTTGC  >  NC_000913/1262952‑1263409

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