Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,871,607	G→C	7.8%	H258Q (CAC→CAG)	PP_3420 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,871,607	0	G	C	7.8%	186.5 / 4.7	61	H258Q (CAC→CAG)	PP_3420	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (24/32);  new base C (2/3);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.73e-01

CGCAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTCCACCAT  >  NC_002947/3871463‑3871753                                                                                                                                                  |                                                                                                                                                   cgcAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGt                                                                                                                                                   <  4:103139/146‑1 (MQ=255)    agcagcGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTCCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTgcagc                                                                                                                                              >  3:221841/1‑148 (MQ=255)    agcagcGTTGGCGGTTTCCCGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCACCTgcagca                                                                                                                                             <  8:257143/149‑1 (MQ=255)      cagcGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCata                                                                                                                                           >  3:373262/1‑149 (MQ=255)      cagcGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCata                                                                                                                                           >  8:99149/1‑149 (MQ=255)       agcGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCatat                                                                                                                                          <  1:190795/149‑1 (MQ=255)       agcGTTGGCGGTTACCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGCCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCACCTgcagcagga                                                                                                                                          <  5:102447/149‑4 (MQ=255)          gTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCACCTgcagcaggacag                                                                                                                                       <  8:10308/149‑7 (MQ=255)           ttGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTgc                                                                                                                                      >  3:228636/1‑149 (MQ=255)                    tCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGgcg                                                                                                                             >  1:302011/1‑149 (MQ=255)                         ggCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGa                                                                                                                        >  3:87260/1‑149 (MQ=255)                         ggCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGa                                                                                                                        <  6:39280/149‑1 (MQ=255)                           cAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATc                                                                                                                      <  5:303471/149‑1 (MQ=255)                                  ggCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAAc                                                                                                                <  5:355305/148‑1 (MQ=255)                                  ggCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTCCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACa                                                                                                               <  4:221841/149‑1 (MQ=255)                                     gcTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGc                                                                                                            <  7:99149/149‑1 (MQ=255)                                           gCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGc                                                                                                            <  1:340567/143‑1 (MQ=255)                                           gCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGc                                                                                                            >  2:340567/1‑143 (MQ=255)                                                tGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGACTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGAt                                                                                                  <  5:12939/148‑1 (MQ=255)                                                 gTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGAGTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGAt                                                                                                  >  1:314110/1‑148 (MQ=255)                                                  ttGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCaa                                                                                               <  6:297096/149‑1 (MQ=255)                                                     gCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCGAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGc                                                                                            <  4:228636/149‑1 (MQ=255)                                                       cAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGccc                                                                                          >  2:72358/1‑149 (MQ=255)                                                        aGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGccc                                                                                          >  7:141328/1‑148 (MQ=255)                                                           gggTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCcccacc                                                                                      <  2:302011/149‑1 (MQ=255)                                                                aCGCATGCGGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCcacccacc                                                                                  <  7:176583/148‑1 (MQ=255)                                                                 cGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACccacca                                                                                 <  2:145626/148‑1 (MQ=255)                                                                     tGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTc                                                                            <  1:207582/149‑1 (MQ=255)                                                                       cTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTc                                                                            >  4:364467/1‑147 (MQ=255)                                                                        tgtgAGCCACTGCGGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAg                                                                         >  2:395041/1‑149 (MQ=255)                                                                        gttcaGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCACCGCCCCCCCCCCCAGTGTCGAg                                                                         >  6:202072/5‑149 (MQ=255)                                                                            aGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGcc                                                                      <  1:72358/148‑1 (MQ=255)                                                                             gCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTTCAGCTGCAGCATATGTGCCGCCCGGCGCACCATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCa                                                                    >  7:279449/1‑149 (MQ=255)                                                                                 cTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATc                                                                <  4:87260/149‑1 (MQ=255)                                                                                       acgacATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGcaca                                                          >  2:48744/1‑149 (MQ=255)                                                                                       acgacATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGcaca                                                          >  2:236330/1‑149 (MQ=255)                                                                                            aTTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGcc                                                     >  3:299093/1‑149 (MQ=255)                                                                                             tttCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCcg                                                    >  1:370567/1‑149 (MQ=255)                                                                                              ttCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCcgc                                                   >  6:150883/1‑149 (MQ=255)                                                                                                   ccAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGAGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTACATCAACGCCCCACCCACCAGTGTCGCGGCCCAGATCAGCACATAGCCGCGATg                                               >  1:83057/1‑148 (MQ=255)                                                                                                         cAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATgc                                          >  1:374155/1‑147 (MQ=255)                                                                                                         cAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATgc                                          <  2:374155/147‑1 (MQ=255)                                                                                                          aaCGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAg                                                                         <  3:201335/115‑1 (MQ=255)                                                                                                          aaCGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGcc                                       <  2:370567/149‑1 (MQ=255)                                                                                                          aaCGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGcc                                       <  3:70032/149‑1 (MQ=255)                                                                                                          aaCGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTCCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAg                                                                         >  4:201335/1‑115 (MQ=255)                                                                                                              tAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGAGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCTAGATCAGCACATAGCCGCGCTGGATGCGCCGGAt                                   <  1:252792/149‑1 (MQ=255)                                                                                                                gg‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGa                                 <  4:299093/149‑1 (MQ=255)                                                                                                                 g‑TTTCGCTTGCCAGCTGGCTGTCCTGCTGCAGGTGCAGGATTTGTGCCGCGCGGCGCACGATCCAGAACAGGCCGAATTAGATCAACGCACCACCCACCAGTGTCGATGGCCAGATGAGCTCATAGCGGCGCTGGATGCGCCGGATGAg                                <  8:279449/149‑1 (MQ=255)                                                                                                                 g‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAg                                <  1:49782/149‑1 (MQ=255)                                                                                                                 g‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAg                                >  2:385798/1‑149 (MQ=255)                                                                                                                           gCCACCAGCCTGTCCTGCTGCATGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGGACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGc                         <  1:48744/147‑1 (MQ=255)                                                                                                                            ccaccaGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGc                      <  2:314110/149‑1 (MQ=255)                                                                                                                                   ccTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCTTTg                >  1:263223/1‑148 (MQ=255)                                                                                                                                      catatgctgcACCTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTg                >  8:95806/5‑145 (MQ=255)                                                                                                                                            ctgcAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCgg                       <  7:350362/132‑1 (MQ=255)                                                                                                                                            ctgcAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCgg                       >  8:350362/1‑132 (MQ=255)                                                                                                                                            ctgcAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTcca      <  3:364467/149‑1 (MQ=255)                                                                                                                                            ctgcAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTc        <  1:385798/147‑1 (MQ=255)                                                                                                                                               cAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTccac     <  4:373262/147‑1 (MQ=255)                                                                                                                                                aGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTCCACCAt  <  8:141328/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CGCAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGG‑TTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTCCACCAT  >  NC_002947/3871463‑3871753

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