Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,912,847	C→G	6.4%	G411A (GGC→GCC)	PP_3453 ←	sensor protein RstB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,912,847	0	C	G	6.4%	202.1 / 3.9	62	G411A (GGC→GCC)	PP_3453	sensor protein RstB
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base G (2/2);  total (31/31)
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.26e-01

ACCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCATC  >  NC_002947/3912703‑3912989                                                                                                                                                 |                                                                                                                                                aCCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCAcc                                                                                                                                             <  2:116547/149‑1 (MQ=255)                             ggCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGCTCGCCCGCCGGCTGATg                                                                                                                 >  5:27238/1‑149 (MQ=255)                                            gCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGcc                                                                                                   >  3:255035/1‑148 (MQ=255)                                               ggCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCaaaa                                                                                               >  2:287135/1‑149 (MQ=255)                                                gCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCAGGCTGATGGCCACCCCCAGGCCAAAac                                                                                              >  3:258756/1‑149 (MQ=255)                                                 cccGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAac                                                                                              >  3:157767/1‑148 (MQ=255)                                                   cGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                           <  2:352213/149‑1 (MQ=255)                                                    gTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAg                                                                                                      >  5:54042/1‑137 (MQ=255)                                                    gTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAg                                                                                                      <  6:54042/137‑1 (MQ=255)                                                     tctgTTCGCGGTTGCTGCGCCCCTTTCCGCCGGCGGGCAAACGATTCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCTGGCTGATGGCCAACCCCAGGCCAAAACCAC‑Cg                                                                                        <  4:258756/149‑1 (MQ=255)                                                         ttCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                           <  1:218921/143‑1 (MQ=255)                                                         ttCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                           >  2:218921/1‑143 (MQ=255)                                                         ttCGCGGTTGCTGCGCCCCTATCCGCCCGCGGGCAAACGAATCCTGAAATGCGCCCCACACAGCCCCGACTTCGCCAGGGTCAACGTCCCACCCTGCCCCACGAACGCCCCCCGGCTGATGGCCCACCCCAGGCAAAAACCAC‑CGgtgt                                                                                    >  6:236388/1‑149 (MQ=255)                                                            gcgGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGc                                                                                 >  1:90191/1‑149 (MQ=255)                                                               gTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                           >  2:144066/1‑137 (MQ=255)                                                               gTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAaccac                                                                                           <  1:144066/137‑1 (MQ=255)                                                                       gcCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGACCCACCCAGCGCGGATTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCTGCTGATGGCCAACCCCAGGCCAAAAACAC‑CGGTGTTGCGATCGCGACt                                                                       <  5:236388/148‑1 (MQ=255)                                                                         cccTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCg                                                                    <  3:215655/149‑1 (MQ=255)                                                                         cccTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCg                                                                    >  4:215655/1‑149 (MQ=255)                                                                           cTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGt                                                                  >  1:303204/1‑149 (MQ=255)                                                                              tCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGTGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTc                                                                 >  6:96053/1‑147 (MQ=255)                                                                              tCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTc                                                                 <  2:303204/147‑1 (MQ=255)                                                                                cgccgGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGa                                                             <  7:106563/149‑1 (MQ=255)                                                                                cgccgGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGa                                                             <  4:255035/149‑1 (MQ=255)                                                                                      cggGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                       <  6:27238/149‑1 (MQ=255)                                                                                      cggGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                       <  5:255476/149‑1 (MQ=255)                                                                                               gAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTg                                               >  8:168263/1‑148 (MQ=255)                                                                                                  tCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGaa                                            >  6:89211/1‑148 (MQ=255)                                                                                                    cTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCCAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAAc                                         <  1:287135/149‑1 (MQ=255)                                                                                                         cTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCACTCGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGc                                     <  1:96617/149‑1 (MQ=255)                                                                                                         cTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTa                                                        >  1:234734/1‑129 (MQ=255)                                                                                                         cTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTa                                                        <  2:234734/129‑1 (MQ=255)                                                                                                            cgcCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTcgc                                 >  6:279942/1‑149 (MQ=255)                                                                                                              ccccACCCAGCGCCGACTGCGCCAGGGTCAAGGTTCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAAc                                         <  3:70520/139‑1 (MQ=255)                                                                                                              ccccACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGTAAAAc                                         >  4:70520/1‑139 (MQ=255)                                                                                                              ccccACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCgcgc                               <  7:168263/149‑1 (MQ=255)                                                                                                               cccacccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGct                              >  1:320005/1‑149 (MQ=255)                                                                                                                ccacccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctc                             <  5:279942/149‑1 (MQ=255)                                                                                                                 cacccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctcc                            >  4:189842/1‑149 (MQ=255)                                                                                                                    ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                       >  5:237055/1‑119 (MQ=255)                                                                                                                    ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAg                                                       <  6:237055/119‑1 (MQ=255)                                                                                                                    ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGACCGGCTGGAAAACCCGCTCGCGctcctc                          >  5:189642/1‑148 (MQ=255)                                                                                                                    ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCtcctcc                         <  2:90191/149‑1 (MQ=255)                                                                                                                    ccaGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCccctcc                         >  8:299590/1‑149 (MQ=255)                                                                                                                         gcCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAgg                    <  2:320005/149‑1 (MQ=255)                                                                                                                          ccGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAggg                   >  7:214756/1‑149 (MQ=255)                                                                                                                                  tcgaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGa                                                             >  5:173484/4‑99 (MQ=38)                                                                                                                                  tcgaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGa                                                             <  6:173484/96‑1 (MQ=38)                                                                                                                                  gcCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAAt           >  4:183002/1‑149 (MQ=255)                                                                                                                                   ccAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATc          <  4:157767/149‑1 (MQ=255)                                                                                                                                   ccAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCCCCTCCACAGCGACGCCAATc          <  7:299590/149‑1 (MQ=255)                                                                                                                                    gaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctcctc                          >  7:37015/2‑132 (MQ=255)                                                                                                                                    gaGGGGCAGGGTGGCACCTTGACCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGctcctc                          <  8:37015/131‑1 (MQ=255)                                                                                                                                    cAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAg                     >  6:311444/1‑137 (MQ=255)                                                                                                                                    cAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAg                     <  5:311444/137‑1 (MQ=255)                                                                                                                                    cAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATcc         >  2:118367/1‑149 (MQ=255)                                                                                                                                         tCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCa    <  5:96053/149‑1 (MQ=255)                                                                                                                                         tCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCa    <  5:89211/149‑1 (MQ=255)                                                                                                                                         tCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACACGGACGCCAATCCCGTCa    >  3:263468/1‑149 (MQ=255)                                                                                                                                           aaGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCAtc  <  8:214756/149‑1 (MQ=255)                                                                                                                                             ggTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATc          <  6:180935/139‑1 (MQ=255)                                                                                                                                             ggTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATc          >  5:180935/1‑139 (MQ=255)                                                                                                                                                 |                                                                                                                                                ACCGCTTTGTGTCGTGGTCCAGCGAGCTGGCAGAGCAGGATTAGCAGGCCCGTCTCTTCGCGGTTGCTGCGCCCCTATCCGCCGGCGGGCAAACGAATCCTGAACTGCGCCCCACCCAGCGCCGACTGCGCCAGGGTCAAGGTGCCACCCTGCCCCTCGATCGCCCGCCGGCTGATGGCCAACCCCAGGCCAAAACCAC‑CGGTGTTGCGATCGCGACTGCGGTCCAGACGGTAGAACGGCTGGAAAACCCGCTCGCGCTCCTCCACAGGGACGCCAATCCCGTCATC  >  NC_002947/3912703‑3912989

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