Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	614,495	A→T	6.9%	L252Q (CTG→CAG)	ispA ←	farnesyl diphosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	614,495	0	A	T	6.9%	225.9 / 7.4	72	L252Q (CTG→CAG)	ispA	farnesyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base T (2/3);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

GGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGCATC  >  NC_002947/614346‑614637                                                                                                                                                      |                                                                                                                                               ggctAGGTGAACCTGTTTACTTAGGGCGTTAGCGCATATAAGGCGCCAAGGGTCGCAGCGCCGCGGCCGTTTCGCCGAACCCTGGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCGTTGGCGGCTTCCAGGCCcagca                                                                                                                                                >  2:12391/4‑150 (MQ=255) ggATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCcagca                                                                                                                                                >  2:185621/1‑150 (MQ=255)          aaGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCgg                                                                                                                                         <  1:402211/148‑1 (MQ=255)          aaGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCgg                                                                                                                                         >  2:402211/1‑148 (MQ=255)           aGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCTGCTGGGcctggaa                                                                                                                                      <  1:6581/150‑6 (MQ=255)              tGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGt                                                                                                                                   <  1:506695/150‑1 (MQ=255)               gTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTc                                                                                                                                  <  2:325577/150‑1 (MQ=255)                tttAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCg                                                                                                                                 >  2:702628/1‑150 (MQ=255)                   aGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGtt                                                                                                                              <  1:987466/150‑1 (MQ=255)                   aGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGtt                                                                                                                              <  1:931603/150‑1 (MQ=255)                   aGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCACAGCAGGGCCgg                                                                                                                                         <  2:280676/139‑1 (MQ=255)                   aGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCACAGCAGGGCCgg                                                                                                                                         >  1:280676/1‑139 (MQ=255)                    gTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGttt                                                                                                                             >  1:735322/1‑150 (MQ=255)                      tACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGt                                                                                                                           <  1:554518/150‑1 (MQ=255)                             ttCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGt                                                                                                                               <  1:69203/139‑1 (MQ=255)                             ttCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGt                                                                                                                               >  2:69203/1‑139 (MQ=255)                                gACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTAt                                                                                                                 >  2:312842/1‑150 (MQ=255)                                gACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTAt                                                                                                                 >  1:1159698/1‑150 (MQ=255)                                   aaTATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTc                                                                                                              <  2:1960/150‑1 (MQ=255)                                         aCGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCAcgc                                                                                                                      >  2:204967/1‑136 (MQ=255)                                         aCGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCAcgc                                                                                                                      <  1:204967/136‑1 (MQ=255)                                           gcgcCAAGGCCCGCAGCGGCTCGGCCTTTTCGCGGAACCCGTGCAGTGCGGCCAGCGCCTGGGCGCGCAGTTCGATTGCATAGGCTTCGGCGGCTTCCAGGCCCAGCAGGGCCGGGTCGGTCGGATTGTCACGCGCTATGTCTGCACCtt                                                                                                      >  2:269892/1‑150 (MQ=255)                                                  ggCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGttt                                                                                               <  1:702628/150‑1 (MQ=255)                                                     ccGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGcc                                                                                            >  1:329860/1‑150 (MQ=255)                                                              cTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGt                                                                                                                           >  2:686660/1‑110 (MQ=255)                                                              cTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGt                                                                                                                           <  1:686660/110‑1 (MQ=255)                                                                 ggCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAgt                                                                                >  2:732789/1‑150 (MQ=255)                                                                  gCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTCTGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAgtg                                                                               >  1:142571/1‑150 (MQ=255)                                                                      tttCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGc                                                                           <  1:185621/150‑1 (MQ=255)                                                                        tCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCGTGGCGTTTGCCCAAGGTAGCAGTGTCGCtt                                                                         <  1:877514/150‑1 (MQ=255)                                                                           ccGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTcc                                                                      >  2:957966/1‑150 (MQ=255)                                                                                  cTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCgg                                                                                                                                         <  1:234653/76‑1 (MQ=255)                                                                                  cTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCgg                                                                                                                                         >  2:234653/1‑76 (MQ=255)                                                                                   ttGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGGTTTGAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAg                                                              <  1:12391/150‑1 (MQ=255)                                                                                   ttGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTc                                                                                                              >  1:950798/1‑102 (MQ=255)                                                                                   ttGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTc                                                                                                              <  2:950798/102‑1 (MQ=255)                                                                                   ttGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAg                                                              >  1:340238/1‑150 (MQ=255)                                                                                       aGTGCGGCCAGCGCCTGGTCGCGCCGTTCGATTGCATTGGCTTTGGCGTCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATg                                                          <  2:735322/150‑1 (MQ=255)                                                                                         tGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATgtc                                                        >  1:1116216/1‑150 (MQ=255)                                                                                                  cgcCTGGTCGCGCAGTTGGATTGCATAGGCTTTGGCGTCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCAc                                               <  2:498214/150‑1 (MQ=255)                                                                                                     cTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTg                                            <  2:142571/150‑1 (MQ=255)                                                                                                      tGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTACGAGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTgg                                           <  1:317156/150‑1 (MQ=255)                                                                                                           gcgcAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGc                                      <  1:732789/150‑1 (MQ=255)                                                                                                             gcAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCa                                    >  1:704203/1‑150 (MQ=255)                                                                                                                     aTTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTc                                                                 >  2:441913/1‑113 (MQ=255)                                                                                                                     aTTGCATAGGCTTGGGCGGCTGCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTc                                                                 <  1:441913/113‑1 (MQ=255)                                                                                                                      ttGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATgg                           <  1:312842/150‑1 (MQ=255)                                                                                                                         cATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCt                        <  1:957966/150‑1 (MQ=255)                                                                                                                          aTAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCg                                                                            >  2:332444/1‑97 (MQ=255)                                                                                                                          aTAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCg                                                                            <  1:332444/97‑1 (MQ=255)                                                                                                                          aTAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCAcc                                              >  2:769471/1‑127 (MQ=255)                                                                                                                          aTAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCAcc                                              <  1:769471/127‑1 (MQ=255)                                                                                                                            aGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCtgtg                     >  1:518176/1‑150 (MQ=255)                                                                                                                                tttGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATa                 <  2:420519/150‑1 (MQ=255)                                                                                                                                    gcggcTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGc                                                                                                   <  2:293793/64‑1 (MQ=255)                                                                                                                                    gcggcTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGc                                                                                                   >  1:293793/1‑64 (MQ=255)                                                                                                                                        ctacccggCCCTGCTGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCtgt                      <  2:872865/131‑1 (MQ=255)                                                                                                                                        ctacccggCCCTGCTGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCtgt                      >  1:872865/7‑137 (MQ=255)                                                                                                                                         ttCCAGGCCCGGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAg        >  1:573835/1‑150 (MQ=255)                                                                                                                                         ttCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAg        >  2:710535/1‑150 (MQ=255)                                                                                                                                             aGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGCa    >  2:993447/1‑150 (MQ=255)                                                                                                                                             aGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGCa    >  2:989236/1‑150 (MQ=255)                                                                                                                                             aGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGCa    <  2:704203/150‑1 (MQ=255)                                                                                                                                              ggCCCTGCTGGGCCTGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCa                                                                                  >  1:821174/1‑71 (MQ=255)                                                                                                                                              ggCCCTGCTGGGCCTGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCa                                                                                  <  2:821174/71‑1 (MQ=255)                                                                                                                                               gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCt                                             >  1:133204/1‑107 (MQ=255)                                                                                                                                               gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCt                                             <  2:133204/107‑1 (MQ=255)                                                                                                                                               gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCa         >  1:378217/1‑143 (MQ=255)                                                                                                                                               gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCa         <  2:378217/143‑1 (MQ=255)                                                                                                                                               gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAgg       >  1:49035/1‑145 (MQ=255)                                                                                                                                               gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAgg       <  2:49035/145‑1 (MQ=255)                                                                                                                                               gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGc     >  1:163451/1‑147 (MQ=255)                                                                                                                                               gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGc     <  2:163451/147‑1 (MQ=255)                                                                                                                                               gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGCATc  >  2:324336/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                               GGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGCATC  >  NC_002947/614346‑614637

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑