Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,611,430	T→C	11.6%	F142L (TTC→CTC)	PP_4931 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,611,430	0	T	C	11.6%	151.7 / 10.4	61	F142L (TTC→CTC)	PP_4931	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/24);  new base C (3/4);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

TACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCTTGCGTGAAGGTAATGCGTTGGTCGGCCTGCGCG  >  NC_002947/5611282‑5611577                                                                                                                                                     |                                                                                                                                                    tACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCCCTGCAGGACCGCGCGt                                                                                                                                                     >  2:300820/1‑149 (MQ=255)    ccGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGTACTGCAGGACCGCGCGTTc                                                                                                                                                   >  8:246254/1‑148 (MQ=255)     cGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGCTCaa                                                                                                                                                 <  4:248675/149‑1 (MQ=255)            tctcCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCCTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAg                                                                                                                                          >  8:103296/1‑149 (MQ=255)              tcCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGCTCAAGGGCAAGgt                                                                                                                                        >  5:124546/1‑149 (MQ=255)                cGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAgg                                                                                                                                         <  1:137064/146‑1 (MQ=255)                      ccACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGCTCAAGGGCAAGGTGTACCGcc                                                                                                                                >  1:125331/1‑149 (MQ=255)                       cACTACCTTTGGTCGCCGGGCACCTTGGCCGGCAAACTCGCAAGCTTATTCGCCAGAAAGGCCGTGCGCGGCCGGGTCGACCATGTCAAGGGTGACCAGTTGCCACCGGCAATGCAGGACCGCGCTGTCAAGGTAAAGGTGTACCGCCt                                                                                                                               >  4:49448/1‑149 (MQ=255)                       cACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCt                                                                                                                               >  6:106649/1‑149 (MQ=255)                              tgtgGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTAc                                                                                                                                    <  3:213562/137‑1 (MQ=255)                              tgtgGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTAc                                                                                                                                    >  4:213562/1‑137 (MQ=255)                                  gTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACtgg                                                                                                                    <  4:46186/149‑1 (MQ=255)                                  gTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACtgg                                                                                                                    >  6:193620/1‑149 (MQ=255)                                   tCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACtggt                                                                                                                   <  6:288156/149‑1 (MQ=255)                                          gCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGAt                                                                                                            >  7:253610/1‑149 (MQ=255)                                            aCCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGCTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGt                                                                                                          <  3:74698/149‑1 (MQ=255)                                             ccTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGCTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTc                                                                                                         <  6:124546/149‑1 (MQ=255)                                                    cGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGcc                                                                                                  >  7:255502/1‑149 (MQ=255)                                                        aaCCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACtgg                                                                                                                    >  6:175707/1‑127 (MQ=255)                                                        aaCCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACtgg                                                                                                                    <  5:175707/127‑1 (MQ=255)                                                        aaCCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCctgct                                                                                              >  3:236094/1‑149 (MQ=255)                                                              acaAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGCTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGc                                                                                           >  5:237406/1‑146 (MQ=255)                                                              acaAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGCTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGc                                                                                           <  6:237406/146‑1 (MQ=255)                                                                aaGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCaa                                                                                                                                                 >  5:192741/1‑90 (MQ=255)                                                                aaGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCaa                                                                                                                                                 <  6:192741/90‑1 (MQ=255)                                                                aaGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAAcc                                                                                      >  8:66739/1‑149 (MQ=255)                                                                  gCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTg                                                                                    <  5:106649/149‑1 (MQ=255)                                                                   cttcttCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTgg                                                                                   >  3:213894/1‑149 (MQ=255)                                                                             aGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGc                                                                         <  5:330039/149‑1 (MQ=255)                                                                                 aGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGc                                                                     <  7:103296/149‑1 (MQ=255)                                                                                       tgcagtgccggtggcaactggtcacccttgacctggtcgacccggcACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGa                                                                                                     >  1:273183‑M2/47‑111 (MQ=255)                                                                                       tgcagtgccggtggcaactggtcacccttgacctggtcgacccggcACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGa                                                                                                     <  2:273183‑M2/65‑1 (MQ=255)                                                                                        gcgcGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCt                                                              <  4:236094/149‑1 (MQ=255)                                                                                        gcagtgccggtggcaactggtcacccttgacctggtcgacccggcACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAgtcgt                                                              >  5:286490‑M2/46‑147 (MQ=255)                                                                                           cgGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGTACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGc                                                           <  7:246254/149‑1 (MQ=255)                                                                                           cgGCCGGGTCGACCAGGTCAAGGGTGACCAGGTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGGCGAGGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGc                                                           <  3:49448/149‑1 (MQ=255)                                                                                            ggCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCt                                                          >  8:282231/1‑149 (MQ=255)                                                                                              CCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTgg                                                        >  6:331356‑M1/1‑149 (MQ=255)                                                                                               CGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGc                                                       >  8:253610‑M1/1‑149 (MQ=255)                                                                                                 GGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAg                                                     <  5:212104‑M1/149‑1 (MQ=255)                                                                                                     gACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAg                                                     >  1:151243/1‑145 (MQ=255)                                                                                                     gACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAg                                                     <  2:151243/145‑1 (MQ=255)                                                                                                      aCCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGaa                                                <  7:282231/149‑1 (MQ=255)                                                                                                       ccAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAAc                                               >  2:48788/1‑149 (MQ=255)                                                                                                           gTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCAt                                           >  8:214427/1‑149 (MQ=255)                                                                                                              aaGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGa                                                                   >  6:165379/1‑122 (MQ=255)                                                                                                              aaGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGa                                                                   <  5:165379/122‑1 (MQ=255)                                                                                                                   tGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCt                                   <  7:66739/149‑1 (MQ=255)                                                                                                                     aCCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCtt                                  <  5:180114/148‑1 (MQ=255)                                                                                                                     aCCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCtt                                  <  1:48788/148‑1 (MQ=255)                                                                                                                         gTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCTTGCGt                              >  2:94814/1‑148 (MQ=255)                                                                                                                         gTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCTTGCGt                              >  1:175333/1‑148 (MQ=255)                                                                                                                          ttGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCTTGCGTg                             >  8:269247/1‑148 (MQ=255)                                                                                                                          gtcgacccggcACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCTTGCGTGa                            <  6:286490‑M2/138‑1 (MQ=255)                                                                                                                            gCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCTTGCGTGaa                           <  5:193620/148‑1 (MQ=255)                                                                                                                                 cggcACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCTTGCGTGAAGGTAAt                     <  4:213894‑M2/145‑1 (MQ=255)                                                                                                                                                gcgcgTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATTGAACCCTTGCGTGAAGGTAATGCGTTGGTCGGCCTg      >  7:132466/1‑149 (MQ=255)                                                                                                                                                gcgcgTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCTTGCGTGaa                           >  4:77911/1‑128 (MQ=255)                                                                                                                                                gcgcgTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCTTGCGTGaa                           <  3:77911/128‑1 (MQ=255)                                                                                                                                                gcgcgTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATCTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCTAGCGTGAAGGTAATGCGTTGGTCGGCCTg      >  7:276005/1‑149 (MQ=255)                                                                                                                                                gcgagTTCAAGGGAAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATCTCCCGATCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGACAGCATAGAACCCTTGCTTTAAGGTAATGCGTTGGTCGGCCTg      >  7:280450/1‑149 (MQ=255)                                                                                                                                                    gTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCTTGCGTGAAGGTAATGCGTTGGTCGGCCTGcgcg  >  4:16889/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGACCAGGTCAAGGGTGACCAGTTGCCACCGGCACTGCAGGACCGCGCGTTCAAGGGCAAGGTGTACCGCCTGGCGGAACTGGTGGTCGATGTCCCGAGCCTGCTGGCCAACCTGGCCGAGCTGGCCGGCGAGTCCTTGCTGGCAGGTGAACGCATCGAACCCTTGCGTGAAGGTAATGCGTTGGTCGGCCTGCGCG  >  NC_002947/5611282‑5611577

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