Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,927,771	C→T	6.4%	S193S (AGC→AGT)	PP_1727 →	GntR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,927,771	0	C	T	6.4%	183.6 / 2.5	62	S193S (AGC→AGT)	PP_1727	GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base T (2/2);  total (30/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.01e-01

CCTCAACCCCCGGTACTTTCCAGGCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTGGATGAGGTGTTGCCTTC  >  NC_002947/1927627‑1927918                                                                                                                                                  |                                                                                                                                                     ccTCAACCCCCGGTACTTTCCAGGCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCt                                                                                                                                                  <  5:181768/149‑1 (MQ=255) ccTCAACCCCCGGTACTTTCCAGACATTCTCCAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGGGTGCCTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGcc                                                                                                                                                   <  1:78960/149‑1 (MQ=255)  cTCAACCCCCGGTACTTTCCAGGCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCtt                                                                                                                                                 >  5:104790/1‑149 (MQ=255)  cTCAACCCCCGGTACTTTCCAGGCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCtt                                                                                                                                                 <  1:237882/149‑1 (MQ=255)        ccccGGTACTTTCCAGGCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAgt                                                                                                                                           >  5:259464/1‑149 (MQ=255)         cccGGTACTTTCCAGGCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAgtg                                                                                                                                          >  5:67308/1‑149 (MQ=255)          ccGGTACTTTCCAGGCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAgtgt                                                                                                                                         >  4:63970/1‑149 (MQ=255)             gTACTTTCCAGGCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTcgg                                                                                                                                      >  5:124903/1‑149 (MQ=255)                        gCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGTGCCTTGAAAgtgt                                                                                                                                         >  5:25957/1‑135 (MQ=255)                        gCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGTGCCTTGAAAgtgt                                                                                                                                         <  6:25957/135‑1 (MQ=255)                        gCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGTGCCTTGAAAgtgt                                                                                                                                         <  6:31659/135‑1 (MQ=255)                        gCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGTGCCTTGAAAgtgt                                                                                                                                         >  5:31659/1‑135 (MQ=255)                          aTTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCggg                                                                                                                                  <  3:140204/140‑1 (MQ=255)                          aTTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCggg                                                                                                                                  >  4:140204/1‑140 (MQ=255)                           ttCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCgggg                                                                                                                        <  7:312944/149‑1 (MQ=255)                           ttCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCgggg                                                                                                                        <  1:190571/149‑1 (MQ=255)                                gAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGca                                                                                                                   <  5:318293/149‑1 (MQ=255)                                   cacGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGca                                                                                                                   >  6:69872/1‑146 (MQ=255)                                   cacGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGca                                                                                                                   <  5:69872/146‑1 (MQ=255)                                        ccTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTAccc                                                                                                           >  8:101445/1‑149 (MQ=255)                                            gcgcAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCgg                                                                                                         <  6:104790/147‑1 (MQ=255)                                            gcgcAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCgg                                                                                                         <  6:124903/147‑1 (MQ=255)                                            gcgcAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCgg                                                                                                         <  6:259464/147‑1 (MQ=255)                                            gcgcAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGt                                                                                                       <  6:67308/149‑1 (MQ=255)                                                aGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCaa                                                                                                    >  8:323062/1‑148 (MQ=255)                                                    aCTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑Gc                                                                                               >  3:207787/1‑148 (MQ=255)                                                          cGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAg                                                                                         >  2:73647/1‑148 (MQ=255)                                                            agaTTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCAt                                                                                      >  3:235511/1‑149 (MQ=255)                                                             gaTTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATg                                                                                     <  7:323062/149‑1 (MQ=255)                                                                     ggCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACtt                                                                              >  2:106419/1‑148 (MQ=255)                                                                             tATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACtg                                                                     >  8:60560/1‑149 (MQ=255)                                                                               tGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACtgt                                                                    >  1:380047/1‑148 (MQ=255)                                                                                 gCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGa                                                                            >  4:354723/1‑138 (MQ=255)                                                                                 gCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGa                                                                            <  3:354723/138‑1 (MQ=255)                                                                                      gAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGa                                                            <  3:63970/149‑1 (MQ=255)                                                                                       aGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGAc                                                                                            >  5:377949/1‑116 (MQ=255)                                                                                       aGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGAc                                                                                            <  6:377949/116‑1 (MQ=255)                                                                                           cGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTAtt                                                       <  7:101445/149‑1 (MQ=255)                                                                                             ggCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACACGCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTg                                                      <  1:73647/149‑1 (MQ=255)                                                                                             ggCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGa                                                                            <  1:178307/126‑1 (MQ=255)                                                                                             ggCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGa                                                                            >  2:178307/1‑126 (MQ=255)                                                                                             ggCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTgg                                                     <  2:380047/149‑1 (MQ=255)                                                                                                  gtgtG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCa                                                >  1:151355/1‑149 (MQ=255)                                                                                                  gtgtG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCa                                                >  7:110932/1‑149 (MQ=255)                                                                                                    gtg‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACg                                              <  7:60560/149‑1 (MQ=255)                                                                                                        cTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATg                                           <  8:131917/149‑1 (MQ=255)                                                                                                         ttCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAg                                                           >  4:364190/1‑132 (MQ=255)                                                                                                         ttCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAg                                                           <  3:364190/132‑1 (MQ=255)                                                                                                         ttCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGAGCGACTGTGACCTCGAGTATTGGCGGCACGATGc                                          >  5:206431/1‑149 (MQ=255)                                                                                                                  cTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATc                                      <  8:262084/144‑1 (MQ=255)                                                                                                                  cTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATc                                      >  7:262084/1‑144 (MQ=255)                                                                                                                       gACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTg                            >  3:122579/1‑149 (MQ=255)                                                                                                                       gACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTg                            >  3:235026/1‑149 (MQ=255)                                                                                                                          aGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGACCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGccc                         <  6:206431/149‑1 (MQ=255)                                                                                                                              ttGCCGGTGGAGGCGGAAAGCGCCTTGAAAGTGTAGGCGGTCAGCCCGGGGGTGCACATTACCCGGGTAAAGA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCTAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGc                     <  1:313286/149‑1 (MQ=255)                                                                                                                                gCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTg                   >  5:160386/1‑149 (MQ=255)                                                                                                                                   ggTGGAGGCGGCAAGCGCCTTGAAAGTGGCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTGGAt                <  4:235026/149‑1 (MQ=255)                                                                                                                                       gAGGCGGCAAGCGCATTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTGGATGAg             <  4:207787/148‑1 (MQ=255)                                                                                                                                          gcggcAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTgg                  <  6:310914/140‑1 (MQ=255)                                                                                                                                          gcggcAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTgg                  >  5:310914/1‑140 (MQ=255)                                                                                                                                               aaGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTGGATGAGGTGtt        <  6:167329/145‑1 (MQ=255)                                                                                                                                               aaGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTGGATGAGGTGtt        >  5:167329/1‑145 (MQ=255)                                                                                                                                                 gcgcCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTGGATGAGGTGTTGCcttc  >  5:458/1‑149 (MQ=255)                                                                                                                                                 gcgcCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTGGATGAGGTGTTGCcttc  <  8:110932/149‑1 (MQ=255)                                                                                                                                                  cgcCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTGGATGAg             >  2:342305/1‑137 (MQ=255)                                                                                                                                                  cgcCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTGGATGAg             <  1:342305/137‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CCTCAACCCCCGGTACTTTCCAGGCATTCTCGAACACGACCTGGCGCAGTCACTGACCGAGATTTATGGGCGGGTGTATGGCATCGAGTACGGGCAGGTGTG‑CTTCGAGATCCTGCCGACAGCGTTGCCGGTGGAGGCGGCAAGCGCCTTGAAAGTGTCGGCGGGCAGCCCGGGGTTGCACATTACCCGGGTCAACA‑GCGACCAGCATGGGCACTTGATCGACTGTGACCTCGAGTATTGGCGGCACGATGCGATCCGCATTCGTGCCCAGGCTGGATGAGGTGTTGCCTTC  >  NC_002947/1927627‑1927918

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