Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	424,649	T→C	7.7%	I221V (ATT→GTT)	PP_0349 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	424,649	0	T	C	7.7%	188.6 / 4.8	65	I221V (ATT→GTT)	PP_0349	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base C (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 = 5.62e-01

GAAACGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTACGGTGCCGAAT  >  NC_002947/424502‑424785                                                                                                                                                       |                                                                                                                                           gAAACGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTACGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATg                                                                                                                                           <  1:263939/149‑1 (MQ=255)  aaaCGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGa                                                                                                                                          <  8:37590/149‑1 (MQ=255)    aCGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCATTGCGGCCTGGCGCAACGAca                                                                                                                                        <  7:339060/149‑2 (MQ=255)        gATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACGATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGt                                                                                                                                    <  5:318402/149‑1 (MQ=255)        gATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGt                                                                                                                                    <  7:153749/149‑1 (MQ=255)           gAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCAGGCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGa                                                                                                                                          >  1:265095/1‑141 (MQ=255)           gAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGa                                                                                                                                          <  2:265095/140‑1 (MQ=255)              gTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGCAGGCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGa                                                                                                                                >  1:107253/1‑149 (MQ=255)               tCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAg                                                                                                                             >  8:87086/1‑149 (MQ=255)               tCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAg                                                                                                                             <  1:49299/149‑1 (MQ=255)               tCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAg                                                                                                                             >  6:334561/1‑149 (MQ=255)               tCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAg                                                                                                                             <  6:221746/149‑1 (MQ=255)                 cGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAgg                                                                                                                            <  3:294687/148‑1 (MQ=255)                  gggCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGc                                                                                                                          >  3:292862/1‑149 (MQ=255)                  gggCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCATTGCGGCCTGGCGCAACGAca                                                                                                                                        >  2:11270/1‑134 (MQ=255)                  gggCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGATGTCATTGCGGCCTGGCGCAACGAca                                                                                                                                        <  1:11270/135‑2 (MQ=255)                   ggCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCGGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCa                                                                                                                         <  3:331259/149‑1 (MQ=255)                         tGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGAt                                                                                                                    >  7:309317/1‑148 (MQ=255)                            ccGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTAcc                                                                                                                >  6:315274/1‑149 (MQ=255)                            ccGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTAcc                                                                                                                >  2:107712/1‑149 (MQ=255)                             cGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTAccc                                                                                                               <  6:325563/149‑1 (MQ=255)                              ggCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTAcccc                                                                                                              <  4:292862/149‑1 (MQ=255)                                 gcgGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGGCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGt                                                                                                           <  5:99204/149‑1 (MQ=255)                                     ccTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACg                                                                                                       >  4:61820/1‑149 (MQ=255)                                      cTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGc                                                                                                      <  3:82666/149‑1 (MQ=255)                                       tGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTAcc                                                                                                                <  6:173082/138‑1 (MQ=255)                                       tGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTAcc                                                                                                                >  5:173082/1‑138 (MQ=255)                                       tGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCt                                                                                                     >  5:287351/1‑149 (MQ=255)                                       tGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCATTGCGGCCTGGCGCAACGAca                                                                                                                                        <  3:319847/114‑2 (MQ=255)                                       tGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCATTGCGGCCTGGCGCAACGAca                                                                                                                                        >  4:319847/1‑113 (MQ=255)                                         ttCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCt                                                                                                     <  5:315274/147‑1 (MQ=255)                                            ggCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCg                                                                                                >  6:330968/1‑149 (MQ=255)                                              cAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGt                                                                                               >  7:29077/1‑148 (MQ=255)                                                  gTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCg                                                                                                                                  <  2:129606/109‑1 (MQ=255)                                                   tCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCg                                                                                                >  2:27064/1‑142 (MQ=255)                                                   tCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCg                                                                                                <  1:27064/142‑1 (MQ=255)                                                    ccA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTc                                                                                                          >  3:10468/1‑131 (MQ=255)                                                    ccA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTc                                                                                                          <  4:10468/131‑1 (MQ=255)                                                     ca‑gcagGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGcc                                                                                        >  1:1962/1‑148 (MQ=255)                                                          agGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCccaggc                                                                                   <  2:107253/149‑1 (MQ=255)                                                                     tGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑a                                                                      >  3:348930/1‑149 (MQ=255)                                                                          ccgccgcTCGTCGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGc                                                                   >  3:34830/5‑149 (MQ=255)                                                                               gcTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTTCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑CCGCCGAt                                                               >  6:197497/1‑148 (MQ=255)                                                                               gcTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGAt                                                               >  5:61348/1‑148 (MQ=255)                                                                                cTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCa                                                             <  5:330968/149‑1 (MQ=255)                                                                                 tCGACGAGCGGCGGGGCGTCGAGGTAGGGGGCAACCATCGCAGACAGGCTGTCGTTGCTCCAGGCCTCTATGACCAGGTCGGAGAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCGGATCAg                                                            <  6:366474/149‑1 (MQ=255)                                                                                   gacgaGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAgg                                                           <  1:247922/148‑1 (MQ=255)                                                                                      gaGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGGACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTtgt                                                       >  7:243188/1‑149 (MQ=255)                                                                                        gcggcgGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGc                                                     <  3:186110/149‑1 (MQ=255)                                                                                        gcggcgGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGc                                                     <  3:61820/149‑1 (MQ=255)                                                                                         cggcggGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCa                                                    <  4:366217/149‑1 (MQ=255)                                                                                           gcggGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAgg                                                  <  4:56309/149‑1 (MQ=255)                                                                                            cggGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCCGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGt                                                 >  3:58056/1‑149 (MQ=255)                                                                                                gCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAg                                             >  5:153130/1‑149 (MQ=255)                                                                                                      cgGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAgcg                                       >  3:154156/1‑149 (MQ=255)                                                                                                                   ccATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGAt                          >  8:6936/1‑149 (MQ=255)                                                                                                                        gCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCtt                     <  6:61348/149‑1 (MQ=255)                                                                                                                          aGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTtgt                   <  2:1962/149‑1 (MQ=255)                                                                                                                           gCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTtgtg                  >  5:58787/1‑149 (MQ=255)                                                                                                                           gCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCAGCGGGATTTCTtgtg                  <  4:350224/149‑1 (MQ=255)                                                                                                                                gCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTa             >  2:22552/1‑149 (MQ=255)                                                                                                                                   gTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGGACTGATTACCCCGGTCACGCTGACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTACgg          >  1:98940/1‑149 (MQ=255)                                                                                                                                    tCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCAAGCTCACCGTCAATGCCCAGGCCAGGGTGACCTTCCGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTACg           <  1:107712/149‑1 (MQ=255)                                                                                                                                      gTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTACGGTg        <  4:154156/148‑1 (MQ=255)                                                                                                                                           gcCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTACGGTGCCGAAt  <  4:361210/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                           GAAACGCGATGAAGTCCGGGCGCATGCCCGGCGCGGCCTGTTCGGCAATGTCCA‑GCAGGCGAGTGGCTGGAGC‑‑GCGCTCGACGAGCGGCGGGGCGTCGCGGTAGGGGGCAACCATCGCAGCCAGGCTGTCGTTGCGCCAGGCCGCAATGACCAGGTCGGACAGGGCACTGATTACCCCGGTCACGCTCACCGTCAATGCCCAGGCCAGGGTGACC‑‑ACGCCGATCAGGTTGTGCAGGTCGAGCCAGCGCAGGCGGCGGGATTTCTTGTGGCGTACGGTGCCGAAT  >  NC_002947/424502‑424785

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