Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,622,284	G→A	7.0%	E400K (GAA→AAA)	PP_1421 →	sensor histidine kinase TctE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,622,284	0	G	A	7.0%	165.0 / 2.3	57	E400K (GAA→AAA)	PP_1421	sensor histidine kinase TctE
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base A (2/2);  total (30/27)
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

CGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTG  >  NC_002947/1622135‑1622424                                                                                                                                                      |                                                                                                                                                 cgctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACC‑GGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtgg                                                                                                                                                  <  1:98/149‑1 (MQ=255)  gctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtgg                                                                                                                                                  <  6:185629/149‑1 (MQ=255)            gAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACTACgg                                                                                                                                        >  7:256551/1‑149 (MQ=255)                 gcCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAg                                                                                                                                   >  7:245202/1‑149 (MQ=255)                 gcCGGTGTGGGGTTAGGGTGCGCCGACATTGATCAACGAGTTGTGGAGCAACCTGGTGGACACTGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGAGTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAg                                                                                                                                   <  8:256551/149‑1 (MQ=255)                   cGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtgg                                                                                                                                                  >  3:220887/1‑132 (MQ=255)                   cGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtgg                                                                                                                                                  <  4:220887/132‑1 (MQ=255)                        tgGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTcg                                                                                                                            >  2:295241/1‑148 (MQ=255)                          gTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTccc                                                                                                                           <  3:63484/148‑1 (MQ=255)                          gTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCg                                                                                                                          >  6:234471/1‑149 (MQ=255)                              aaGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGc                                                                                                                      >  3:252778/1‑149 (MQ=255)                                gggggAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCg                                                                                                                    >  4:103581/1‑149 (MQ=255)                                   ggAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgc                                                                                                                 <  3:149052/149‑1 (MQ=255)                                   ggAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgc                                                                                                                 <  6:136096/149‑1 (MQ=255)                                   ggAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgc                                                                                                                 <  3:144454/149‑1 (MQ=255)                                           cATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACgg                                                                                                         >  1:142458/1‑149 (MQ=255)                                                 tCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc                                                                                                   >  4:123080/1‑149 (MQ=255)                                                 tCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGAAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc                                                                                                   >  5:38392/1‑149 (MQ=255)                                                 tCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGAAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc                                                                                                   >  5:42929/1‑149 (MQ=255)                                                   aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGa                                                                                                 <  4:252778/149‑1 (MQ=255)                                                   aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCAGGCAGTGCTGGAAGTGGAAgacgt                                                                                                                                           >  8:115364/1‑106 (MQ=255)                                                      gAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCg                                                                                              >  1:77656/1‑149 (MQ=255)                                                        ctcgttgAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGt                                                                                                                         >  7:167856/4‑120 (MQ=255)                                                        ctcgttgAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGt                                                                                                                         <  8:167856/117‑1 (MQ=255)                                                          tgtttcgCAACCTGGTTGAGAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGTGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGtt                                                                                           <  6:212919/142‑1 (MQ=255)                                                          tgttgAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGtt                                                                                           <  2:131961/148‑1 (MQ=255)                                                           gttgAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAAGTATTCCCGTGGCCGAGCGCGAACGGgtgt                                                                                                     >  3:3494/1‑137 (MQ=255)                                                           gttgAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAAGTATTCCCGTGGCCGAGCGCGAACGGgtgt                                                                                                     <  4:3494/137‑1 (MQ=255)                                                              gAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCGCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTAc                                                                                       <  1:295241/148‑1 (MQ=255)                                                               aGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgc                                                                                                                 <  7:288656/121‑1 (MQ=255)                                                               aGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTAcc                                                                                      <  3:103581/148‑1 (MQ=255)                                                               aGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGCCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgc                                                                                                                 >  8:288656/1‑121 (MQ=255)                                                                 cAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCATGGCCGAGCGCGAACGGgt                                                                                                       >  1:122913/1‑129 (MQ=255)                                                                    ccTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCCGCGaa                                                                                <  8:245202/149‑1 (MQ=255)                                                                      tggtggACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGt                                                                              >  7:255350/1‑149 (MQ=255)                                                                         tggACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCg                                                                                  <  4:86311/142‑1 (MQ=255)                                                                         tggACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCg                                                                                  >  3:86311/1‑142 (MQ=255)                                                                           gACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCa                                                                         <  8:255350/149‑1 (MQ=255)                                                                            aCAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAg                                                                        <  2:142458/149‑1 (MQ=255)                                                                               aTGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAggg                                                                      >  4:116593/1‑148 (MQ=255)                                                                                  cgctGGCGTATGCGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGt                                                                  >  5:295051/4‑149 (MQ=255)                                                                                      ggCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGAAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTgg                                                                <  6:42929/147‑1 (MQ=255)                                                                                      ggCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGAAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTgg                                                                <  6:38392/147‑1 (MQ=255)                                                                                        cgcATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGt                                                              >  5:287827/1‑147 (MQ=255)                                                                                         gcATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAg                                                           >  2:189040/1‑149 (MQ=255)                                                                                                   gCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAt                                                 <  2:77656/149‑1 (MQ=255)                                                                                                           cAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCg                                     <  6:287827/149‑1 (MQ=255)                                                                                                            aaTGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGc                                      >  8:53184/1‑147 (MQ=255)                                                                                                            aaTGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGc                                      >  5:24881/1‑147 (MQ=255)                                                                                                             aTGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGc                                      >  6:272768/1‑146 (MQ=255)                                                                                                             aTGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCgag                                   <  3:123080/149‑1 (MQ=255)                                                                                                               gtgATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGc                                      >  4:234497/1‑144 (MQ=255)                                                                                                               gtgATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGc                                      <  3:234497/144‑1 (MQ=255)                                                                                                                tgATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCGAGATc                                <  5:234471/149‑1 (MQ=255)                                                                                                                 gATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCGAGATCt                               >  7:78992/1‑149 (MQ=255)                                                                                                                        cGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCGAGATCTGTCggg                         <  6:295051/148‑1 (MQ=255)                                                                                                                         gggTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCGAGATCTGTCGGgcg                       >  3:272174/1‑149 (MQ=255)                                                                                                                           gTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCGAGATCTGTCGGGCGCa                     >  5:274626/1‑149 (MQ=255)                                                                                                                                      ggCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAg          <  3:116593/149‑1 (MQ=255)                                                                                                                                       gCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGt         >  1:65933/1‑149 (MQ=255)                                                                                                                                       gCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGt         <  1:189040/149‑1 (MQ=255)                                                                                                                                               ggaagtggaagACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTgg                                       <  1:295842/111‑1 (MQ=255)                                                                                                                                               ggaagtggaagACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTgg                                       >  2:295842/1‑111 (MQ=255)                                                                                                                                               ggaagtggaagACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAAAGGGTGTTCGAGCGGTTATACCGGAGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTg  >  7:193255/1‑148 (MQ=255)                                                                                                                                                gaagtggaagACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGCTTATAAGACGTTGGCGAGATCTGTCGGGCGCACCTGGCGAAAGtcagcg    >  1:257002/1‑144 (MQ=255)                                                                                                                                                      |                                                                                                                                                 CGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCGGTTCTACCGGCGAAGTGCGCAGGGCAGTGGGTTAGGGCTGGCGAT‑‑‑‑CGTTGGCGAGATCTGTCGGGCGCACCTGGCGCAAGTGACCCTG  >  NC_002947/1622135‑1622424

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