Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,928,996	G→A	5.2%	R366C (CGC→TGC)	PP_1728 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,928,996	0	G	A	5.2%	230.4 / 2.2	77	R366C (CGC→TGC)	PP_1728	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (38/35);  new base A (2/2);  total (40/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 = 5.34e-01

CCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAGGGTACGTCGGGTCAGGCGTCGCGACAGTGGCT  >  NC_002947/1928852‑1929142                                                                                                                                                 |                                                                                                                                                   ccAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAgg                                                                                                                                                >  6:191667/1‑149 (MQ=255)  cagcagGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAggg                                                                                                                                               >  6:328218/1‑149 (MQ=255)  cagcagGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAggg                                                                                                                                               >  6:322671/1‑149 (MQ=255)           gCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTg                                                                                                                                       <  5:191667/148‑1 (MQ=255)           gCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTg                                                                                                                                       <  7:169053/148‑1 (MQ=255)           gCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTg                                                                                                                                       <  8:42725/148‑1 (MQ=255)             cGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGgtg                                                                                                                                    >  2:97628/1‑149 (MQ=255)             cGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGgtg                                                                                                                                    >  6:245069/1‑149 (MQ=255)              gAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGCCGGTGCTGgtgt                                                                                                                                   >  8:262044/1‑149 (MQ=255)                 gTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGAGGTGCTGGTGTCGc                                                                                                                                >  7:140269/1‑149 (MQ=255)                      gcAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCcggcg                                                                                                                           >  6:86183/1‑149 (MQ=255)                        aGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGga                                                                                                                         <  4:257829/149‑1 (MQ=255)                                ggTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAgg                                                                                                                 >  8:303160/1‑149 (MQ=255)                                 gtgtACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCACAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGt                                                                                                                >  2:200842/1‑149 (MQ=255)                                 gtgtACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCACAGGGAGGTGCTGGTGTCGCCGGCGGAGACGAAGGt                                                                                                                >  6:196571/1‑149 (MQ=255)                                      cacCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGa                                                                                                           <  7:303160/149‑1 (MQ=255)                                       acCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGAc                                                                                                                      >  3:208029/1‑137 (MQ=255)                                       acCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGAc                                                                                                                      <  4:208029/137‑1 (MQ=255)                                        ccGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGaa                                                                                                                   <  8:328443/139‑1 (MQ=255)                                        ccGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGaa                                                                                                                   <  8:323744/139‑1 (MQ=255)                                        ccGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGaa                                                                                                                   >  7:323744/1‑139 (MQ=255)                                        ccGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGaa                                                                                                                   >  7:328443/1‑139 (MQ=255)                                        ccGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGaa                                                                                                          >  4:298732/1‑148 (MQ=255)                                         cGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGaa                                                                                                          >  1:109570/1‑147 (MQ=255)                                         cGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACa                                                                                                        <  5:337583/149‑1 (MQ=255)                                         cGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACa                                                                                                        <  4:330257/149‑1 (MQ=255)                                           ccTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGc                                                                                                      <  1:97628/149‑1 (MQ=255)                                                 gcTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACAGCAGCGCCTGGCGCAGGGCGGTGCTGGTGTCGCCGTCGGAGACGCAGGTGGCGGACAGCGAAATg                                                                                                >  7:116959/1‑149 (MQ=255)                                                     gCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGc                                                                                              <  5:245069/147‑1 (MQ=255)                                                        ggCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGCGCACCATGCACCCCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGCCGACGGTGGCGAACAGCGAAATGGccac                                                                                           >  8:277505/1‑147 (MQ=255)                                                           ccGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGaa                                                                                      >  4:216520/1‑149 (MQ=255)                                                                tGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAgc                                                                                 >  5:55459/1‑149 (MQ=255)                                                                 ggggTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAgcg                                                                                >  4:143454/1‑149 (MQ=255)                                                                     ttGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCg                                                                             >  7:83023/1‑148 (MQ=255)                                                                        gcggcATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAAt                                                                         >  6:10331/1‑149 (MQ=255)                                                                            cATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCACAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAAc                                                                     <  5:196571/149‑1 (MQ=255)                                                                            cATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCACAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAAc                                                                     <  1:200842/149‑1 (MQ=255)                                                                             aTGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAAcc                                                                    <  7:290004/149‑1 (MQ=255)                                                                                       aCCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGcc                                                          >  3:130309/1‑149 (MQ=255)                                                                                          aGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAg                                                       <  8:132873/149‑1 (MQ=255)                                                                                           ggCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAgg                                                      <  6:55459/149‑1 (MQ=255)                                                                                           ggCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAgg                                                      <  5:86183/149‑1 (MQ=255)                                                                                           ggCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAg                                                       >  8:107650/1‑148 (MQ=255)                                                                                            gCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGa                                                     <  2:109570/149‑1 (MQ=255)                                                                                                   aGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGAt                                              <  7:256761/149‑1 (MQ=255)                                                                                                    gCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATg                                             <  7:262044/149‑1 (MQ=255)                                                                                                    gCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATg                                             <  3:298732/149‑1 (MQ=255)                                                                                                    gCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATg                                             <  3:143454/149‑1 (MQ=255)                                                                                                        cATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTAt                                         >  6:20097/1‑149 (MQ=255)                                                                                                        cATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTAt                                         >  6:226017/1‑149 (MQ=255)                                                                                                           gCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGaa                                                                <  4:181292/123‑1 (MQ=255)                                                                                                           gCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGaa                                                                >  3:181292/1‑123 (MQ=255)                                                                                                           gCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGc                                      >  8:63290/1‑149 (MQ=255)                                                                                                           gCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGc                                      >  8:68276/1‑149 (MQ=255)                                                                                                           gCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGc                                      >  1:129383/1‑149 (MQ=255)                                                                                                            cacaCCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGc                                      <  5:328218/148‑1 (MQ=255)                                                                                                            cacaCCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGc                                      <  5:322671/148‑1 (MQ=255)                                                                                                             acacCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATTAACCAGAAGCTGCCGAGGAACTCGATGGTATCGCt                                     >  5:344848/1‑148 (MQ=255)                                                                                                              caccaccacCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGa                                   <  5:10331/149‑1 (MQ=255)                                                                                                              caccaccacCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGa                                   <  4:130309/149‑1 (MQ=255)                                                                                                               accaccaccAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAg                                  <  8:83023/149‑1 (MQ=255)                                                                                                                 caccaccaGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGaa                                                                >  6:371481/1‑117 (MQ=255)                                                                                                                 caccaccaGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGaa                                                                <  5:371481/117‑1 (MQ=255)                                                                                                                 caccaccaGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGa                                   >  4:80236/1‑146 (MQ=255)                                                                                                                 caccaccaGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGa                                   <  3:80236/146‑1 (MQ=255)                                                                                                                   ccaccaGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGAGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAGGGTa                              >  1:156299/1‑149 (MQ=255)                                                                                                                    caccaGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAGCCAGAAGCGGCCGAGGATCTCGATGGTATCGCGGAGGGTAc                             >  5:367907/1‑149 (MQ=255)                                                                                                                       caGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAgg                                 <  2:173783/142‑1 (MQ=255)                                                                                                                       caGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAgg                                 >  1:173783/1‑142 (MQ=255)                                                                                                                       caGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAGGGTACGTc                          >  4:273233/1‑149 (MQ=255)                                                                                                                        aGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAGGGTACGTCg                         <  7:107650/149‑1 (MQ=255)                                                                                                                         gCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTc                                                                                                                                  <  8:230083/43‑1 (MQ=255)                                                                                                                            cGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAGGGTACGTCGGGTc                     <  3:216520/149‑1 (MQ=255)                                                                                                                                    aTCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAGGGTACGTCGGGTCAGGCGTcg             <  3:273233/149‑1 (MQ=255)                                                                                                                                       aGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAGGGTACGTCGGGTCAGGCGTCGCGa          <  2:129383/149‑1 (MQ=255)                                                                                                                                        gCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAGGGTACGTCGGGTCAGGCGTCGCGa          >  5:112594/1‑148 (MQ=255)                                                                                                                                         cGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAGGGTACGTCGGGTCAGGCGTCGCGa          >  6:142429/1‑147 (MQ=255)                                                                                                                                               gcgcAGGGCGGTGCTGGTGTCGCCGGCGGACACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAGGGTACGTCGGGTCAGGCGTCGCGACAGTGGCt  >  6:220087/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAACCAGAAGCTGCCGAGGATCTCGATGGTATCGCTGAGGGTACGTCGGGTCAGGCGTCGCGACAGTGGCT  >  NC_002947/1928852‑1929142

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