Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	522,571	G→T	6.3%	L52L (CTC→CTA)	PP_0435 ←	M23/M37 family peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	522,571	0	G	T	6.3%	220.7 / 3.4	70	L52L (CTC→CTA)	PP_0435	M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (34/31);  new base T (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

GCTCAGCGGCCGGGGCCTGTTGATCGGCCTTCTGGGTTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCGGGGGCGCTTTAGGCGGTTCGTTG  >  NC_002947/522426‑522715                                                                                                                                                  |                                                                                                                                                 gCTCAGCGGCCGGGGCCTGTTGATCGGCCTTCTGGGTTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGc                                                                                                                                               <  3:36075/149‑1 (MQ=255)     aGCGGCCGGGGCCTGTTGATCGGCCTTCTGGGTTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGa                                                                                                                                           <  5:181595/149‑1 (MQ=255)       cGGCCGGGGCCTGTTGATCGGCCTTCTGGGTTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGTGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAgg                                                                                                                                         <  1:97971/149‑1 (MQ=255)          ccGGGGCCTGTTGATCGGCCTTCTGGGTTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGtt                                                                                                                                       >  5:139314/1‑148 (MQ=255)           cGGGGCCTGATGATCGGCCTTCTGGGTTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGAGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGa                                                                                                                                     <  6:150337/149‑1 (MQ=255)                      gATCGGCCTTCTGGGTTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTc                                                                                                                          >  3:162923/1‑149 (MQ=255)                         cGGCCTTCTGGGTTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCtt                                                                                                                        <  5:136866/148‑1 (MQ=255)                              ttCTGGGTTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTc                                                                                                                  >  3:93819/1‑149 (MQ=255)                              ttCTGGGTTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTc                                                                                                                  >  3:100101/1‑149 (MQ=255)                                cTGGGTTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTc                                                                                                                  >  5:257254/1‑147 (MQ=255)                                cTGGGGTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTc                                                                                                                  <  6:257254/147‑1 (MQ=255)                                     tttCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCTAGCTCGAGGTTGAg                                                                                                                                    <  1:33452/124‑1 (MQ=255)                                     tttCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCTAGCTCGAGGTTGAg                                                                                                                                    >  2:33452/1‑124 (MQ=255)                                     tttCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCg                                                                                                           >  6:222577/1‑149 (MQ=255)                                     tttCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCg                                                                                                           <  7:166072/149‑1 (MQ=255)                                         gggggTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTg                                                                                                        <  3:190470/148‑1 (MQ=255)                                          ggggTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGTGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTgga                                                                                                      <  5:273379/149‑1 (MQ=255)                                           ggggAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTggaa                                                                                                     <  4:99484/149‑1 (MQ=255)                                                gCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGaaa                                                                                                >  2:52494/1‑149 (MQ=255)                                                gCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGaaa                                                                                                >  1:42981/1‑149 (MQ=255)                                                gCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGaaa                                                                                                >  5:300444/1‑149 (MQ=255)                                                  ctctTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATa                                                                                              >  7:116165/1‑149 (MQ=255)                                                   tctTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATAc                                                                                             >  7:127464/1‑149 (MQ=255)                                                   tctTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATAc                                                                                             >  6:193864/1‑149 (MQ=255)                                                   tctTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATAc                                                                                             >  2:166370/1‑149 (MQ=255)                                                     ttCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTggaagg                                                                                                   >  1:233607/1‑141 (MQ=255)                                                     ttCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTggaagg                                                                                                   <  2:233607/141‑1 (MQ=255)                                                       cGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGCTTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCa                                                                                           <  7:200283/147‑1 (MQ=255)                                                       cGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAg                                                                                          >  3:108928/1‑148 (MQ=255)                                                              ggCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGcc                                                                                   <  1:166370/148‑1 (MQ=255)                                                              ggCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGc                                                                                    >  4:274313/1‑147 (MQ=255)                                                                         aGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAgg                                                                       <  2:42981/149‑1 (MQ=255)                                                                             ttGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGggcgg                                                                   >  6:374524/1‑149 (MQ=255)                                                                               gCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAgg                                                                       >  1:356892/1‑143 (MQ=255)                                                                               gCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAgg                                                                       <  2:356892/143‑1 (MQ=255)                                                                               gCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGggcgg                                                                   <  1:52494/147‑1 (MQ=255)                                                                                   aCCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGc                                                             >  8:327962/1‑149 (MQ=255)                                                                                        cGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGc                                                          <  6:139314/147‑1 (MQ=255)                                                                                         gCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTTCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTg                                                        <  8:160956/148‑1 (MQ=255)                                                                                          cTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTggc                                                      <  5:222577/149‑1 (MQ=255)                                                                                           tttGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTggcg                                                     >  3:54375/1‑149 (MQ=255)                                                                                              gCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTggcggc                                                   <  6:300444/148‑1 (MQ=255)                                                                                               cACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGcc                                                  >  3:372404/1‑148 (MQ=255)                                                                                                aCCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTggaagg                                                                                                   <  8:69303/98‑1 (MQ=255)                                                                                                aCCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTggaagg                                                                                                   >  7:69303/1‑98 (MQ=255)                                                                                                  cAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCaa                                                >  7:85650/1‑147 (MQ=255)                                                                                                   aGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCaa                                                >  4:114061/1‑146 (MQ=255)                                                                                                    ggggTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGa                                            <  4:93819/149‑1 (MQ=255)                                                                                                      ggTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGCGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGCCGCAGATGCCGCTGGAGGCCAACAGATg                                          >  3:239751/1‑149 (MQ=255)                                                                                                             tgcgGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCtt                                    <  7:236441/147‑1 (MQ=255)                                                                                                              gcgGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCg                                  <  4:108928/149‑1 (MQ=255)                                                                                                               cgGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCAAGCTCGAGGTTGAGGGTAGTTTTATTGGCTTCAACTTCCCTGGAAGGAAATACCAAAAAACCAAGGATGAGCAGGGCGGCAATGCCGCTGCCGACAAAAAGATggatataaa                                  >  8:236441/1‑141 (MQ=255)                                                                                                                  gATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATa                              <  4:162923/149‑1 (MQ=255)                                                                                                                  gATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCAAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGAt                               >  5:278332/1‑148 (MQ=255)                                                                                                                    ttCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCGGGCGGCCAACAGATGGCTCTTCGGATaaa                            <  4:100101/149‑1 (MQ=255)                                                                                                                       gTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCg                         >  2:44974/1‑149 (MQ=255)                                                                                                                       gTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCg                         >  1:351070/1‑149 (MQ=255)                                                                                                                       gTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCg                         >  4:249851/1‑149 (MQ=255)                                                                                                                            tCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCGGGGgc                    >  5:161677/1‑149 (MQ=255)                                                                                                                            tCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCGGGGgc                    >  6:139154/1‑149 (MQ=255)                                                                                                                             cAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCGGGGgc                    <  3:274313/148‑1 (MQ=255)                                                                                                                             cAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCGGGGgc                    <  3:249851/148‑1 (MQ=255)                                                                                                                                 tgctCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCGGGGGCGCtt                >  8:343541/1‑148 (MQ=255)                                                                                                                                   ctCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAg                                                                                 <  3:2933/81‑1 (MQ=255)                                                                                                                                   ctCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAg                                                                                 >  4:2933/1‑81 (MQ=255)                                                                                                                                         aGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCGGGGGCGCTTTAGGCGGtt       <  3:114061/149‑1 (MQ=255)                                                                                                                                            gCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCGGGGGCGCTTTAGGCGGTTCGt    <  5:193864/149‑1 (MQ=255)                                                                                                                                              tttCTAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCGGGGgcgc                  <  3:287920/133‑1 (MQ=255)                                                                                                                                              tttCTAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCGGGGgcgc                  >  4:287920/1‑133 (MQ=255)                                                                                                                                              tttCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCGGGGGCGCTTTAGGCGGTTCGTTg  >  2:139377/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                 GCTCAGCGGCCGGGGCCTGTTGATCGGCCTTCTGGGTTTCGGGGGTAGCCTCTTCGATCTGGGCGAAAGGCGAGCCTTGCTCACCTTCGCTTTGCACCAGGGGTGCGGCGCGGGATTCGTCTTTCAGCTGCTCGGCAGGGCTTTCGAGCTCGAGGTTGAGGGTAGTTTTCTTGGCTTCCACTTCGCTGGAAGGAAATACCAGCAAAGCCAGGCTGAGCAGGGCGGCGATGCCGCTGGCGGCCAACAGATGGCTCTTCGGATAAAGCGGGGGCGCTTTAGGCGGTTCGTTG  >  NC_002947/522426‑522715

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