Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,587,112	C→T	A461V (GCA→GTA)	narQ →	sensor histidine kinase NarQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,587,112	0	C	T	94.6%	105.7 / NA	38	A461V (GCA→GTA)	narQ	sensor histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (15/20);  minor base G (1/1);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGA  >  NC_000913/2586879‑2587348                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             gTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCTTCTGCCAACCCTGGTACTGGATGCGCAAATGc                                                                                                                                                                                                                             <  1:182682/251‑1 (MQ=255) gTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGc                                                                                                                                                                                                                             <  1:132632/251‑1 (MQ=255)            ggAAGCGTGCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTGCCCAGGCGTTGAATGATGCTTAGCGGCAGTTACGCGAGCTGTTGACTACTGTCCGCCTGACGCTGCAGCAGTCGGATCTCCCCTCCGCATTGATGGAAATGCTGGAGACGTTACAAAATCAAACCAGCGCCAAACGTG‑CCCTCGACTGCCGTGTGCCAACCCTGTTACTGGATGCGCAAATGCTGGTGCATttg                                                                                                                                                                                                                  <  2:146066/250‑1 (MQ=255)              aaGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATttgtt                                                                                                                                                                                                                <  1:93649/251‑1 (MQ=255)                       aTACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTGACCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTTATGGAAATTCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAattat                                                                                                                                                                                                       <  1:371970/251‑1 (MQ=255)                           cGGAAGATAACGCCACCGCACAACGTATCATGACCGATTTTTCACAGGCGTTGAATGATGCTTATCGGCAGATACGCGAGCTGTTGCCTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCACTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGGACTGGATGCGAAAATGCAGGAGCATTTGATGCAAATTAATcgc                                                                                                                                                                                                   >  1:416908/1‑251 (MQ=255)                                 aTAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCCAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGAGCATTTGTTGCAAATTATTCGCGAAGCg                                                                                                                                                                                             >  1:176378/1‑251 (MQ=255)                                  tAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAATTACGCGAGCGGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTTCCCTCCGCCTTGAGGGGAATGGTGGGTACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCCAATGCAGGTGCATTTTGTGCAAatta                                                                                                                                                                                                        >  2:286773/1‑239 (MQ=255)                                  tAACGCCACCGCACAAAGTATCATGGCCGATTTTGCCCAGGCGTTGAATGATGCTTATCGGCAGTTCCGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGTCTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAACGCAGGTGCAATTGTTGTAAatta                                                                                                                                                                                                        <  1:286773/239‑1 (MQ=255)                                   aaCGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGt                                                                                                                                                                                           <  1:346739/251‑1 (MQ=255)                                     cGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCCTTTGTTGCAAATTATTCGCGAAGCGGTgg                                                                                                                                                                                         >  2:362157/1‑250 (MQ=255)                                      gCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCCGCGCCCAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCCGGGGCCTTTGTTGCAAATTATTCGCGAAGCGGTGCt                                                                                                                                                                                        >  2:321415/1‑251 (MQ=255)                                                   gTATCACGGCAGATATTTCACAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGAAGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑AGGCCCTCGACTGCCGTCTGCCAAACCTGGTACTGGATGCGCACATGCAGGGGCTTTTGTTGCCAATTATTAGAGAAGCGGTGGGGAATGCGATGa                                                                                                                                                                             >  2:452759/1‑249 (MQ=255)                                                     aTCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCCGGGGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGacgc                                                                                                                                                                          >  2:306759/1‑247 (MQ=255)                                                                     cccAGGAGTTGAATGATGCTTATCGGCAGTTACCCGAGCTGTTGAATACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCCACCCTGGTACTGGATGCGCAAATGCAGGGGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGa                                                                                                                                                          >  2:62233/1‑250 (MQ=255)                                                                        aGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TTACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCGGAATGCGATGAAGCACGCCAACGCCAGCGAAATc                                                                                                                                                      <  1:62233/251‑1 (MQ=255)                                                                             ttGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGACATGCTGGGTACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAAACCTGGTACTGGATGCGCAAATGCAGGTGCCTTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGt                                                                                                                                                 >  1:78838/1‑251 (MQ=255)                                                                              tGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGt                                                                                                                                                 <  1:449419/250‑1 (MQ=255)                                                                                   gatgCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGc                                                                                                                                                                          <  2:363510/220‑1 (MQ=255)                                                                                   gatgCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGc                                                                                                                                                                          >  1:363510/1‑220 (MQ=255)                                                                                                aGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCg                                                                                                                              <  1:301797/251‑1 (MQ=255)                                                                                                           cTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATc                                                                                                                    >  2:241541/1‑250 (MQ=255)                                                                                                               tGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCacac                                                                                                                <  2:426770/250‑1 (MQ=255)                                                                                                                    actTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAATTGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGtt                                                                                                            <  2:356357/249‑1 (MQ=255)                                                                                                                            ccTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCGGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTg                                                                                                  <  2:181352/251‑1 (MQ=255)                                                                                                                                             cGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATc                                                                                 >  1:111599/1‑251 (MQ=255)                                                                                                                                                         ccGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGaaa                                                                     <  2:176378/251‑1 (MQ=255)                                                                                                                                                               tGAGGGAAATGCTGGATACGTTACAAACTCAAACCAGCGCCAAAC‑TGACACTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCACATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAAAGCCAGCGAAAACGCCGTCAGTTGCGTCACCGCGCCGGACGGCCATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGACCCGAAAGAAccc                                                               >  2:399871/1‑251 (MQ=255)                                                                                                                                                                      aaTGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTc                                                        <  1:306759/251‑1 (MQ=255)                                                                                                                                                                                                  aGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAAcg                            <  1:241541/251‑1 (MQ=255)                                                                                                                                                                                                          aC‑TGACCCTCGACTGGCGTCTGCCAACCCTGGGACTGTATGCGCAAAGGCAGGTGCATTTGTTGCAAATTATTCGCTAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAAGCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTGTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATTCGCGAACGCGCGGAAc                    <  2:416908/251‑1 (MQ=255)                                                                                                                                                                                                              gACCCTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATAACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGc                 >  2:249052/1‑251 (MQ=255)                                                                                                                                                                                                                cccTCGACTGCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCt                >  1:92756/1‑250 (MQ=255)                                                                                                                                                                                                                         gCCGTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAGCGCGCGGAACGGCTAGGTGGGa        <  1:117119/249‑1 (MQ=255)                                                                                                                                                                                                                            gTCTGCCAACCCTGGTACTGGATGCGCAAATGCAGCCGCACTTCTTGCAAATTCTTCGCGAAGCGCTGCTGCATGCGATGCAGCCCGAAAACCTCAACGAAATAGACCCCAGTTGCAACCCCGCGCCGCACCGAACTAACACGGTTTATATCCGTCATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGACATCATGGTCAGAATATCATGCGCGAACGCGCGGAACGGCTTGGTGGGCCGCTg   >  1:262335/1‑251 (MQ=255)                                                                                                                                                                                                                                gCCAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGa  <  2:79927/248‑1 (MQ=255)                                                                                                                                                                                                                                  cAACCCTGGTACTGGATGCGCAAATGCAGTTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAAGGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCTAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAg              <  1:443912/234‑1 (MQ=255)                                                                                                                                                                                                                                  cAACCCTGGTACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGAGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGCAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGACCGGCTAg              >  2:443912/1‑234 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             GTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAAC‑TGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGA  >  NC_000913/2586879‑2587348

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