Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,672,934	T→G	11.8%	Q25H (CAA→CAC)	PP_4136 ←	LuxR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,672,934	0	T	G	11.8%	156.9 / 11.2	59	Q25H (CAA→CAC)	PP_4136	LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (29/23);  new base G (2/5);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

ATGGCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGCCACCTG  >  NC_002947/4672788‑4673082                                                                                                                                                   |                                                                                                                                                      aTGGCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTg                                                                                                                                                     >  4:313893/1‑149 (MQ=255)         gATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCa                                                                                                                                              <  5:157375/148‑1 (MQ=255)         gATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttg                                                                                                                                               <  2:228530/147‑2 (MQ=255)           tCACAGCGGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCAAGTGCTGttgtggt                                                                                                                                           <  1:253981/149‑6 (MQ=255)           tCACAGCAGCCACCCTGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATTAGCACACCACAACAGCAAGTGCTGttgtggt                                                                                                                                           <  5:181161/149‑6 (MQ=255)           tCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCTCGt                                                                                                                                           >  8:58965/1‑149 (MQ=255)                 cagcCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCCCACCACAACAGCACTTGCTGTTCACGTGGctt                                                                                                                                      >  6:177495/1‑146 (MQ=255)                    ccACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAgg                                                                                                                                  <  2:173893/149‑1 (MQ=255)                                gCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGc                                                                                                                       >  3:238323/1‑148 (MQ=255)                                       aaaTTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGc                                                                                                                 >  4:154873/1‑147 (MQ=255)                                          tttGCGGAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAggg                                                                                                            <  8:173696/149‑1 (MQ=255)                                          tttGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACt                                                                                                                                                       <  3:298658/106‑1 (MQ=255)                                          tttGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACt                                                                                                                                                       >  4:298658/1‑106 (MQ=255)                                          tttGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAggg                                                                                                            >  4:315353/1‑149 (MQ=255)                                          tttGCGCAGGATGATGGATACGTGAAAATTCCACGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGTCATAAGCACACCACAACAGCACTAGCAGTACACGAGGCGTCAGGCGAAGCGGGCCGTCCGGCAggg                                                                                                            >  2:153566/1‑149 (MQ=255)                                           ttGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGc                                                                                                                                        <  6:223953/120‑1 (MQ=255)                                           ttGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGc                                                                                                                                        >  5:223953/1‑120 (MQ=255)                                               gcAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTtgc                                                                                                       >  6:75517/1‑149 (MQ=255)                                                   gATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGt                                                                                                                                                 >  7:308824/1‑103 (MQ=255)                                                   gATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGt                                                                                                                                                 <  8:308824/103‑1 (MQ=255)                                                    aTGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCTCGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGa                                                                                                  <  7:58965/149‑1 (MQ=255)                                                     tGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAg                                                                                                 <  3:313893/149‑1 (MQ=255)                                                     tGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAg                                                                                                 <  3:154873/149‑1 (MQ=255)                                                          gATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTc                                                                                             >  4:134402/1‑148 (MQ=255)                                                                tGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGtt                                                                                      >  1:191941/1‑149 (MQ=255)                                                                 gAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTc                                                                                                                                     >  7:15322/1‑101 (MQ=255)                                                                 gAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTc                                                                                                                                     <  8:15322/101‑1 (MQ=255)                                                                      ttCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGtt                                                                                      <  4:20081/143‑1 (MQ=255)                                                                      ttCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGtt                                                                                      >  3:20081/1‑143 (MQ=255)                                                                            ctCGACTCCTTGCAATCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGtt                                                                         <  4:159481/148‑1 (MQ=255)                                                                             tCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACCACAGGACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGttt                                                                        <  8:3961/149‑1 (MQ=255)                                                                             tCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGttt                                                                        <  5:177495/149‑1 (MQ=255)                                                                                aCTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGaa                                                                     <  4:238323/149‑1 (MQ=255)                                                                                    cTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTc                                                                  >  2:215960/1‑148 (MQ=255)                                                                                     ttGCATTCAAGGATCTGGACGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCCTCAGGCGAACCGGGGAGTCCGCCAGCGGTTCCTGCGAGTGTCGAGAGTTGGCAGTC‑ATCTTTGAACTcaa                                                                >  4:203477/1‑147 (MQ=255)                                                                                         aTTCAAGGATCTGGCCGATCTCCGACGAGCTCTTGTCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACgg                                                             <  1:153566/148‑1 (MQ=255)                                                                                           tCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGt                                                           >  8:161118/1‑148 (MQ=255)                                                                                             aaGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGAt                                                        >  8:240666/1‑149 (MQ=255)                                                                                             aaGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGCTGAt                                                        >  8:191122/1‑149 (MQ=255)                                                                                               ggATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTCTCAGACAATGTTT‑AACTCCACGGCTGAtt                                                       >  8:279914/1‑148 (MQ=255)                                                                                                      gCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCgg                                                 <  5:75517/147‑1 (MQ=255)                                                                                                      gCCGATCTCCCACGAGATCTTGCCATAGGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGATGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGAATGACCgggg                                               >  6:141862/1‑149 (MQ=255)                                                                                                       ccGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGc                                              <  7:279914/149‑1 (MQ=255)                                                                                                               cctcGAGCTGTTGCCATAATCACACCACCTCAGGACTTGCTGTTCACGTGGCGTGAGGCGAACCTGGGCGTCCGCCACGGGTTGCTGGGAGTGTCTAGCGTTGTCTGTC‑ATCTTTGAACTCCACGGGTGATTGAGCGGGGCATtgctgc                                      <  3:203477/146‑1 (MQ=255)                                                                                                               ccACGAGCTCTTGCCTTAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGATCTACACGGGTGATTGAGCGGGGCATtgct                                        <  5:141862/147‑1 (MQ=255)                                                                                                               ccACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCa                                             <  2:191941/142‑1 (MQ=255)                                                                                                               ccACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATtgctgc                                      >  2:89612/1‑149 (MQ=255)                                                                                                                                  gCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTtg                                                                                                        >  1:65475/1‑65 (MQ=255)                                                                                                                                  gCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTtg                                                                                                        <  2:65475/65‑1 (MQ=255)                                                                                                                                      accacAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAg               >  1:10014/1‑149 (MQ=255)                                                                                                                                        gtgaacaGCAAGTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTAc                                  <  6:152346/125‑1 (MQ=255)                                                                                                                                        gtgaacaGCAAGTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGaa                <  8:105036/143‑1 (MQ=255)                                                                                                                                        gtgaacaGCAAGTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACATGCAGTCACCTTACGaa                >  7:105036/4‑146 (MQ=255)                                                                                                                                        cacaACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGca             >  7:278150/1‑149 (MQ=255)                                                                                                                                        cacaACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGca             >  5:336435/1‑149 (MQ=255)                                                                                                                                        cacaACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGca             >  3:6216/1‑149 (MQ=255)                                                                                                                                         acaacaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATtgctgc                                      <  6:263882/123‑1 (MQ=255)                                                                                                                                         acaacaGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATtgctgc                                      >  5:263882/1‑123 (MQ=255)                                                                                                                                             caGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGc        <  7:240666/149‑1 (MQ=255)                                                                                                                                             caGCACGTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGCGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGc        >  5:232213/1‑149 (MQ=255)                                                                                                                                                cACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGCCAc     >  6:126400/1‑149 (MQ=255)                                                                                                                                                   ttGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGCCACCTg  >  4:224511/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      ATGGCCTTGATCACAGCAGCCACCCGGGTAGGCACTTCAAATTTGCGCAGGATGTTGGATACGTGAAAATTCACCGTCGACTCCTTGCATTCAAGGATCTGGCCGATCTCCCACGAGCTCTTGCCATAAGCACACCACAACAGCACTTGCTGTTCACGTGGCGTCAGGCGAACCGGGGCGTCCGCCAGGGGTTGCTGCGAGTGTCGAGCGTTGTCAGTC‑ATGTTTGAACTCCACGGGTGATTGAGCGGGGCATTGCTGCCTACCTGCAGTCACCTTACGAAGCACAGGCCACCTG  >  NC_002947/4672788‑4673082

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