Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,248,772	G→A	6.5%	P419S (CCG→TCG)	PP_3724 ←	acyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,248,772	0	G	A	6.5%	182.8 / 2.7	62	P419S (CCG→TCG)	PP_3724	acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base G (25/32);  new base A (2/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGG  >  NC_002947/4248629‑4248920                                                                                                                                                |                                                                                                                                                     cTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAg                                                                                                                                                 <  7:220557/149‑1 (MQ=255)     cccTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGa                                                                                                                                             <  4:238804/149‑1 (MQ=255)     cccTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGa                                                                                                                                             <  6:86223/149‑1 (MQ=255)       cTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCTCGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAgg                                                                                                                                           <  5:134200/149‑1 (MQ=255)       cTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAgg                                                                                                                                           <  1:208795/149‑1 (MQ=255)       cTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAgg                                                                                                                                           <  3:376629/149‑1 (MQ=255)       cTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAgg                                                                                                                                           <  2:246024/149‑1 (MQ=255)           tCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGAc                                                                                                                                       <  4:283220/149‑1 (MQ=255)           caaaagaaCCGCCAGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGAc                                                                                                                                       >  6:62254/7‑149 (MQ=255)                  aCCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGcc                                                                                                                                 >  5:304787/1‑148 (MQ=255)                   ccgccagCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGa                                                                                                                               <  5:62254/143‑1 (MQ=255)                    cGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGAt                                                                                                                              <  2:179310/149‑1 (MQ=255)                     gCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATc                                                                                                                             >  2:96270/1‑149 (MQ=255)                         gCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGc                                                                                                                         >  4:120382/1‑149 (MQ=255)                           ggTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCg                                                                                                                       <  1:133708/149‑1 (MQ=255)                             ttCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTg                                                                                                                     <  7:519/149‑1 (MQ=255)                                ttttGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCg                                                                                                                  >  1:376431/1‑149 (MQ=255)                                     gTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCaa                                                                                                             <  1:223556/149‑1 (MQ=255)                                     gTAAGCGACGATGCGCCCCCAGGCAGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCaa                                                                                                             >  1:360468/1‑149 (MQ=255)                                        cgCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACg                                                                                                           <  2:18433/147‑1 (MQ=255)                                           gacgaTGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCAt                                                                                                       >  7:119763/1‑149 (MQ=255)                                              gaTGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGc                                                                                                    <  4:68722/149‑1 (MQ=255)                                               aTGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCa                                                                                                   <  6:304787/149‑1 (MQ=255)                                                  cgcCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGt                                                                                                 <  5:144454/148‑1 (MQ=255)                                                   gcCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTa                                                                                                >  6:84642/1‑148 (MQ=255)                                                          ggccgccaggccgTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACg                                                                                                           >  4:657/1‑130 (MQ=255)                                                          ggccgccaggccgTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACg                                                                                                           <  3:657/130‑1 (MQ=255)                                                           gccgccaggccgTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGtt                                                                                                                >  4:189330/1‑124 (MQ=255)                                                           gccgccaggccgTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGcccc                                                                                         >  4:128038/1‑147 (MQ=255)                                                           gccgccaggccgTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGcccc                                                                                         <  3:128038/147‑1 (MQ=255)                                                            ccgccaggccgTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCcagc                                                                                      <  1:96270/149‑1 (MQ=255)                                                                           ccACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCAtt                                                                         <  7:256945/147‑1 (MQ=255)                                                                             acacCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCtt                                                                     <  2:376431/149‑1 (MQ=255)                                                                               acCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGc                                                                   >  6:8889/1‑149 (MQ=255)                                                                                 ccTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGt                                                                 >  1:27466/1‑149 (MQ=255)                                                                                         gCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGGAGCAGGTCGCATTCCTTGCGTTCGCCg                                                           >  5:340577/1‑147 (MQ=255)                                                                                         gCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCg                                                           <  6:340577/147‑1 (MQ=255)                                                                                             cGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAAc                                                     <  1:324760/149‑1 (MQ=255)                                                                                                 gCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAAcgccg                                                 >  4:263945/1‑149 (MQ=255)                                                                                                  ccccGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAAcgccgc                                                >  2:264607/1‑149 (MQ=255)                                                                                                      gcgcATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCggg                                            <  3:263945/149‑1 (MQ=255)                                                                                                       cgcATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGAGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAACCGCCGCCggg                                            >  5:231634/1‑148 (MQ=255)                                                                                                       cgcATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGAGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCggg                                            >  5:226284/1‑148 (MQ=255)                                                                                                        gcATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTa                                          <  5:84642/149‑1 (MQ=255)                                                                                                              aGGCGCTGACTCGCACTGGCTGTGAACTGTTCCGTGAATGCCGAGGTGACGAAGCCGATCTTGCCGTGGAGTTCAACTGCATCGCAGTACAGCCCCAGCAGCAGGt                                                                               <  2:271838/106‑1 (MQ=255)                                                                                                              aGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGt                                                                               >  1:271838/1‑106 (MQ=255)                                                                                                               ggCGCTCACGCGCAGTGGCAGGCAACAGTTCCGGGAATGCGGAGGTGACGAAGCCGATGTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTGGCCGGTAAACGCCGCCGGGTAGTccagc                                   <  2:360468/149‑1 (MQ=255)                                                                                                                   cTCACGCGCACTGGCAGGCAACAGTTCCGAGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTccagcacca                               <  6:226284/149‑1 (MQ=255)                                                                                                                   cTCACGCGCACTGGCAGGCAACAGTTCCGAGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTccagcacca                               <  6:231634/149‑1 (MQ=255)                                                                                                                     cACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTc                                                                        <  8:242523/106‑1 (MQ=255)                                                                                                                      aCGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCtt                                                                     >  2:283828/1‑108 (MQ=255)                                                                                                                      aCGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCtt                                                                     <  1:283828/108‑1 (MQ=255)                                                                                                                      aCGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTcca                                     <  2:93907/140‑1 (MQ=255)                                                                                                                      aCGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTcca                                     >  1:93907/1‑140 (MQ=255)                                                                                                                      aCGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTcca                                     <  4:352983/140‑1 (MQ=255)                                                                                                                      aCGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTACTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTcca                                     >  3:352983/1‑140 (MQ=255)                                                                                                                      aCGCGCACTCCCAGGCAACAGTTACGGGAAAGCCGAGGTGACGAAGCAGATCTTGCAGTGGCTTTCAACTGCATCGAAGTACAGAAACAGCAGCAGTTCGAATTCCTTGCGTTCGCCGGTAAACGCAGACAGATAGAcaagcaccagc                             >  8:341860/1‑148 (MQ=255)                                                                                                                        gcgcACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAGACGCCGCCGGGTAGTCCAGCACCAGCATc                          >  5:340835/1‑149 (MQ=255)                                                                                                                                 caggcaACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGccc                 <  4:112419/149‑1 (MQ=255)                                                                                                                                  aggcaACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCcg                >  1:287645/1‑149 (MQ=255)                                                                                                                                      aaCAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCgcggc            >  5:26064/1‑149 (MQ=255)                                                                                                                                           ttCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGcgccgc       >  2:36637/1‑149 (MQ=255)                                                                                                                                                ggAAAGCCGAGGTGACGAAGCCGATCTTGCCGTCGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTACTCCCGAACCAGCAACGCCGCGCCCGCGGCGCCGCTCAgg  >  2:326602/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     CTCGCCCTGGTCCATCAACCGCTGGCGGTTCTTTTGGTAAGCGACGATGCGCCCCCAGGCCGCCAGGCCGTCTTCCACACCCTGCAGGGCAGCGATGCCCCGCGCATGCAGGCGCTCACGCGCACTGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGG  >  NC_002947/4248629‑4248920

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