Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,033,545	A→T	8.8%	D29V (GAC→GTC)	PP_2646 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,033,545	0	A	T	8.8%	160.6 / 4.7	57	D29V (GAC→GTC)	PP_2646	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base T (3/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

GCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACGC  >  NC_002947/3033406‑3033693                                                                                                                                            |                                                                                                                                                       gCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcg                                                                                                                                               >  2:61691/1‑149 (MQ=255)  cccAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgc                                                                                                                                              <  8:279799/149‑1 (MQ=255)   ccAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgc                                                                                                                                              >  8:84111/1‑148 (MQ=255)      ggCCATGAAGACCCTCTATATCCGGCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGTGAGCAGAACATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGTCGCTGCCGCCAGGCCCTGGACGACGTccagcgcc‑agc                                                                                                                                         <  8:214911/149‑1 (MQ=255)          aTGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCGCTGGACGTCGTccagggcctggcg                                                                                                                                        >  4:2804/1‑142 (MQ=255)          aTGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCGCTGGACGTCGTccagggcctggcg                                                                                                                                        <  3:2804/147‑6 (MQ=255)               gACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATcc                                                                                                                                >  3:57646/1‑149 (MQ=255)                  ccTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCggg                                                                                                                             >  6:130650/1‑149 (MQ=255)                    tctATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGgtg                                                                                                                           >  6:103980/1‑149 (MQ=255)                    tctATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGATCGGCCG‑CTGGCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGgtg                                                                                                                           >  2:249826/1‑149 (MQ=21)                     ctATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGgtgt                                                                                                                          >  3:185895/1‑149 (MQ=255)                      tataTCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGgtg                                                                                                                           <  5:103980/147‑1 (MQ=255)                       atatCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGtt                                                                                                                         >  8:275064/1‑148 (MQ=255)                        tatCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGgtgt                                                                                                                          <  4:185895/146‑1 (MQ=255)                        tatCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGATCGGCCG‑CTGGCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGgtgtcg                                                                                                                        <  1:249826/148‑3 (MQ=18)                         atCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGa                                                                                                                       >  8:31470/1‑148 (MQ=255)                            cgccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCa                                                                                                                   >  7:274985/1‑149 (MQ=255)                              ccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACg                                                                                                                 >  7:120894/1‑149 (MQ=255)                              ccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACg                                                                                                                 <  4:57646/149‑1 (MQ=255)                                 tttCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGa                                                                                                                       >  3:319723/1‑140 (MQ=255)                                 tttCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGa                                                                                                                       <  4:319723/140‑1 (MQ=255)                                        cAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCa                                                                                                       <  3:340964/149‑1 (MQ=255)                                                cTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCtt                                                                                                <  1:61691/148‑1 (MQ=255)                                                  gAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCa                                                                                                           <  5:261533/135‑1 (MQ=255)                                                  gAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCa                                                                                                           >  6:261533/1‑135 (MQ=255)                                                   aGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                               <  7:275064/146‑1 (MQ=255)                                                   aGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                               <  7:31470/146‑1 (MQ=255)                                                   aGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                               <  8:120894/146‑1 (MQ=255)                                                   aGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                               <  5:32423/146‑1 (MQ=255)                                                          cAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGaaa                                                                                      <  5:130650/148‑1 (MQ=255)                                                           aaGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATc                                                                                    >  7:106809/1‑149 (MQ=255)                                                           aaGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATc                                                                                    >  6:302082/1‑149 (MQ=255)                                                                gagCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAg                                                                               >  6:344931/1‑149 (MQ=255)                                                                   cAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtggtgg                                                                                           <  6:142638/134‑1 (MQ=255)                                                                   cAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtggtgg                                                                                           >  5:142638/1‑134 (MQ=255)                                                                   cAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTg                                                                            >  7:235679/1‑149 (MQ=255)                                                                    aGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTgg                                                                           <  5:302082/149‑1 (MQ=255)                                                                      aCCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTggg                                                                          >  1:189082/1‑148 (MQ=255)                                                                      aCCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTggg                                                                          <  2:189082/148‑1 (MQ=255)                                                                        cATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGg‑cgc                                                                      >  8:278771/1‑149 (MQ=255)                                                                         aTGCGAGTACTGATATGTGCAGGGCGCAACTACGACGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGg‑cgcg                                                                     >  7:300484/1‑149 (MQ=255)                                                                          tGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAg                                                                               >  1:93929/1‑139 (MQ=255)                                                                          tGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAg                                                                               <  2:93929/139‑1 (MQ=255)                                                                                     atatGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGgcgc                                                         <  5:344931/149‑1 (MQ=255)                                                                                          tgCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATc                                                    <  4:287302/149‑1 (MQ=255)                                                                                            cAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATCTg                                                  <  7:278771/149‑1 (MQ=255)                                                                                            cAGGGCGCAACTACGACGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATCt                                                   <  8:300484/148‑1 (MQ=255)                                                                                               ggCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATCTGGTg                                               >  2:270349/1‑149 (MQ=255)                                                                                               ggCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATCTGGTg                                               <  8:274985/149‑1 (MQ=255)                                                                                                         cGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGc                                     >  7:302210/1‑149 (MQ=255)                                                                                                              acacCCGGATCGGCCG‑CTGGCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAAc                                >  2:59089/1‑149 (MQ=18)                                                                                                                  ccGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCAGGGGGGCAGCCAGTATCTTGGTgggg                                                                                         >  4:268049/1‑89 (MQ=255)                                                                                                                    ggCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATc                                                    <  1:254796/123‑1 (MQ=255)                                                                                                                    ggCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTAGATCAACGTTGTCACCCAATAAATTAAAAAAAAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATc                                                    >  2:254796/1‑123 (MQ=255)                                                                                                                              ccAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGACGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGc                 <  1:270349/149‑1 (MQ=255)                                                                                                                              ccAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCa                <  8:106809/149‑1 (MQ=255)                                                                                                                                     cTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGc         <  8:235679/149‑1 (MQ=255)                                                                                                                                            acgTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACt                               >  5:92394/1‑120 (MQ=255)                                                                                                                                            acgTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACt                               <  6:92394/120‑1 (MQ=255)                                                                                                                                            acgTCCAGCGCC‑AGCGGCCGATCCGGGTGTAGATCAACGTTGTCACCCTATAAATAAAAAAAAAAATCGAAAGATGGG‑AGCGCGAGCATGGCGCCGATCTGGTGCGTTAACCGCCAAACTGTCAGCTGCAAGGCAAGCAGGCCGAACGc  >  8:177257/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                       GCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGG‑CGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGGCCGAACGC  >  NC_002947/3033406‑3033693

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑