Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,884,980	T→G	6.2%	A53A (GCA→GCC) ‡	PP_5159 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,884,980	0	T	G	6.2%	227.8 / 2.9	74	A53A (GCA→GCC) ‡	PP_5159	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base G (3/2);  total (37/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 = 1.00e+00

CCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAG  >  NC_002947/5884840‑5885125                                                                                                                                             |                                                                                                                                                  ccTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGTGtt                                                                                                                                           <  1:138913/149‑1 (MQ=255) ccTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGtt                                                                                                                                           <  8:184407/149‑1 (MQ=255)   ttGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCCAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGtt                                                                                                                                           >  4:332702/1‑147 (MQ=255)      gTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCGCAGGGGTTgggc                                                                                                                                       <  7:329098/148‑4 (MQ=255)          aTTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGc                                                                                                                                   <  5:103539/148‑1 (MQ=255)           ttGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCt                                                                                                                                 >  4:247698/1‑149 (MQ=255)              ccAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCGCAGGGGTTgggcg                                                                                                                                      <  8:273904/141‑5 (MQ=255)                aGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGc                                                                                                                            <  3:62401/149‑1 (MQ=255)                 ggCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCa                                                                                                                           >  2:38027/1‑149 (MQ=255)                        cTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGgc                                                                                                                    >  1:341049/1‑149 (MQ=255)                                  gccgccAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCCCTc                                                                                                          >  4:221163/1‑149 (MQ=255)                                  gccgccAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTc                                                                                                          >  8:296272/1‑149 (MQ=255)                                        aGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTg                                                                                                     >  7:98985/1‑148 (MQ=255)                                                gTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgc                                                                                              <  3:36577/147‑1 (MQ=255)                                                 tCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgcgcg                                                                                           <  6:111980/149‑1 (MQ=255)                                                 gggCAGGCCGAGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTgcgcg                                                                                           >  1:162120/3‑149 (MQ=255)                                                   gCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCa                                                                                         >  2:17018/1‑149 (MQ=255)                                                   gCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCa                                                                                         >  3:267311/1‑149 (MQ=255)                                                   gCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCGGGGTAGCTCACCGCCAAACCCCGGCGGGGTTGAGTTAGCCGTAGGCATCCCGGATCCGGCACTCGGCTGCAAGTGCGCGCa                                                                                         >  3:246443/1‑149 (MQ=255)                                                        ccTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGtt                                                                                                                                           <  7:225809/94‑1 (MQ=255)                                                        ccTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGtt                                                                                                                                           >  8:225809/1‑94 (MQ=255)                                                            ggCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCCGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCActttctgc                                                                                 >  4:235670/1‑141 (MQ=255)                                                              ccTGCCCAGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCgg                                                                                >  1:272195/1‑147 (MQ=255)                                                                 gCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggcagg                                                                            >  2:245846/1‑148 (MQ=255)                                                                  ccccGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCggcaggc                                                                           <  3:247698/148‑1 (MQ=255)                                                                            aTCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTgcgc                                                                >  5:212256/1‑149 (MQ=255)                                                                            aTCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTgcgc                                                                <  1:38027/149‑1 (MQ=255)                                                                             tCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGcc                                                               <  8:98985/149‑1 (MQ=255)                                                                              ccACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCct                                                              >  3:173159/1‑149 (MQ=255)                                                                               cACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCctg                                                             >  6:43430/1‑149 (MQ=255)                                                                               cACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCctg                                                             >  3:179722/1‑149 (MQ=255)                                                                               cACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGAATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCctg                                                             >  1:360302/1‑149 (MQ=255)                                                                                 cTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCctgct                                                           >  5:327055/1‑149 (MQ=255)                                                                                 cTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCctgct                                                           >  8:88337/1‑149 (MQ=255)                                                                                   ggATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTgg                                                         >  5:134185/1‑149 (MQ=255)                                                                                    gATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCctgctg                                                          <  7:296272/147‑1 (MQ=255)                                                                                    gATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTgg                                                         <  7:88337/148‑1 (MQ=255)                                                                                    gATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTgg                                                         <  1:245846/148‑1 (MQ=255)                                                                                       cTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGAc                                                      <  2:341049/148‑1 (MQ=255)                                                                                        tGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGAc                                                      >  5:348899/1‑147 (MQ=255)                                                                                          gTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGa                                                       <  1:148589/144‑1 (MQ=255)                                                                                          gTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCa                                                    <  4:163148/147‑1 (MQ=255)                                                                                          gTCACCGCCAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCa                                                    <  3:332702/147‑1 (MQ=255)                                                                                           tCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGc                                                  >  7:34394/1‑148 (MQ=255)                                                                                             aCCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCAt                                               >  2:259474/1‑149 (MQ=255)                                                                                             aCCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCAt                                               <  6:348899/149‑1 (MQ=255)                                                                                                     tGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGcc                                       <  1:259474/149‑1 (MQ=255)                                                                                                             cgcACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGAtcgtc                               <  2:272195/149‑1 (MQ=255)                                                                                                                     ctGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGacgacg                        <  2:360302/148‑1 (MQ=255)                                                                                                                       gggTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCAt                                               >  3:344381/1‑123 (MQ=255)                                                                                                                       gggTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCAt                                               <  4:344381/123‑1 (MQ=255)                                                                                                                       gggTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                             >  5:86834/1‑125 (MQ=255)                                                                                                                       gggTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                             <  6:86834/125‑1 (MQ=255)                                                                                                                       gggTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGcc                     >  1:181190/1‑149 (MQ=255)                                                                                                                        ggTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGcc                     <  4:173159/148‑1 (MQ=255)                                                                                                                         gTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                             <  5:10246/123‑1 (MQ=255)                                                                                                                         gTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                             <  1:168819/123‑1 (MQ=255)                                                                                                                         gTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                             >  2:168819/1‑123 (MQ=255)                                                                                                                         gTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                             >  6:10246/1‑123 (MQ=255)                                                                                                                            gCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCaa                 <  3:221163/148‑1 (MQ=255)                                                                                                                            gCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCCGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAg                >  6:251160/1‑149 (MQ=255)                                                                                                                              tCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCa              >  7:356087/1‑149 (MQ=255)                                                                                                                               ccccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGcc                     <  4:342605/141‑1 (MQ=255)                                                                                                                               ccccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGcc                     >  3:342605/1‑141 (MQ=255)                                                                                                                                cccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCa              <  8:34394/147‑1 (MQ=255)                                                                                                                                 gaactcAACCCCGCAGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGGCCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCg           >  6:68127/6‑149 (MQ=255)                                                                                                                                 ccGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCg           <  5:43430/149‑1 (MQ=255)                                                                                                                                    ctcAACCCCGCAGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGgt        >  8:332959/3‑149 (MQ=255)                                                                                                                                    cccAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGgt        <  4:267311/149‑1 (MQ=255)                                                                                                                                     ccAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGgtg       <  6:212256/149‑1 (MQ=255)                                                                                                                                       aaCCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGtt     <  8:356087/149‑1 (MQ=255)                                                                                                                                        aCCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                             >  1:97949/1‑108 (MQ=255)                                                                                                                                        aCCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATgg                                             <  2:97949/108‑1 (MQ=255)                                                                                                                                          cccTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAg  <  4:179722/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                  CCTTGGTTGATTGCCAGGCGCAGCTGATCCATGGCCGCCAGGTACTGGTCGCAGGCCTCGGCCTGCCCCGGGTCGATCCACTGGATCTGGTCACCGCGAATGGTTTCGCGCACCTCCTGGGTAGCTCCCCGCCCAACCCCTGCGGGGTTGAGTTCGCCTTCGGCATCACGGCGCCGGCACTCGGCTGCAAGTGCGCGCACCAGGTCGGCAGGCAGAAAGTGCGCCTGCTGGGACCAGCCATGGGTGGCCAGATCGTCGACGACGGCCGCAAGCATCGGGTGTTCAG  >  NC_002947/5884840‑5885125

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