Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,230,033	G→C	6.1%	P61R (CCG→CGG)	PP_4607 ←	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,230,033	0	G	C	6.1%	209.3 / 3.6	65	P61R (CCG→CGG)	PP_4607	transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base C (2/2);  total (35/30)
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.51e-01

ACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCAGCCTGGAA  >  NC_002947/5229895‑5230181                                                                                                                                           |                                                                                                                                                     aCGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCCCGGCGgcgcg                                                                                                                                            <  5:18025/149‑2 (MQ=255)    cccACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCtttt                                                                                                                                         <  7:264252/149‑1 (MQ=255)     ccACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTc                                                                                                                                        >  5:359321/1‑149 (MQ=255)     ccACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCCCCGgcggcgcggctggcgga                                                                                                                                        >  2:326479/1‑142 (MQ=255)     ccACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCCCCGgcggcgcggctggcgga                                                                                                                                        <  1:326479/149‑8 (MQ=255)            cAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCaaa                                                                                                                                 <  6:359321/149‑1 (MQ=255)            cAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCa                                                                                                                                   >  7:122324/1‑147 (MQ=255)               cACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCa                                                                                                                              >  2:287291/1‑149 (MQ=255)               cACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCa                                                                                                                              >  3:249444/1‑149 (MQ=255)                   aggcaggcTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGc                                                                                                                               <  2:330307/144‑1 (MQ=255)                   aggcaggcTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGc                                                                                                                               >  1:330307/1‑144 (MQ=255)                     gcaggcTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTggg                                                                                                                        >  6:175530/1‑149 (MQ=255)                       aggcTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGcc                                                                                                                      >  1:130097/1‑149 (MQ=255)                             ccAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTc                                                                                                                 <  5:200660/148‑1 (MQ=255)                               aaCGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTc                                                                                                                 <  4:136798/146‑1 (MQ=255)                               aaCGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTc                                                                                                                 >  3:136798/1‑146 (MQ=255)                                  gCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGc                                                                                                           <  2:130097/149‑1 (MQ=255)                                  gCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGc                                                                                                           >  4:124036/1‑149 (MQ=255)                                    cAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGt                                                                                                               <  4:96068/143‑1 (MQ=255)                                    cAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGt                                                                                                               >  3:96068/1‑143 (MQ=255)                                            ccACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGgt                                                                                                   <  8:122324/147‑1 (MQ=255)                                             cACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTgg                                                                                                                         >  1:316573/1‑124 (MQ=255)                                             cACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTgg                                                                                                                         <  2:316573/124‑1 (MQ=255)                                                     gccagccaCCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCg                                                                                        >  1:132987/1‑149 (MQ=255)                                                      ccagccaCCACGCGCCGGTGGTCTGGCCATGGCATAGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCATTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCgg                                                                                       >  5:52409/1‑149 (MQ=255)                                                       cagccaCCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGc                                                                                      <  2:124934/149‑1 (MQ=255)                                                       cagccaCCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGc                                                                                      >  3:244842/1‑149 (MQ=255)                                                                    ccGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAgcgc                                                                         <  2:132987/149‑1 (MQ=255)                                                                       gtggtCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTg                                                                      <  1:287291/149‑1 (MQ=255)                                                                       agatcttGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCTGTGTGCGGGGTCGGCCAGCAACCAGCGCtt                                                                       <  1:283429/142‑1 (MQ=255)                                                                          gtCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAAt                                                                   >  5:38973/1‑149 (MQ=255)                                                                            cTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCaaa                                                                                                                                 >  2:318184/1‑85 (MQ=255)                                                                            cTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCaaa                                                                                                                                 <  1:318184/85‑1 (MQ=255)                                                                            cTGGCAATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGAGAGGTGTGCGGGGTCGGCCAGAAACCAGCGCTTGAATgc                                                                 >  4:328887/1‑149 (MQ=255)                                                                              ggCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGACAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCATCGCtt                                                                       >  2:283429/1‑141 (MQ=255)                                                                               gCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgcc                                                                <  5:175530/147‑1 (MQ=255)                                                                                     gCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCa                                                                             >  4:192070/1‑128 (MQ=255)                                                                                     gCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCtt                                                         <  4:249444/148‑1 (MQ=255)                                                                                       aTCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCgg                                                      >  1:261603/1‑149 (MQ=255)                                                                                               aGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCg                                              >  8:191236/1‑149 (MQ=255)                                                                                               aGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCAAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCg                                              >  4:347464/1‑149 (MQ=255)                                                                                                  gcTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGt                                           <  3:124036/149‑1 (MQ=255)                                                                                                    tGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGt                                           <  7:78736/147‑1 (MQ=255)                                                                                                    tGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGt                                           >  8:78736/1‑147 (MQ=255)                                                                                                     gCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGc                                        <  4:244842/149‑1 (MQ=255)                                                                                                      cAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCt                                       >  3:36926/1‑149 (MQ=255)                                                                                                          gCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGtt                                   <  6:38973/149‑1 (MQ=255)                                                                                                                 gCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCgg                            >  3:253549/1‑149 (MQ=255)                                                                                                                  ctctTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGt                           >  8:154159/1‑149 (MQ=255)                                                                                                                  ctctTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGt                           >  2:29340/1‑149 (MQ=255)                                                                                                                  ctctTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGt                           <  3:347464/149‑1 (MQ=255)                                                                                                                    ctTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGtt                          >  7:190902/1‑148 (MQ=255)                                                                                                                       cGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGaa                       <  8:190902/148‑1 (MQ=255)                                                                                                                       cGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGaa                       <  7:191236/148‑1 (MQ=255)                                                                                                                        gTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAAc                      >  5:248477/1‑148 (MQ=255)                                                                                                                        gTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAAc                      >  1:315698/1‑148 (MQ=255)                                                                                                                            gCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTg                 <  1:29340/149‑1 (MQ=255)                                                                                                                            gCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTg                 >  5:118043/1‑149 (MQ=255)                                                                                                                             ccAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTg                 >  1:17914/1‑148 (MQ=255)                                                                                                                             ccAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTg                 >  4:186905/1‑148 (MQ=255)                                                                                                                             ccAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGc                <  2:17914/149‑1 (MQ=255)                                                                                                                               aaaagcgccgcCCCGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAg              >  2:192826/5‑149 (MQ=255)                                                                                                                                 ccgcgccgcCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAggg            <  4:36926/149‑1 (MQ=255)                                                                                                                                   gcgccgcCGTGGCGGCGCTTTGCTGCCTAAGAAGCTGGGACTGGTCTTAGGCTGCGCGGTGTGCGTTTTCGTCCAGCAACCAGCGCTTGATTGCCGCACCTTCGGCTTGCCATGGTTGCTCGTTGATCCGTTTTAACCATTTCAGTGCa          <  6:52409/149‑1 (MQ=255)                                                                                                                                   gcgccgcCCCGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAggg            >  1:25936/1‑147 (MQ=255)                                                                                                                                   gcgccgcCCCGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAggg            <  2:25936/147‑1 (MQ=255)                                                                                                                                     gccgccGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCAGc        <  4:253549/149‑1 (MQ=255)                                                                                                                                       cgccgGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCAGCCt      <  2:261603/149‑1 (MQ=255)                                                                                                                                           ggggCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCAGCCTGGaa  >  4:274655/1‑149 (MQ=255)                                                                                                                                           ggggAGGCGCTTTTCTGCAAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCTGCCAGCAACCAGAGCTTTAATGCCGCACCTTCGGCTTGCACGGGTTTCTCGTTGATCCGTTTGAACCATTGCAGGGCAGCATGGaa  >  5:42608/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                     ACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCAGCCTGGAA  >  NC_002947/5229895‑5230181

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