Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,717	T→C	L8L (TTG→CTG)	creC →	sensory histidine kinase CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,717	0	T	C	94.6%	105.1 / ‑3.9	37	L8L (TTG→CTG)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (16/19);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

TCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAA  >  NC_000913/4636474‑4636959                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    ggggTATGATTTTTTATTGCTGAAACCGTTTCTCAAGTCGCCAGGCCGCGTTTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTggg                                                                                                                                                                                                                                              <  1:58798/249‑1 (MQ=255)    gTATGAGTTTTTATTGCTGAAGCCGTTCCTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTa                                                                                                                                                                                                                                           <  1:83627/251‑1 (MQ=255)             tttATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACt                                                                                                                                                                                                                                  >  1:271535/1‑251 (MQ=255)                       aggacgttcttcaagcggccaggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTggtggcggtg                                                                                                                                                                                                                        <  1:163376/230‑1 (MQ=255)                            gTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTggtggcggtggcggc                                                                                                                                                                                                                   >  1:83872/1‑248 (MQ=255)                             ttACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGcc                                                                                                                                                                                                                  >  1:30696/1‑251 (MQ=255)                               attcaaatcgccaggctgcgTCTGGCCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGTCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTg                                                                                                                                                                                                                <  1:263395/234‑1 (MQ=255)                                         ccaggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGGTCGTACTg                                                                                                                                                                                                      >  1:173635/4‑251 (MQ=255)                                          caggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACCCCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGGCTGGTTCGGACTgg                                                                                                                                                                                                     >  2:80363/3‑251 (MQ=255)                                          caggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTgg                                                                                                                                                                                                     <  1:231802/249‑1 (MQ=255)                                                               aGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTa                                                                                                                                                                               <  1:80363/251‑1 (MQ=255)                                                                 cAACTGATGTATAGCATATGGGAAGACGCGCAGTACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGGCAAGCTGCGCGCCAGCAGCCCCGATCTTTCCCCGAATAATACTCAACGCTGCATGGGATATATTCTGAGGGGCCTGTAATGGTTATCCGCAGGCGGTTGTTGCTGGGCTAGTTGTTACTGGTGGCGGGGGCAGCCGGGTTCTTACTGGCCATTTTTGTCAAAG‑AAGTTaaa                                                                                                                                                                             <  1:15879/251‑1 (MQ=255)                                                                         ggATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCgg                                                                                                                                                                     >  1:137866/1‑250 (MQ=255)                                                                                   tGGGAAGATGCGCAGGACAACTACGATCGCACCGATGATACCCCCATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGAATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTTAAAGTAAGT‑AAACCGGGCGTGCGAGGAGc                                                                                                                                                            >  2:253636/1‑250 (MQ=255)                                                                                             cgcAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAggggg                                                                                                                                                 >  1:250591/1‑250 (MQ=255)                                                                                              gcAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAtttttt                                                                                                                                                                                                                                     >  1:253586/1‑167 (MQ=255)                                                                                              gcAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTAtttttt                                                                                                                                                                                                                                     <  2:253586/167‑1 (MQ=255)                                                                                                     ccctccgaccgccccgtgtaTACCCACATTAAAACGTGGCGTGCCAAGTTGCGCGCCATCAACCCCGATCTGTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGTGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGTGGACGTTAATc                                                                                                                                         <  2:250591/232‑1 (MQ=255)                                                                                                                      tgATACCCACATTAAAACG‑TTCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTg                                                                                                                       >  1:74948/2‑251 (MQ=255)                                                                                                                                  tAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGCTGTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTgg                                                                                                             <  2:67123/250‑1 (MQ=255)                                                                                                                                     aaCGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGcgcg                                                                                                         >  1:245171/1‑251 (MQ=255)                                                                                                                                            cGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTTGTACTGGCTATTTTTGTCAAAG‑AAGTTTAAACGGGCGTGCGAAGGGCAACGGAGGGGGCGGTAAATGATACCGCCACGTTGCTGGGGGGGCTGGggcgggccgct                                                                                                  >  2:364210/1‑244 (MQ=255)                                                                                                                                                       ccgcgcCATCAACCCCGATCTTTCACCGTTTAATACTCATCGCGGCATGGGATACAGCCTGAGGGGCCTGTAATGCGTAGCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTGCGTACTGGCTATTGTTGTCAAAG‑AAGTTAAACCGGTCATGCTAAGAGCAACGGATGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGctctct                                                                                         <  1:216340/248‑1 (MQ=255)                                                                                                                                                                aaaacccGATCTGTCACCTAGTAATTCTATTCGCGGCATGGGATATAGCCTGAGGGGCCTTGAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTGTTACTGGTTGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGGCAAAG‑AAGGTAAACCGGGCGGGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAAc                                                                              <  2:302516/247‑1 (MQ=255)                                                                                                                                                                              cACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCTTGTAATGCGTATCGGCATGCGGTTGCTGCTTGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTg                                                                 <  1:379396/250‑1 (MQ=255)                                                                                                                                                                                         gttgaTCGCGGTTTGGGATATAGTCTGGGGTGGCTGTAATGCGTATCGGCAGTCGGTTGCTGCCGGGATAGTTTTGACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGGCCAAG‑AAGTTAAACCGGGCGGGCGAAGAGCAACTGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCGGGCGTGTCCCGAGTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGAGCTGGCGttt                                                     <  2:207666/247‑1 (MQ=255)                                                                                                                                                                                               cgcgGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTTCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATc                                                 >  1:373166/1‑249 (MQ=255)                                                                                                                                                                                                  ggCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTa                                            <  1:364210/251‑1 (MQ=255)                                                                                                                                                                                                   gCATGGGATATATCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGTCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTac                                           <  1:194777/251‑1 (MQ=255)                                                                                                                                                                                                   gCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAAATGGCGCAGGCGTTTAATCAGcaac                                           >  1:109948/1‑248 (MQ=255)                                                                                                                                                                                                               gCCTTAGGTGCCTGTTATGCGTGTCGGCATGCGGTTGTTGCTGGGCTATTTTTGACTGGTGGCGGTGGCGGCCTGGGTTGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGATCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGTGCTGGCGCGTCCCGATTTGCTCTCTGGGTACCCACCGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGt                               <  2:9772/251‑1 (MQ=255)                                                                                                                                                                                                                     ggggCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCTGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTcgcg                         <  1:209500/251‑1 (MQ=255)                                                                                                                                                                                                                     ggggCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCCACGCAAGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTcgcc                         >  1:16459/1‑250 (MQ=255)                                                                                                                                                                                                                     ggggCCTGTAAGGCGTACCGGCATGCGGTTGCTGCTTGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCGGGCGGAGCTGGCGCGGCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTcgcg                         <  1:52470/251‑1 (MQ=255)                                                                                                                                                                                                                           tGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACt                                                                  <  1:223895/204‑1 (MQ=255)                                                                                                                                                                                                                              aaTGCGTATCGGCATGCGGTTGCTGCTGGGCTAGTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCg                <  1:323345/251‑1 (MQ=255)                                                                                                                                                                                                                                             catgttgctGCTTGGCTATTATTTACGGGAGGCGGTGGCGGCCTGGTTCGTACTGGCTATTGTTGTCAAAG‑AAGTTAAACCGGTCGTGCGAAGAGCAACGGAGGGGACGTTAAGCGATACCTCAACGTTGCTGGCGGAGCGGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTGCGCGCCAATATCGGTGGCATTAACaaa  <  2:282371/247‑1 (MQ=255)                                                                                                                                                                                                                                               ggTTGCTGCTGGGCTCTTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGGACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATAGGTGGCATTAACaaa  >  1:282371/1‑248 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    TCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAA  >  NC_000913/4636474‑4636959

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