Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	150,829	G→T	8.2%	L41M (CTG→ATG) ‡	PP_0142 ←	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	150,829	0	G	T	8.2%	191.9 / 6.6	61	L41M (CTG→ATG) ‡	PP_0142	ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (23/33);  new base T (2/3);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

GCCTGGCTGGCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGG  >  NC_002947/150684‑150970                                                                                                                                                  |                                                                                                                                               gcCTGGCTGGCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGgcg                                                                                                                                             >  7:270514/1‑149 (MQ=255)   ctggctggCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGcgct                                                                                                                                           >  6:315411/1‑149 (MQ=255)      gctggctTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGcgctc                                                                                                                                          <  7:90254/147‑1 (MQ=255)        tggctTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGTTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCgctcgctc                                                                                                                                      >  1:357252/1‑149 (MQ=255)          gctTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCCTGcgctcgcatt                                                                                                                                     <  4:122351/148‑4 (MQ=255)          gctTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCCTGcgctcgcat                                                                                                                                      <  6:126466/147‑3 (MQ=255)           ctTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTcg                                                                                                                                   >  7:254203/1‑149 (MQ=255)            ttCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCtctc                                                                                                                                    <  5:211669/147‑1 (MQ=255)            ttCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCtctc                                                                                                                                    >  6:211669/1‑147 (MQ=255)                   aGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCCTGcgctcgcatt                                                                                                                                     >  6:47784/1‑136 (MQ=255)                   aGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCCTGcgctcgcatt                                                                                                                                     <  5:47784/139‑4 (MQ=255)                       cgTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTg                                                                                                                       <  3:218423/149‑1 (MQ=255)                       cgTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTg                                                                                                                       <  3:229038/149‑1 (MQ=255)                           ggTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCg                                                                                                                    <  4:346879/148‑1 (MQ=255)                            gTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTa                                                                                                                  <  1:228157/149‑1 (MQ=255)                              gCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAgg                                                                                                                                               >  4:8506/1‑118 (MQ=255)                              gCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAgg                                                                                                                                               <  3:8506/118‑1 (MQ=255)                              gCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAgt                                                                                                                >  7:208193/1‑149 (MQ=255)                                  tGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGc                                                                                                            <  8:99260/149‑1 (MQ=255)                                         ggCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGTTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTc                                                                                                     <  2:357252/149‑1 (MQ=255)                                         ggCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTc                                                                                                     <  4:141928/149‑1 (MQ=255)                                           cgcTCGGAGCCTAGCAGTTCTGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCa                                                                                                   <  6:112602/149‑1 (MQ=255)                                           cgcTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCa                                                                                                   <  8:270514/149‑1 (MQ=255)                                                    ccTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCgg                                                                                           >  3:98994/1‑148 (MQ=255)                                                      tAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAgtgt                                                                                                              <  6:252863/127‑1 (MQ=255)                                                      tAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAgtgt                                                                                                              >  5:252863/1‑127 (MQ=255)                                                           gTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACg                                                                                    >  5:33486/1‑148 (MQ=255)                                                             tCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTAAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACgca‑g                                                                                >  8:225472/1‑149 (MQ=255)                                                              cGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACgca‑g                                                                                <  8:208193/148‑1 (MQ=255)                                                                  cagcagCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGca‑ggcagg                                                                           >  3:131233/1‑149 (MQ=255)                                                                   agcagcGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCa‑ggcaggc                                                                          >  3:226368/1‑149 (MQ=255)                                                                   agcagcGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCa‑ggcaggc                                                                          <  4:131233/149‑1 (MQ=255)                                                                          aaGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTgg                                                                   >  2:263681/1‑149 (MQ=255)                                                                               ccGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGgctgct                                                              >  4:265390/1‑149 (MQ=255)                                                                                cGGCAGTGTCCAGTCGGCCCAGCTGGCTGAGGTCGTCGACGTTGTCGGCAGGGCAATGCGAGCGCAGGCTCTCGCTCTCGCGCTTGAGTTTTGCGTAGTGTGCCAGGGTCCAGTCACAGGTAATACGCA‑GTCAGGGCGGTTGGctgctg                                                             <  4:226368/149‑1 (MQ=255)                                                                                cGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGctgctg                                                             <  8:254203/149‑1 (MQ=255)                                                                                cGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGctgctg                                                             >  6:138887/1‑149 (MQ=255)                                                                                  gCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGgt                                                           <  5:138887/149‑1 (MQ=255)                                                                                   cAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGgtg                                                          <  4:98994/149‑1 (MQ=255)                                                                                             gcggcCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTgg                                                >  3:324874/1‑149 (MQ=255)                                                                                                gcCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGgcgc                                             >  2:9222/1‑149 (MQ=255)                                                                                                  ccAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTg                                           <  5:315411/149‑1 (MQ=255)                                                                                                    aGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTAAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTggg                                         <  7:225472/149‑1 (MQ=255)                                                                                                         gctgAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTg                                     <  1:9222/148‑1 (MQ=255)                                                                                                          ctgAGGTCGGCGACGGTGTCGGCAGGGCAATGCGTGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCACGGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTgg                                    <  1:263681/149‑1 (MQ=255)                                                                                                           tgAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAAGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGt                                   >  7:93460/1‑148 (MQ=255)                                                                                                              ggTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAAGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATa                               <  8:93460/149‑1 (MQ=255)                                                                                                                  ggCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCt                           <  6:33486/149‑1 (MQ=255)                                                                                                                    cGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTc                          >  3:120156/1‑148 (MQ=255)                                                                                                                     gACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCa                                                                                                          <  3:97874/68‑1 (MQ=255)                                                                                                                     gACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCa                                                                                                          >  4:97874/1‑68 (MQ=255)                                                                                                                       cGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATg                      <  7:205369/149‑1 (MQ=255)                                                                                                                            tCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCTGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAAtt                 <  7:71052/149‑1 (MQ=255)                                                                                                                            tCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAgg                             >  5:349025/1‑137 (MQ=255)                                                                                                                            tCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAgg                             <  6:349025/137‑1 (MQ=255)                                                                                                                            tCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAAtt                 <  4:120156/149‑1 (MQ=255)                                                                                                                            tCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGGTGGCGTAGTGTGCCAGGGTCCAGTCACCGGGAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAAtt                 <  2:308870/149‑1 (MQ=255)                                                                                                                               gCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAg                              <  3:15660/133‑1 (MQ=255)                                                                                                                               gCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAg                              >  4:15660/1‑133 (MQ=255)                                                                                                                                    gCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGgc         <  4:324874/149‑1 (MQ=255)                                                                                                                                     gagagCGAGCGCCTGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGgcg        >  2:307499/5‑149 (MQ=255)                                                                                                                                          cGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTaa    <  3:265390/148‑1 (MQ=255)                                                                                                                                           gAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTCAGACGCA‑GGAAGGCCGGTTGGCTGCGGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGCCTCGGAGCGAAgg  >  5:76038/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                               GCCTGGCTGGCTTCAGGCAGCGCGTGGGTGCAGTGCGACAGGCGCTCGGAGCCTAGCAGTTCGGCCAGCAGCGAAGCGCCGGCAGTGTCCAGGCGGCCCAGCTGGCTGAGGTCGGCGACGGTGTCGGCAGGGCAATGCGAGCGCAGGCGCTCGCTCTCGCGCTTGAGGTTGGCGTAGTGTGCCAGGGTCCAGTCACCGGTAATACGCA‑GGCAGGCCGGTTGGCTGCTGGTGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGG  >  NC_002947/150684‑150970

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