Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,484,827	C→T	L152F (CTT→TTT)	evgS →	sensor histidine kinase EvgS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,484,827	0	C	T	97.7%	152.2 / ‑3.6	45	L152F (CTT→TTT)	evgS	sensor histidine kinase EvgS
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (19/24);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

TTGCATACCGATTCGCAGCAACGGGTTCGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTCACTCCTTAAATGTA  >  NC_000913/2484584‑2485063                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             ttgCATACCGATTCGCAGCAACGGGTTCGGGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCACACCATTGATAATTACCTTTCCGGCTCTGGTAACCACCTTTCACg                                                                                                                                                                                                                                        >  2:375536/1‑250 (MQ=255)  tgCATACCGATTCGCAGCAACGGGTTCGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGa                                                                                                                                                                                                                                       <  1:1087107/250‑1 (MQ=255)    cATACCGATTCGCAGCAACGGGTTCGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTc                                                                                                                                                                                                                                    >  1:1225133/1‑251 (MQ=255)     aTACCGATTCGCAGCAACGGGTTCGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCa                                                                                                                                                                                                                                   <  1:444767/251‑1 (MQ=255)       aCCGATTCGCAGCAACGGGTTCGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTGGCAGAAGGTGAAGTCGATAGAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCGTTCGCGACTCAAt                                                                                                                                                                                                                                 <  2:551920/251‑1 (MQ=255)           aTTCGCAGCAACGGGTTCGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACAGTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCg                                                                                                                                                                                                                              >  1:1026671/1‑250 (MQ=255)                             gTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCAc                                                                                                                                                                                                                >  1:1089536/1‑246 (MQ=255)                             gTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCAc                                                                                                                                                                                                                <  2:1089536/246‑1 (MQ=255)                              tGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAAc                                                                                                                                                                                                          >  1:69051/1‑251 (MQ=255)                                     aaTGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTaa                                                                                                                                                                                                    <  2:223434/250‑1 (MQ=255)                                         cTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATAtt                                                                                                                                                                                               >  1:486693/1‑251 (MQ=255)                                           gATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATAtt                                                                                                                                                                                               <  1:676799/249‑1 (MQ=255)                                            attattTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCt                                                                                                                                                                                            <  2:486693/251‑1 (MQ=255)                                               attTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCgg                                                                                                                                                                                         >  2:384412/1‑251 (MQ=255)                                                           tAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTa                                                                                                                                                                              <  2:1225133/250‑1 (MQ=255)                                                                    cGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGAc                                                                                                                                                                    <  1:790950/251‑1 (MQ=255)                                                                     gatAAATATCAAATTAACACTGCGGGAATACGCAGATCATCAAAAAGCAATGGATGCGCTTGCAGAAGGAGAAGCCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTAACTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACg                                                                                                                                                                   >  2:563472/3‑251 (MQ=255)                                                                        aanTATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAg                                                                                                                                                                 <  2:1111718/250‑1 (MQ=255)                                                                          atatCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTa                                                                                                                                                              <  2:863406/251‑1 (MQ=255)                                                                                    aaCACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTGCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCTCTTCCCTCACCAAAACCAGGAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCatcaat                                                                                                                                                    <  2:189925/251‑1 (MQ=255)                                                                                     acacTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATtcat                                                                                                                                                        >  1:1243283/1‑246 (MQ=255)                                                                                     acacTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATtcat                                                                                                                                                        <  2:1243283/246‑1 (MQ=255)                                                                                          ccGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTc                                                                                                                                              >  2:446219/1‑251 (MQ=255)                                                                                                     cAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGcaaca                                                                                                                                   <  1:97047/251‑1 (MQ=255)                                                                                                             caAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAtct                                                                                                                            <  2:1026671/250‑1 (MQ=255)                                                                                                                   gCAATGGTCGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTGCGCCGCCTCTTAATAAGGACAGTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTGCCTCACCAAAAGCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGATGTAATTGATCAATCATTTCCAAAAGCAACAAATATCTCGTTTAc                                                                                                                     <  2:718557/251‑1 (MQ=255)                                                                                                                         gACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAAttt                                                                                                               <  1:400004/251‑1 (MQ=255)                                                                                                                            gcgcTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTAGCTCTTTTACAAATTtata                                                                                                            <  1:903180/251‑1 (MQ=255)                                                                                                                                     gAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCAtt                                                                                                   >  1:732245/1‑251 (MQ=255)                                                                                                                                            aaGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCttt                                                                                                                        <  1:1016386/223‑1 (MQ=255)                                                                                                                                            aaGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCttt                                                                                                                        >  2:1016386/1‑223 (MQ=255)                                                                                                                                                               acaTTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCaaaa                                                                                                                                         >  2:364497/1‑187 (MQ=255)                                                                                                                                                               acaTTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCaaaa                                                                                                                                         <  1:364497/187‑1 (MQ=255)                                                                                                                                                                                  ccTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATtcatca                                                                                                                                                      >  1:472561/1‑155 (MQ=255)                                                                                                                                                                                  ccTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATtcatca                                                                                                                                                      <  2:472561/155‑1 (MQ=255)                                                                                                                                                                                        aataatGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCAGAATGATTACTTTATTGGTAGTAAcatcat                                                <  1:384412/251‑1 (MQ=255)                                                                                                                                                                                                cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACcagc                                         >  2:684753/1‑250 (MQ=255)                                                                                                                                                                                                cATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATGAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACcagca                                        <  1:446219/251‑1 (MQ=255)                                                                                                                                                                                                                  aTTGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTAtt                      >  1:918576/1‑251 (MQ=255)                                                                                                                                                                                                                   ttGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTAt                                    <  1:509907/236‑1 (MQ=255)                                                                                                                                                                                                                   ttGATAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTAt                                    >  2:509907/1‑236 (MQ=255)                                                                                                                                                                                                                       tAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTtcact                 >  1:465808/1‑251 (MQ=255)                                                                                                                                                                                                                       tAATTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTtcact                 >  2:75130/1‑251 (MQ=255)                                                                                                                                                                                                                         aTTACCTTTCCGGCGCTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTtcactca               <  2:69051/251‑1 (MQ=255)                                                                                                                                                                                                                                      cgcTGGTAACCACCTTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTATTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTCACTCCTTAAATGTa  <  1:684753/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             TTGCATACCGATTCGCAGCAACGGGTTCGTGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCTCTTAATAATGACATTGCTGCAACCAAACCATTGATAATTACCTTTCCGGCGCTGGTAACCACCCTTCACGACTCAATGCGACCGCTTACCTCACCAAAACCAGTAAATATTGCTCGGGTAGCAAATTACCCCCCAGACGAGGTAATTCATCAATCATTTCCAAAAGCAACAATTATCTCTTTTACAAATTTATATCAGGCATTAGCATCCGTCTCAGCTGGGCACAATGATTACTTTATTGGTAGTAACATCATTACCAGCAGTATGATTTCCCGCTATTTCACTCACTCCTTAAATGTA  >  NC_000913/2484584‑2485063

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