Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,027,751	T→A	8.1%	I227F (ATC→TTC)	PP_3552 ←	two‑component system sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,027,751	0	T	A	8.1%	179.4 / 5.6	65	I227F (ATC→TTC)	PP_3552	two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (27/31);  new base A (2/3);  total (30/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 = 9.06e-01

CAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCC  >  NC_002947/4027608‑4027896                                                                                                                                                |                                                                                                                                                      cAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATgcgc                                                                                                                                                   >  5:364182/1‑148 (MQ=255)  aGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCATCGAAGCGCag                                                                                                                                                 <  7:350392/149‑2 (MQ=255)   ggCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAg                                                                                                                                                <  6:364182/149‑1 (MQ=255)    gCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAgg                                                                                                                                               <  6:180767/149‑1 (MQ=255)       cTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGa                                                                                                                                            <  5:329684/149‑1 (MQ=255)        tGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGAt                                                                                                                                           >  2:53036/1‑149 (MQ=255)             gTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTcg‑‑‑‑c                                                                                                                                  >  2:300753/1‑149 (MQ=255)             gTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGACGCGCAAGGTGATGTcg‑‑‑‑c                                                                                                                                  >  2:185181/1‑149 (MQ=255)              tAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCATCGAAGCGCa                                                                                                                                                  >  2:34162/1‑136 (MQ=255)              tAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGAGCGCCAGGCTTTCCTGCGCATCGAAGCGCa                                                                                                                                                  <  1:34162/136‑1 (MQ=255)                aCGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑Cgg                                                                                                                                <  6:104485/148‑1 (MQ=255)                      aCTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCACGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACg                                                                                                                          >  6:356079/1‑148 (MQ=255)                            gTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGacca                                                                                                                   >  6:365070/1‑149 (MQ=255)                                gCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTc                                                                                                               >  1:116503/1‑149 (MQ=255)                                  aTAGTTGAGCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCAGAGGCCAGCTCGCCGAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGGGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACGACCTcca                                                                                                             <  4:257935/149‑1 (MQ=255)                                  aTAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTcca                                                                                                             >  7:260878/1‑149 (MQ=255)                                   tAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTccac                                                                                                            <  5:365070/149‑1 (MQ=255)                                    aGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTccac                                                                                                            >  7:35353/1‑148 (MQ=255)                                           caccGCGGCCAGTGGCTGATTGTTCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACcaccgcac                                                                                                    >  1:251689/1‑149 (MQ=255)                                           caccGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTc                                                                                                               >  6:243221/1‑138 (MQ=255)                                           caccGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTc                                                                                                               <  5:243221/138‑1 (MQ=255)                                                 ggCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGgtgt                                                                                              <  3:299053/149‑1 (MQ=255)                                                   ccAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAg                                                                                            <  2:220734/149‑1 (MQ=255)                                                    cAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAgg                                                                                           >  8:59874/1‑149 (MQ=255)                                                           tgattgatCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAgg                                                                                                                                               <  1:254904/94‑1 (MQ=255)                                                           tgattgatCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAgg                                                                                                                                               >  2:254904/1‑94 (MQ=255)                                                                gatCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACAACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCg                                                                               >  8:100419/1‑149 (MQ=255)                                                                  tctcGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACAAGCTGACAACCTCCACCACCGCACCGGTCTAGGTATCACGGATCGcc                                                                             >  5:48632/1‑149 (MQ=255)                                                                         ggCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTg                                                                       <  8:35353/148‑1 (MQ=255)                                                                          gCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGc                                                                     <  5:356079/149‑1 (MQ=255)                                                                          gCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGc                                                                     <  1:76732/149‑1 (MQ=255)                                                                               gCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCTCCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGAttt                                                                >  4:243733/1‑149 (MQ=255)                                                                                cccGAGGCGAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTc                                                               <  4:144012/149‑1 (MQ=255)                                                                                cccGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTc                                                               <  2:116503/149‑1 (MQ=255)                                                                                 ccGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTc                                                               >  2:250063/1‑148 (MQ=255)                                                                                   gAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCTCCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGaa                                                            <  3:243733/149‑1 (MQ=255)                                                                                   gAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGaa                                                            <  1:53036/149‑1 (MQ=255)                                                                                      gCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAAcc                                                          >  3:239824/1‑148 (MQ=255)                                                                                       ccAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCCTCCTGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGaa                                                            <  1:300753/149‑1 (MQ=255)                                                                                         aGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACg                                                      <  2:251689/149‑1 (MQ=255)                                                                                          gCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACgg                                                     <  1:250063/149‑1 (MQ=255)                                                                                          gCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACgg                                                     <  1:136178/149‑1 (MQ=255)                                                                                            tCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACgg                                                     >  1:41734/1‑147 (MQ=255)                                                                                             cgcgCAGGGTAGAGAGTCGGGCTCTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCAGCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTa                                                   <  6:48632/148‑1 (MQ=255)                                                                                              gCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGACGATTTCGAACCAACGGTa                                                   <  5:50784/147‑1 (MQ=255)                                                                                                    ggTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGcctcct                                            <  8:260878/148‑1 (MQ=255)                                                                                                        gCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGACGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGc                                       <  1:185181/149‑1 (MQ=255)                                                                                                           aGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTg                                    <  3:166311/149‑1 (MQ=255)                                                                                                            ggCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTgg                                   >  5:330854/1‑149 (MQ=255)                                                                                                                        gcgcCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCa                        <  4:239824/148‑1 (MQ=255)                                                                                                                          gcCCGGCTTTCCTGCGCTTCGATGCGCCCGGTGATGTCG‑‑‑‑CGGGACACGCTGCCCACCTCCACCACCGCACCGGTGTCGATCTCACTGATCGCCCTAAAATAAAAAAAAAAAAAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAgg                      >  6:143968/1‑148 (MQ=255)                                                                                                                          gcCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGc                                                                     <  4:253941/101‑1 (MQ=255)                                                                                                                            cAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGGAGCCTCCTGTCGCCTGGCGTACCCGGCAGGGCa                   >  8:333026/1‑149 (MQ=255)                                                                                                                            cAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGAAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCa                   >  5:270178/1‑149 (MQ=255)                                                                                                                              ggCTTTCCTGCGCTTCGATGCGCAAGGTGATGGCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATc                 <  8:252553/149‑1 (MQ=255)                                                                                                                               gCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGc                                  >  2:138555/1‑131 (MQ=255)                                                                                                                               gCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGc                                  <  1:138555/131‑1 (MQ=255)                                                                                                                                   ccttGCGCATCGAAGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGtgtg            >  4:17920/4‑149 (MQ=255)                                                                                                                                     cTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATc                 <  6:45548/142‑1 (MQ=255)                                                                                                                                     cTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATc                 >  5:45548/1‑142 (MQ=255)                                                                                                                                        cgcATCGAAGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAg        <  3:17920/148‑1 (MQ=255)                                                                                                                                          cTTCGATGCGCAAGGAGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACgg                                                     >  5:232904/1‑101 (MQ=255)                                                                                                                                          cTTCGATGCGCAAGGAGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACgg                                                     <  6:232904/101‑1 (MQ=255)                                                                                                                                           ttCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGt    >  4:84695/1‑149 (MQ=255)                                                                                                                                             cGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTcc  >  6:103272/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                      CAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCG‑‑‑‑CGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGTAGCCGTCC  >  NC_002947/4027608‑4027896

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