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	97.4%	119.3 / ‑4.9	39	coding (1382/1701 nt)	narQ	sensor histidine kinase NarQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/24);  minor base . (1/0);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCT  >  NC_000913/2586866‑2587334                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                acttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑Tg                                                                                                                                                                                                                            <  1:316485/251‑1 (MQ=255)       gTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATgcg                                                                                                                                                                                                                      >  1:67259/1‑251 (MQ=255)        tATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATgcgc                                                                                                                                                                                                                     >  1:426981/1‑251 (MQ=255)             aGTTGGCGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACTAAGTATCTTGGCCGATTTTTCCCCTGCTTTGAATGATGCTTATCGGCAGTTACGCGATCTGTTGACGACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCTCATTGAGTGAAATGCTGGATACGTTTCAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAGCCCGGTTAC‑TGGATGCGTAAATg                                                                                                                                                                                                                <  2:462566/251‑1 (MQ=255)                 gACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAgg                                                                                                                                                                                                            <  1:474398/251‑1 (MQ=255)                      tACTGAAGCGTTCGATCCCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCGGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATTCGCAAATGCAGGTGCa                                                                                                                                                                                                        <  2:502197/250‑1 (MQ=255)                       aCTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTGCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACGTGCCGCCTGACGCTGTAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGACTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGATGCAt                                                                                                                                                                                                       <  2:82419/250‑1 (MQ=255)                          gAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATttg                                                                                                                                                                                                    <  1:157474/250‑1 (MQ=255)                            aGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATttgttg                                                                                                                                                                                                 <  1:396275/251‑1 (MQ=255)                                    aTACCGGAAGCTAACGACCCCGCGCAAAGGATCCTGTCCGATTTTTCCCAGGCGCTGAATGATGCTTATCGGCAGTTACGCGAACTCATGTCTACGTTCCGCCTGAAGCTGCAGCGGGAGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACACCATCACACCCGCGCCAAACTTACCCGCGACTGGCGTCTGCCCACCCCGG‑ACATGGATGCGAAAATGCAGGTGCCTTTGTTGCCAattgt                                                                                                                                                                                         >  2:69858/1‑249 (MQ=255)                                                  cGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGc                                                                                                                                                                           >  1:460935/1‑251 (MQ=255)                                                     cACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCCACCCTGGTAC‑TGGATGCGCAAATGCAGTTGCATTTGTTGCACATTATTCGCTAAGCGTTGCTg                                                                                                                                                                         >  1:133820/1‑250 (MQ=255)                                                             aaaGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCGCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATGTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATg                                                                                                                                                                <  2:184733/251‑1 (MQ=255)                                                                         ccGATTTTTCCCATGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGCGGGAAATTCTGGATACGTTACAAAATCAAACCAGCGCCACACTTACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATTCGCACATGCAGGTGCATTTTTTGCAAATTATTCGCGACGCGGTGCTGAATGCTGTGCAGCacgccacc                                                                                                                                                    >  2:84247/1‑249 (MQ=255)                                                                                 tCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGCTGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTTAATGCGATGAAGCACGCCAACGCCAGCg                                                                                                                                             >  1:197163/1‑250 (MQ=255)                                                                                  cccAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCCGCAGGCGGATCTCCCCTCCGCATTGGGGGAAATGCCGGGTACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGACTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGATAAATGCGATGAAGCacgccaacgcacgcga                                                                                                                                            >  2:231694/1‑244 (MQ=255)                                                                                                   gCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCg                                                                                                                           <  2:391826/250‑1 (MQ=255)                                                                                                       aTCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCCAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgccaacgcc                                                                                                                                                 >  1:501790/1‑224 (MQ=255)                                                                                                       aTCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCacgccaacgcc                                                                                                                                                 <  2:501790/224‑1 (MQ=255)                                                                                                          ggCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGTCAGCGAAATCGCCGTCAGTTGCGTCACcgcg                                                                                                                   <  2:498137/251‑1 (MQ=255)                                                                                                                     gAGCTGTTGAATACTTTCCGCCTGACGCGGCAGCCGGCGGATCTCCCCTCAGCATTGAGCGACATGCTGGATACGATACAAAATCAAACCAGCGCCAAACTGACCATCGACTGCAGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCa                                                                                                         >  2:342288/1‑250 (MQ=255)                                                                                                                            tGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCacac                                                                                                  <  1:57354/250‑1 (MQ=255)                                                                                                                             gACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACg                                                                                                 <  1:401764/250‑1 (MQ=255)                                                                                                                                 actTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTta                                                                                            <  2:176889/251‑1 (MQ=255)                                                                                                                                      ccGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGTGAAATGCTGGATACGTTACAAAAACAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGGTGCGTCGCCGCGCCGGACGTCAATCACACGGTTTATATc                                                                                        <  2:459385/250‑1 (MQ=255)                                                                                                                                        gCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGt                                                                                     <  1:10568/251‑1 (MQ=255)                                                                                                                                             aCGCTGCAGCATGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGAGACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAAGCTCCGTCAGTGGCGTCACCGCGCCGGACGGCAATCTCACGGTTTATATCCGTGATaa                                                                                <  2:67212/251‑1 (MQ=255)                                                                                                                                                 tGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAACCCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACggg                                                                            >  1:63050/1‑251 (MQ=255)                                                                                                                                                                      ccGCATTGAGGGAAATGCTGGATACGGTACAAAATCAAACCAGCGCCAAACGGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGACGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCATCGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGaaa                                                       <  2:224018/251‑1 (MQ=255)                                                                                                                                                                        gCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGa                                                     <  1:66403/251‑1 (MQ=255)                                                                                                                                                                                ggAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCTTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCCCGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATACCGGGATTGGTATCGGTGAACCGAAAGAACCCGa                                               >  1:268052/1‑249 (MQ=255)                                                                                                                                                                                      gCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCAtt                                       >  2:128216/1‑251 (MQ=255)                                                                                                                                                                                         ggATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTAt                                     <  1:401427/250‑1 (MQ=255)                                                                                                                                                                                               gTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAAc                                                   >  2:3880/1‑230 (MQ=255)                                                                                                                                                                                               gTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAAc                                                   <  1:3880/230‑1 (MQ=255)                                                                                                                                                                                                       aTCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATg                      <  2:67259/251‑1 (MQ=255)                                                                                                                                                                                                         cAAAACAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGTTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGTGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATgc                     <  2:365029/250‑1 (MQ=255)                                                                                                                                                                                                                   ccAAACTGACCCGCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCTGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGATAGACCCCGAAGGTCATTATGGTCTGAATATCATGCGCGAAcgcgcg          <  2:52107/251‑1 (MQ=255)                                                                                                                                                                                                                             ccTCGACTGCCGTCTGCCAACCCTGGTAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCt  >  1:281999/1‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                ACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCAC‑TGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGTGCTGAATGCGATGAAGCACGCCAACGCCAGCGAAATCGCCGTCAGTTGCGTCACCGCGCCGGACGGCAATCACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCT  >  NC_000913/2586866‑2587334

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