Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,678,019	A→G	6.0%	Q120Q (CAA→CAG)	PP_2348 →	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,678,019	0	A	G	6.0%	196.6 / 2.4	66	Q120Q (CAA→CAG)	PP_2348	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base G (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

CTATCAGCGGCCGTTTTCCGGTTACTACTTCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATGGGCCTGGTGGCGCTGATCCTGATTCT  >  NC_002947/2677872‑2678165                                                                                                                                                    |                                                                                                                                                   cTATCAGCGGCCGTTTTCCGGTTACTACTTCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAAc                                                                                                                                                   <  5:266282/149‑1 (MQ=255)          gCCGTTTTCCGGTTACTACTTCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGccc                                                                                                                                          >  7:248054/1‑149 (MQ=255)           ccGTTTTCCGGTTACTACTTCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCt                                                                                                                                         <  3:166767/149‑1 (MQ=255)                ttCCGGTTACTACTTCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTg                                                                                                                                    >  2:151748/1‑149 (MQ=255)                 tCCGGTGACTACTTCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCCGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTg                                                                                                                                    <  8:284538/148‑1 (MQ=255)                   cGGTTACTACTTCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGcc                                                                                                                                  >  8:244485/1‑148 (MQ=255)                     gTTACTACTTCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGAc                                                                                                                               >  3:360713/1‑149 (MQ=255)                         ctactTCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTAc                                                                                                                            <  7:244485/148‑1 (MQ=255)                            ctTCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGgc                                                                                                                        >  8:168060/1‑149 (MQ=255)                            ctTCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGgc                                                                                                                        <  2:184863/149‑1 (MQ=255)                             ttCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGgcg                                                                                                                       >  8:118144/1‑149 (MQ=255)                                  cATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTgg                                                                                                                  >  3:304455/1‑149 (MQ=255)                                       aTTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAg                                                                                                             >  3:238302/1‑149 (MQ=255)                                        tttCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCa                                                                                                              >  1:8083/1‑147 (MQ=255)                                        tttCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCa                                                                                                              <  2:8083/147‑1 (MQ=255)                                             aCCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACAtc                                                                                                       >  5:20823/1‑149 (MQ=255)                                                 gggTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGa                                                                                                   >  6:160639/1‑149 (MQ=255)                                                     aCCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAg                                                                                               >  4:69096/1‑149 (MQ=255)                                                               ctcGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCTCGGAGGGCCAGCAGTTGCTGGccttc                                                                                                                                        <  4:197198/98‑4 (MQ=38)                                                               ctcGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCTCGGAGGGCCAGCAGTTGCTGGccttc                                                                                                                                        >  3:197198/1‑95 (MQ=38)                                                                tcGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAg                                                                                     >  5:130141/1‑148 (MQ=255)                                                                tcGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAg                                                                                     >  7:26293/1‑148 (MQ=255)                                                                tcGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAg                                                                                     >  8:50862/1‑148 (MQ=255)                                                                  gTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGa                                                                                   <  4:238302/148‑1 (MQ=255)                                                                  gTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGa                                                                                   <  4:304455/148‑1 (MQ=255)                                                                  gTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGa                                                                                   <  7:168060/148‑1 (MQ=255)                                                                   ttCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGAc                                                                                  >  2:1530/1‑148 (MQ=255)                                                                     cGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTgg                                                                                                                  <  6:333997/114‑1 (MQ=255)                                                                         gtgGGACCTGGACATGCCCAAGCCAGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCAc                                                                           >  5:162003/1‑149 (MQ=255)                                                                          tgGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCAc                                                                           <  4:360713/148‑1 (MQ=255)                                                                             gACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGt                                                                       <  8:248054/149‑1 (MQ=255)                                                                               ccTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTgc                                                                     >  6:120157/1‑149 (MQ=255)                                                                               ccTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTgc                                                                     <  5:160639/149‑1 (MQ=255)                                                                                     cATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAg                                                               >  7:340866/1‑149 (MQ=255)                                                                                        gCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACAtctc                                                                                                     <  8:5178/108‑1 (MQ=255)                                                                                        gCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACAtctc                                                                                                     >  7:5178/1‑108 (MQ=255)                                                                                         cccAAGCGCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAgtg                                                                                             >  7:243617/1‑115 (MQ=255)                                                                                         cccAAGCCCGCAGCGCCCGGGATTTACGATTGCCATGAACTCGTCCCGTAAGTCCAGCAACTGCTTGCCCTGCTTGCCGACTACCTGCGCCTGGGCCAGGACATCTCGATCAgtg                                                                                             <  8:243617/115‑1 (MQ=255)                                                                                          ccAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCg                                                                         <  8:242566/134‑1 (MQ=255)                                                                                          ccAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCg                                                                         >  7:242566/1‑134 (MQ=255)                                                                                          ccAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTc                                                          <  8:26293/149‑1 (MQ=255)                                                                                                 cGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCt                                                   >  1:225833/1‑149 (MQ=255)                                                                                                     gcgcCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTgca                                               <  6:20823/149‑1 (MQ=255)                                                                                                      cgcCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTgcag                                              <  7:50862/149‑1 (MQ=255)                                                                                                       gcCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTgcagc                                             <  1:1530/149‑1 (MQ=255)                                                                                                        cccGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTgcag                                              >  4:71578/1‑147 (MQ=255)                                                                                                         ccGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTgcag                                              >  4:68449/1‑146 (MQ=255)                                                                                                                 tCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGc                                     <  6:130141/147‑1 (MQ=255)                                                                                                                   cGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTgg                                 <  2:225833/149‑1 (MQ=255)                                                                                                                     aTAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTggg                                >  3:102615/1‑148 (MQ=255)                                                                                                                        gCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATg                            <  5:120157/149‑1 (MQ=255)                                                                                                                         ccATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATgg                           <  1:151748/149‑1 (MQ=255)                                                                                                                          cATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATggg                          >  4:71457/1‑149 (MQ=255)                                                                                                                           aTGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCtt                                                  >  1:147830/1‑124 (MQ=255)                                                                                                                           aTGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATGGGc                         >  8:22088/1‑149 (MQ=255)                                                                                                                           aTGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATGGGc                         >  1:256086/1‑149 (MQ=255)                                                                                                                               aCTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATGGGCCtgg                     <  2:256086/149‑1 (MQ=255)                                                                                                                                cTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATGGGCCtggt                    <  8:70074/149‑1 (MQ=255)                                                                                                                                   ggCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATGGGCCtggtg                   <  3:69096/147‑1 (MQ=255)                                                                                                                                        gcagggCCAGCAGTTGCTGGCCTTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTgg                                 >  7:341830/5‑128 (MQ=255)                                                                                                                                        gcagggCCAGCAGTTGCTGGCCTTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTgg                                 <  8:341830/124‑1 (MQ=255)                                                                                                                                              ccAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATGGGCCTGGTGGCGCTGATCCTGa      <  3:71578/149‑1 (MQ=255)                                                                                                                                               cAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATGGGCCTGGTGGCGCTGATCCTGAt     >  5:85821/1‑149 (MQ=255)                                                                                                                                                aGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATGGGCCTGGTGGCGCTGATCCTGAtt    <  3:68449/149‑1 (MQ=255)                                                                                                                                                 gCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATGGGCCTGGTGGCGCTGATCCTGATTc   >  5:199291/1‑149 (MQ=255)                                                                                                                                                  cAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATGGGCCTGGTGGCGCTGATCCTGATTCt  <  6:162003/149‑1 (MQ=255)                                                                                                                                                   aaCTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCtt                                                  >  4:19939/1‑100 (MQ=255)                                                                                                                                                   aaCTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCtt                                                  <  3:19939/100‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CTATCAGCGGCCGTTTTCCGGTTACTACTTCCGCATCGATTTCGACCAGGGTACCTGGAGGTCTCGTTCGCTGTGGGACCTGGACATGCCCAAGCCCGCAGCGCCCGGGCTTTCCGATAGCCATGAACTCGGCCCGGAAGGCCAGCAACTGCTGGCCCTGCGTGCCGACTACCGGCGCCTGGGCCAGGACATCTCGATCAGTGTTGCTCAGGACTATTCACCGGTGCGCGAAGGTTTCCGACGCTTGCAGCAGATCGGCCTGGGCATGGGCCTGGTGGCGCTGATCCTGATTCT  >  NC_002947/2677872‑2678165

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: ‑