Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,363,695	A→T	5.6%	V61D (GTC→GAC)	PP_2076 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,363,695	0	A	T	5.6%	198.4 / 2.1	71	V61D (GTC→GAC)	PP_2076	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/37);  new base T (2/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

ATGGCGCTTGATCAGCAGGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGTTTCTTTCAACGAGTCCATGTATAGCTCCAGCGAG  >  NC_002947/2363550‑2363841                                                                                                                                                   |                                                                                                                                                     aTGGCGCTTGATCAGCAGGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGAcgg                                                                                                                                                   >  6:56263/1‑149 (MQ=34)    gcgcTTGATCAGCAGGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGGGTAGCTGGAATGGGGCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGGCCGGCAATTTGGG‑GAGGCGGGCACGTTGTGCGGGGTTGCCGCCGAcggcg                                                                                                                                                 <  5:109227/148‑1 (MQ=11)    gcgcTTGATCAGCAGGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGGCCAGGCGGGCACGTTGTGCCGGGTTGCCGCCGAcggcg                                                                                                                                                 <  1:66674/149‑1 (MQ=21)    gcgcTTGATCAGCAGGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGAcggcgg                                                                                                                                                <  4:129754/149‑1 (MQ=34)        ttGATCAGCAGGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCGCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGGGCCGGGGTGCCGCCGACGGCGGCGTa                                                                                                                                            <  6:374514/149‑1 (MQ=255)        atGATCAGCAGGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTa                                                                                                                                            <  4:62366/148‑1 (MQ=21)          gATCAGCAGGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGTCggcggcaacc                                                                                                                                           <  6:299275/148‑5 (MQ=21)           aTCAGCAGGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAgc                                                                                                                                          <  1:347117/148‑1 (MQ=34)            tCAGCAGGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGAcg                                                                                                                                                    >  8:155190/1‑137 (MQ=34)            tCAGCAGGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGAcg                                                                                                                                                    <  7:155190/137‑1 (MQ=34)               gcagGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGTTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAgg                                                                                                                                     <  3:268261/149‑1 (MQ=21)                   gTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAgggg                                                                                                                                   >  5:212877/1‑147 (MQ=34)                   gTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAgggg                                                                                                                                   >  1:331279/1‑147 (MQ=34)                   gTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAgggg                                                                                                                                   >  3:250519/1‑147 (MQ=34)                    tAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGgcag                                                                                                                                <  2:38485/149‑1 (MQ=34)                        gaAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGt                                                                                                                            >  5:140064/1‑149 (MQ=34)                        gaAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGt                                                                                                                            <  6:140064/149‑1 (MQ=34)                         aaaGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAgg                                                                                                                                     <  1:130404/139‑1 (MQ=34)                         aaaGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAgg                                                                                                                                     >  2:130404/1‑139 (MQ=34)                             cccTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGaa                                                                                                                       <  8:125031/149‑1 (MQ=34)                               cTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑gc                                                                                                                    >  3:108419/1‑149 (MQ=34)                                    gcAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑Gc‑ggcgg                                                                                                              <  5:56263/149‑1 (MQ=34)                                    gcAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑Gc‑ggcgg                                                                                                              <  8:128841/149‑1 (MQ=34)                                    gcAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑Gc‑ggcgg                                                                                                              <  6:212877/149‑1 (MQ=34)                                       aaGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTg                                                                                                           <  7:297730/149‑1 (MQ=34)                                           gCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGt                                                                                                       >  7:323658/1‑149 (MQ=34)                                            cGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTc                                                                                                      >  2:140710/1‑149 (MQ=34)                                               aTGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAAGGCAGGCGGGTGCAGTCg                                                                                                     >  1:311853/1‑149 (MQ=255)                                               aTGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGc                                                                                                   <  7:236484/149‑1 (MQ=34)                                                  gCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑Gc‑ggcgg                                                                                                              >  2:185914/1‑135 (MQ=34)                                                  gCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑Gc‑ggcgg                                                                                                              <  1:185914/135‑1 (MQ=34)                                                    gTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACgg                                                                                              >  7:7431/1‑149 (MQ=34)                                                    gTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACgg                                                                                              <  4:302367/149‑1 (MQ=34)                                                             aTGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCGGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTAc                                                                                     <  7:162033/149‑1 (MQ=25)                                                             aTGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGCTGCAGTCGGCGACGGGGCTGTTAc                                                                                     >  4:117783/1‑149 (MQ=25)                                                                  tCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGgc                                                                                                                                  <  4:234994/101‑1 (MQ=35)                                                                  tCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGgc                                                                                                                                  >  3:234994/1‑101 (MQ=35)                                                                     gCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGcc                                                                             >  6:78269/1‑149 (MQ=37)                                                                         gCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTAc                                                                                     >  1:277834/1‑137 (MQ=37)                                                                         gCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTAc                                                                                     <  2:277834/137‑1 (MQ=37)                                                                         gCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGc                                                                         >  1:181301/1‑149 (MQ=37)                                                                         gCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGc                                                                         >  7:286169/1‑149 (MQ=37)                                                                                 ccTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGCTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCCGGCGGTAGccc                                                                 >  8:335540/1‑148 (MQ=21)                                                                                   tGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCg                                                                 >  5:149248/1‑147 (MQ=34)                                                                                     gAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAg                                                              <  5:60150/148‑1 (MQ=34)                                                                                             tGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATgcgc                                                     >  1:79320/1‑149 (MQ=34)                                                                                              gagaTAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGcc                                                    >  3:242140/1‑149 (MQ=34)                                                                                               agaTAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCt                                                   <  2:331279/149‑1 (MQ=34)                                                                                               agaTAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCt                                                   <  8:286169/149‑1 (MQ=34)                                                                                               agaTAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCt                                                   <  4:250519/149‑1 (MQ=34)                                                                                                gaTAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTg                                                  >  1:248303/1‑149 (MQ=34)                                                                                                 aTAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGTCGGGCGGTAGCCGGAGGTGATGCGCCTg                                                  >  6:231058/1‑148 (MQ=21)                                                                                                 aTAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGc                                                 <  8:323658/149‑1 (MQ=34)                                                                                                  tAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGca                                                <  4:108419/149‑1 (MQ=34)                                                                                                            aTTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAgg                                      >  8:21734/1‑149 (MQ=34)                                                                                                            aTTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAgg                                      <  2:248303/149‑1 (MQ=34)                                                                                                                 gg‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGGAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGtttc                                 <  3:117783/149‑1 (MQ=21)                                                                                                                  g‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGtttct                                >  5:315227/1‑149 (MQ=35)                                                                                                                         gggCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGTCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGTTTCTTTCAAc                          <  5:231058/149‑1 (MQ=255)                                                                                                                           gCAGGTTGTGCCGGGTGGCGGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGTGTGCAGTCGTCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGTTGATGCGCCTGCACAGTTCGAGGGTTTCTTTCAACGa                        <  7:335540/149‑1 (MQ=14)                                                                                                                             aCGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGTTTCTTTCAACGAGt                      >  1:262892/1‑149 (MQ=255)                                                                                                                               gTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGTTTCTTTCAACGAGTcc                    <  2:79320/149‑1 (MQ=255)                                                                                                                                  gtgCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGTTTCTTTCAACGAGTCCAt                  <  2:181301/148‑1 (MQ=255)                                                                                                                                   tgCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGTTTCTTTCAACGAGTCCATGt                <  5:78269/149‑1 (MQ=255)                                                                                                                                     gcGGGTGGCCGCCGACGGCGGCGTAGCGCAGGCGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGGTACGCAGGGGGGGGCGGTAGCCGGAGGTTATGCGCCTGCACAGTTCGAGGGTTTCTTTAAACGAGTCCATGtat              <  4:242140/148‑1 (MQ=14)                                                                                                                                       gctacgccgccGTCGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAgg                                      <  2:19972/118‑1 (MQ=255)                                                                                                                                       gctacgccgccGTCGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAgg                                      >  1:19972/6‑123 (MQ=255)                                                                                                                                       gctacgccgccGTCGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCCGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGTTTCTTTCAACGAGTCCATGTATAg            >  6:29917/6‑149 (MQ=255)                                                                                                                                        ggTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGTTTCTTTCAACGAGTCCATGTATAGc           >  5:105820/1‑149 (MQ=255)                                                                                                                                           tGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGTTTCTTTCAACGAGTCCATGTATAGCTcc        <  2:121984/149‑1 (MQ=255)                                                                                                                                                  gACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGTTTCTTTCAACGAGTCCATGTATAGCTCCAGCGAg  >  8:10813/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                     ATGGCGCTTGATCAGCAGGTAGAGAAAGCCCTCGCGCAAAGGGCGGATGGCGTAGCTGGAATGGGTCAGCTCGCCGACATCCTGGAACGGGCTGAGATAGCCCGGCAATTTGGG‑CAGGCGGGCACGTTGTGCCGGGTTGCCGCCGACGGCGGCGTAGCGCAGGGGCAGCAGTTCGAA‑GC‑GGCGGGTGCAGTCGGCGACGGGGCTGTTACGCAGGGCCGGGCGGTAGCCGGAGGTGATGCGCCTGCACAGTTCGAGGGTTTCTTTCAACGAGTCCATGTATAGCTCCAGCGAG  >  NC_002947/2363550‑2363841

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