Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	614,492	A→T	7.0%	L253Q (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,492	0	A	T	7.0%	222.4 / 5.5	72	L253Q (CTG→CAG)	ispA	farnesyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (32/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 = 9.66e-01

GGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGCATC  >  NC_002947/614345‑614637                                                                                                                                                    |                                                                                                                                                  gggATAGGTGAAGCTGGTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCGGGACGCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCcagc                                                                                                                                                 >  1:420519/1‑150 (MQ=255)  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:931603/150‑1 (MQ=255)                    aGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGtt                                                                                                                              <  1:987466/150‑1 (MQ=255)                    aGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCACAGCAGGGCCgg                                                                                                                                         >  1:280676/1‑139 (MQ=255)                    aGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCACAGCAGGGCCgg                                                                                                                                         <  2:280676/139‑1 (MQ=255)                     gTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGttt                                                                                                                             >  1:735322/1‑150 (MQ=255)                       tACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGt                                                                                                                           <  1:554518/150‑1 (MQ=255)                              ttCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGt                                                                                                                               >  2:69203/1‑139 (MQ=255)                              ttCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGt                                                                                                                               <  1:69203/139‑1 (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                                                                                                                           <  1:686660/110‑1 (MQ=255)                                                               cTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGt                                                                                                                           >  2:686660/1‑110 (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                                                                                                                                         >  2:234653/1‑76 (MQ=255)                                                                                   cTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCgg                                                                                                                                         <  1:234653/76‑1 (MQ=255)                                                                                    ttGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGGTTTGAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAg                                                              <  1:12391/150‑1 (MQ=255)                                                                                    ttGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTc                                                                                                              <  2:950798/102‑1 (MQ=255)                                                                                    ttGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTc                                                                                                              >  1:950798/1‑102 (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                                             <  2:133204/107‑1 (MQ=255)                                                                                                                                                gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCt                                             >  1:133204/1‑107 (MQ=255)                                                                                                                                                gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCa         <  2:378217/143‑1 (MQ=255)                                                                                                                                                gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCa         >  1:378217/1‑143 (MQ=255)                                                                                                                                                gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAgg       <  2:49035/145‑1 (MQ=255)                                                                                                                                                gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAgg       >  1:49035/1‑145 (MQ=255)                                                                                                                                                gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGc     <  2:163451/147‑1 (MQ=255)                                                                                                                                                gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGc     >  1:163451/1‑147 (MQ=255)                                                                                                                                                gCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGCATc  >  2:324336/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTGCATAGGCTTTGGCGGCTTCCAGGCCCAGCAGGGCCGGGTAGGTCGGTTTGTCACGCGCTATGTCTGCACCTTGGCGTTTGCCCAAGGTAGCAGTGTCGCTTTCCACGTCGAGGATGTCGTCCTGCACCTGGAATGCCAGGCCAATGGCCTGTGCATAAGTTTGCAGGGCATC  >  NC_002947/614345‑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: ‑