Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,827	C→T	D44D (GAC→GAT)	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,827	0	C	T	93.3%	113.8 / ‑1.6	45	D44D (GAC→GAT)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (24/18);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

CCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCT  >  NC_000913/4636596‑4637069                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   cccACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACGGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAGGAAGTTAAACCAGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCgg                                                                                                                                                                                                                                 >  1:47234/1‑251 (MQ=255)                     gCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGAGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGCCTCGGCCTGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGAGGTTAATCGATACCGCCACGTTTCGGGGAGAGCTGGGTCGTCCCGCTTTg                                                                                                                                                                                                             >  1:279435/1‑251 (MQ=255)                      ccAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGc                                                                                                                                                                                                            >  1:233709/1‑251 (MQ=255)                                cgcCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCCTGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACAGGGCGTGCGAAGAGCAACGGAGGGGACGTTAAACGATACCGCAACGTTGCTGGCGGAGCTGGCGCATCCCGCTTTGCTCTCTtgggg                                                                                                                                                                                                  >  1:397021/1‑250 (MQ=255)                                       aaCCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACg                                                                                                                                                                                           <  1:256099/251‑1 (MQ=255)                                        aCCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGCTACCGCAACGTTGCTGGCGGAGATGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGc                                                                                                                                                                                          >  1:446109/1‑251 (MQ=255)                                             gATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGGACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGGCGTTAATCGCTAACGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCCTggg                                                                                                                                                                                     >  1:64133/1‑251 (MQ=255)                                                                  cATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAAt                                                                                                                                                                <  1:596414/251‑1 (MQ=255)                                                                        ggCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGCTTTGCTCTCTGGGGGCCCAACGCATGGGCAACTGGCGCAGGGGTTTAATCAGCt                                                                                                                                                           >  1:529616/1‑250 (MQ=255)                                                                                   taGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGGGGTGGCGCGTCCCGATTTGCTCTCTGGGGGGCCAACGCATGGGCAACTGGCGCCGGGGGTTAATCAGCTACAACATCGccc                                                                                                                                               >  1:97660/1‑251 (MQ=255)                                                                                                     aTGCGTATCGGCGTGCGGTTGCTGCTGTGCTTGTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGTGTCGTTAATCGATACCGCTACGTTGCTGGCGGAGCTTGCGCGTCCCGATTTGCTCTCTGGGTTCCCAACGCAGGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCgg                                                                                                                             <  2:469456/251‑1 (MQ=255)                                                                                                         gTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGc                                                                                                                         <  1:36553/251‑1 (MQ=255)                                                                                                                cATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACa                                                                                                                  <  1:675297/251‑1 (MQ=255)                                                                                                                      gttgctGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACCCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACCTCGCCCGTTTCGCGCCAATATCGGGGGCATTTACAAAATgc                                                                                                            >  1:291978/1‑251 (MQ=255)                                                                                                                       ttgctgCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGACGTTAAACCGGGCGTGCGAAGAGCAAAGGAGGGGACGTTAAACGATACCGCAACGTTGCTGGCGGAGCTGGCGCGGCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCCGATCAAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACCAAGTgcg                                                                                                           >  1:602867/5‑251 (MQ=255)                                                                                                                       ttgctgCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCCATATCGGGTGCATTAACAAAGTgcg                                                                                                           >  1:596264/5‑251 (MQ=255)                                                                                                                           tgCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAAAAAAGTGCGCAAt                                                                                                       >  1:656271/1‑251 (MQ=255)                                                                                                                            gCTGGGTTATTTTTTACGGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTGCGCGCCAATATCGGTTGCATTAACAAAGTGCGCAATg                                                                                                      <  2:446109/251‑1 (MQ=255)                                                                                                                               gggCTTTTTTTTACGTGTGGCGGGGGCGGCCTGGTTCGTACTGGCAATTTTTGGCAAAGAAGTTAATCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGTGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGTCGTAGGCGTTTAATCAGCTACAACATCGCCCGTTGCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGaa                                                                                                    <  2:379891/250‑1 (MQ=255)                                                                                                                               gggCTATTTTTTACTGGGGACGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGGCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATTCCGCAACGTTGCTGGCTGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGaa                                                                                                    <  1:60120/250‑1 (MQ=255)                                                                                                                               gggCGATTTGTTACTGGTGGCGGTGGGAGGCTGGTTCGTACTGGCCATTTTTGTCAAGGAAGGTAATCCGGGCGTGCGAAGAGCAACGGAGGGGACGGTGATCGACACCGCGACGTTGCTGGCGGAGCTGGCGCGTACCGAGTTGCTCTCTGGGGACCCAACGCATGGGCAACGGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTGCGCGCCAATATCGGTGGGATTAACAAAGTGCGCAATGaag                                                                                                   <  2:47234/251‑2 (MQ=255)                                                                                                                                ggCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGGACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAata                                                                                                  <  1:482146/251‑1 (MQ=255)                                                                                                                                    aTTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCCGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCCATATCGGTGGCATTAACAAAGTGCGCAATGAATATCAt                                                                                              >  1:255465/1‑251 (MQ=255)                                                                                                                                      tttttACTGGTGGCGGTGGCGGCCTGGTTCGTACTGTCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTa                                                                                                                                                                  <  2:379879/181‑1 (MQ=255)                                                                                                                                      tttttACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTa                                                                                                                                                                  >  1:379879/1‑181 (MQ=255)                                                                                                                                              gtgggcggtggcGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTTGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTatat                                                                                      <  1:275639/246‑1 (MQ=255)                                                                                                                                                 ggcggtggcGGCCTGGATCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTACCGATTTGCTCTCTGGGGACCCAACGCAGTGGCAAACGGCGCCGGCGTTTAATCAGCTACAAACTTGCCCCTTTCGCGCCAAAATCGGTGGCATTAACAAAAGGGGCCATGACTATCATGTCTTAATGACCg                                                                                 >  2:374813/1‑251 (MQ=255)                                                                                                                                                  gcggtggcGGCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAATTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGa                                                                                <  1:492521/251‑1 (MQ=255)                                                                                                                                                     gtggcgtccTGGTTGGTACTGGGTATTTGTGTGATAGAAGTTAAACCGGTGGTTCGAAGAGCAACGGAGTGGCCTTTAAGGGATACCGCAACGTTTCTGGCGTAGCTTGCGCGGCCCGATTTGCGCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCGCGTTTCGCGCCAAGATCGTTGGTATTAACAAAGTGCGCAAGGAATGTCATGTCTATATGCCCggttc                                                                             <  2:111256/244‑5 (MQ=255)                                                                                                                                                          ggCCTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACGGGGAGTGCGAAGAGAAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGGTTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAACGGCGCAATGAATATCATGTCTATATGACCGATTCGCAgg                                                                        >  1:565799/2‑251 (MQ=255)                                                                                                                                                            ccTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAgggt                                                                      <  1:20422/251‑2 (MQ=255)                                                                                                                                                             cTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAgggt                                                                      >  1:183773/1‑249 (MQ=255)                                                                                                                                                             cTGGTTCGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGGTGCTGGCGGAGCTGGCGCGGCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGGTTAATCAGCTACAACATCGACCGTTTCCCGCCAATAACGGTGGCATTAACAACATCCGCAATGAATATCATgacaataaaaccgattcgccgggt                                                                      >  2:134192/1‑227 (MQ=255)                                                                                                                                                                     tACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCCACTGGCGCCGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCCATGAAAATCATGTCTCTATGACCGATTCGCAGGGTAAAGTATTg                                                             >  2:229325/1‑251 (MQ=255)                                                                                                                                                                             aTTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGAtt                                                      <  1:374813/250‑1 (MQ=255)                                                                                                                                                                             aTTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGAAGTTAATCGATACCGCAACGTTGCTGGCGGCGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGAAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGGGGCATTAACAAAGTGCGGAATGAATATCAAGTCTATATGACCGATTCGCCGGGTAAAGAATTGTTCGATTc                                                     >  1:565927/1‑251 (MQ=255)                                                                                                                                                                              gtgtgttCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAATGGAGGGGACGGTAATCGATACCTGAACGTTGCTTGCGGAGCTGGCGCGGCCCGATTTTCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCATCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGGGCGCAATGAATATCATGTCTATATGACCGTTTCGCAGGGTAAAGTATTGGTCGATTc                                                     <  2:333231/244‑1 (MQ=255)                                                                                                                                                                                     cAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAAt                                              >  1:489152/1‑250 (MQ=255)                                                                                                                                                                                      aaaGAAGTTAAACCGGGCGGGCGAAGATCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATaa                                            <  2:480453/251‑1 (MQ=255)                                                                                                                                                                                       aagaagTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATaa                                            >  1:126663/1‑250 (MQ=255)                                                                                                                                                                                              tAAACCGGGCGTGCGAAGAGCATCGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCgg                                                                                                                             <  2:266902/162‑1 (MQ=255)                                                                                                                                                                                              tAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCgg                                                                                                                             >  1:266902/1‑162 (MQ=255)                                                                                                                                                                                                              agagCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTAGCAGGGTAAAGTATTGTTCGATTAGTCAAATAAAGCCGTTGGGCAGGATTATTggcg                    >  2:572861/1‑251 (MQ=255)                                                                                                                                                                                                                               gttaaTCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATg                                                                                                      <  1:390626/148‑1 (MQ=255)                                                                                                                                                                                                                                 taaTCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGGTTTGCTCTCTGGGGACCCAACGCATGGTCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGGTTCGCGCCAATATCGGTTGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTTTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCACGCTGGAATGACGTCTGGcg  >  1:519507/3‑249 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   CCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCT  >  NC_000913/4636596‑4637069

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