Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,287,688	T→A	8.9%	L272* (TTG→TAG)	PP_2016 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,287,688	0	T	A	8.9%	183.9 / 6.9	67	L272* (TTG→TAG)	PP_2016	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/36);  new base A (4/2);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

GCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGTTCAC  >  NC_002947/2287552‑2287831                                                                                                                                            |                                                                                                                                                gCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCa                                                                                                                                     <  6:24951/149‑1 (MQ=255) gCGTACCCGAAGAAGATTCGCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCa                                                                                                                                     <  6:265682/149‑1 (MQ=255)   gTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCAACTGGTAGGCCGTgtacttgc                                                                                                                                   <  1:187721/149‑8 (MQ=255)     aCCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTa                                                                                                                                  >  1:55255/1‑148 (MQ=255)       ccGAAGTTGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTaaaa                                                                                                                               <  6:83439/149‑1 (MQ=255)         gaagaaGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAg                                                                                                                             <  3:78292/149‑1 (MQ=255)         gaagaaGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAg                                                                                                                             >  8:135475/1‑149 (MQ=255)         gaagaaGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATATTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAg                                                                                                                             >  3:228450/1‑149 (MQ=255)               gaTTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATTTTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTg                                                                                                                       <  1:93231/149‑1 (MQ=255)                 ttGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGca                                                                                                                     <  4:59878/149‑1 (MQ=255)                       gCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCg                                                                                                               >  6:52924/1‑149 (MQ=255)                          tCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTAc                                                                                                            >  2:304084/1‑149 (MQ=255)                             aCCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCg                                                                                                          >  1:39/1‑148 (MQ=255)                               cTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGaa                                                                                                        >  8:21398/1‑148 (MQ=255)                               cTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGaa                                                                                                        <  7:21398/148‑1 (MQ=255)                                   aCAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGGCCTGCACATGCGTACCGAACAAcc                                                                                                   <  6:70983/149‑1 (MQ=255)                                     aaGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTaaaa                                                                                                                               <  7:63673/119‑1 (MQ=255)                                     aaGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTaaaa                                                                                                                               >  8:63673/1‑119 (MQ=255)                                       ggACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCaa                                                                                               <  7:135475/149‑1 (MQ=255)                                                      gggCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAg                                                                                <  6:178627/149‑1 (MQ=255)                                                      gggCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAg                                                                                <  5:179086/149‑1 (MQ=255)                                                       ggCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAgg                                                                               <  8:110926/149‑1 (MQ=255)                                                        gCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGa                                                                              >  1:22474/1‑149 (MQ=255)                                                         cACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGAt                                                                             >  5:28537/1‑149 (MQ=255)                                                          aCGCTGCTGGGCATCAAGGGCCATTAACTTATCTTCGACACTGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCATTTGGCCGTGCTCCAGTAAAAATGACCTGCACATTTGTACCGAACAACCGCAAGTGATc                                                                                         <  7:56088/136‑1 (MQ=255)                                                          aCGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACTGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATc                                                                                         >  8:56088/1‑136 (MQ=255)                                                          aCGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGAt                                                                             >  7:175791/1‑148 (MQ=255)                                                                   ggCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGgc                                                                     <  6:365267/147‑1 (MQ=255)                                                                     cATCAAGGGCCAGTACCTGATCTTCGACTCCTGTGTGATCAATATTAGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACCTGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCg                                                                 <  4:228450/149‑1 (MQ=255)                                                                     cATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCTCAAGCCTCTCACGTAGGGCGAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGgcg                                                                    <  1:312112/149‑1 (MQ=255)                                                                      aTCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCg                                                                 >  4:368142/1‑148 (MQ=255)                                                                           gggCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTAc                                                           <  3:226510/149‑1 (MQ=255)                                                                               cAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGa                                                       >  5:368626/1‑149 (MQ=255)                                                                                 gTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATc                                                     <  6:368626/149‑1 (MQ=255)                                                                                  tACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATc                                                     >  4:82699/1‑148 (MQ=255)                                                                                     cTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATcgac                                                  >  2:363634/1‑148 (MQ=255)                                                                                      tGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGaacaa                                                                                                     >  5:276356/1‑96 (MQ=255)                                                                                      tGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGaacaa                                                                                                     <  6:276356/96‑1 (MQ=255)                                                                                        aTCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATcgacga                                                <  5:52924/147‑1 (MQ=255)                                                                                        aTCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATcgacga                                                <  2:55255/147‑1 (MQ=255)                                                                                              gACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCAcc                                        >  2:79817/1‑149 (MQ=255)                                                                                                cacCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTCCACATCCGTACCGAACAACCGCAAGTGATCTACCTAAAGGATTACCAGGCGCCGGAGAACCTGATCAACGACACGCACCt                                       >  7:245953/1‑148 (MQ=255)                                                                                                 acCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGa                                     <  1:362525/149‑1 (MQ=255)                                                                                                  ccGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGAc                                    <  3:182952/149‑1 (MQ=255)                                                                                                   cGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGa                                     <  4:325095/147‑1 (MQ=255)                                                                                                   cGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGa                                     >  3:325095/1‑147 (MQ=255)                                                                                                   cGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGAcc                                   >  5:274340/1‑149 (MQ=255)                                                                                                        tgATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCg                               <  3:297461/148‑1 (MQ=255)                                                                                                        tgATCAATATACGCAAG‑‑T‑AGACGGCCTACCTTTTGGCCGTGCTCCAGTAATAAGGACCTGCACATGCGTACCGAACAACCGCTAGTGATCTACCTCATGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGa                              <  7:170182/149‑1 (MQ=255)                                                                                                          aTCATTATTCGGAAG‑‑T‑ACTCGGTCTAGCAGTTGGCCGTGCTCCATTAAAAAGGAGCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGTATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAg                             <  8:245953/148‑1 (MQ=255)                                                                                                             aaTATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGt                         <  3:368142/149‑1 (MQ=255)                                                                                                               tatTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTc                       <  6:28537/149‑1 (MQ=255)                                                                                                                   cGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAg                    <  8:175791/148‑1 (MQ=255)                                                                                                                   cGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAg                    <  2:22474/148‑1 (MQ=255)                                                                                                                     cAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGAc                 >  5:84803/1‑149 (MQ=255)                                                                                                                       aG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACgaga                                              >  1:179588/1‑118 (MQ=255)                                                                                                                       aG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACgaga                                              <  2:179588/118‑1 (MQ=255)                                                                                                                        ctggagcacGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGAc                 >  2:30662/7‑149 (MQ=255)                                                                                                                         tggagcacGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCAAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGATGTTCGAGGAc                 >  1:109682/6‑148 (MQ=255)                                                                                                                             gcacGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATTCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCacac            >  8:143455/2‑149 (MQ=255)                                                                                                                              cacGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGCCGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACg           >  6:242463/1‑149 (MQ=255)                                                                                                                               acGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGa                     >  5:240056/1‑138 (MQ=255)                                                                                                                               acGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGa                     <  6:240056/138‑1 (MQ=255)                                                                                                                               acGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGc          <  1:30662/149‑1 (MQ=255)                                                                                                                                    cTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGtt     <  1:304084/149‑1 (MQ=255)                                                                                                                                      aCCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGTTCa   >  6:327807/1‑149 (MQ=255)                                                                                                                                       ccAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGTTCAc  <  2:39/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                GCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGTTCAC  >  NC_002947/2287552‑2287831

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