Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,867,228	G→A	6.5%	Q160* (CAG→TAG)	thyA ←	thymidylate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,867,228	0	G	A	6.5%	181.5 / 2.3	61	Q160* (CAG→TAG)	thyA	thymidylate synthase
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (2/2);  total (30/31)
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.18e-01

AGGGAAATTTCCCTGGTCTCGACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCG  >  NC_002947/5867080‑5867375                                                                                                                                                     |                                                                                                                                                    aGGGAAATTTCCCTGGTCTCGACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTg                                                                                                                                                     >  8:79641/1‑149 (MQ=255) aGGGAAATTTCCCTGGTCTCGACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTg                                                                                                                                                     >  4:291999/1‑149 (MQ=255)  gggAAATTTCCCTGGTCTCGACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGAGCGCAGTTCCAGCCGTGGAACAGTATGCGGCTGCTACCAGCGTCGTTGGAGATGGTGTCAATGCACTgg                                                                                                                                                    >  8:107146/1‑149 (MQ=255)       aTTTCCCTGGTCTCGACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTAc                                                                                                                                               <  7:82036/149‑1 (MQ=255)        tttCCCTGGTCTCGACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTAcc                                                                                                                                              <  7:79641/149‑1 (MQ=255)                 tctcGACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGt                                                                                                                                          <  1:317656/144‑1 (MQ=255)                 tctcGACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGt                                                                                                                                          >  2:317656/1‑144 (MQ=255)                 tctcGACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATc                                                                                                                                     <  7:31316/149‑1 (MQ=255)                   tcGACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGc                                                                                                                                   >  3:289573/1‑149 (MQ=255)                     gACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGcc                                                                                                                                  <  7:92800/148‑1 (MQ=255)                     gACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGc                                                                                                                                   <  3:291999/147‑1 (MQ=255)                     gACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTCCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTg                                                                                                                                                     <  7:137099/129‑1 (MQ=255)                     gACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTCCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTg                                                                                                                                                     >  8:137099/1‑129 (MQ=255)                              gTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGa                                                                                                                        <  7:107146/149‑1 (MQ=255)                              gTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGa                                                                                                                        >  4:115610/1‑149 (MQ=255)                                ggAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGAt                                                                                                                       >  6:161799/1‑148 (MQ=255)                                   aCTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACa                                                                                                                           <  3:41498/141‑1 (MQ=255)                                   aCTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACa                                                                                                                           >  4:41498/1‑141 (MQ=255)                                         aCAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGt                                                                                                                                          >  1:261440/1‑120 (MQ=255)                                         aCAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGt                                                                                                                                          <  2:261440/120‑1 (MQ=255)                                         aCAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAAcgc                                                                                                              >  7:282505/1‑148 (MQ=255)                                            gCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACg                                                                                                          >  6:62866/1‑149 (MQ=255)                                             cAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACa                                                                                                                           >  6:243784/1‑131 (MQ=255)                                             cAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACa                                                                                                                           <  5:243784/131‑1 (MQ=255)                                             cAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGc                                                                                                         >  3:301482/1‑149 (MQ=255)                                                gTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGt                                                                                                      <  5:161799/149‑1 (MQ=255)                                                    cacGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCctc                                                                                                  <  4:289573/149‑1 (MQ=255)                                                       ggTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTAcc                                                                                                                                              <  5:270444/102‑1 (MQ=255)                                                       ggTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTAcc                                                                                                                                              >  6:270444/1‑102 (MQ=255)                                                       ggTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTc                                                                                               <  8:282505/149‑1 (MQ=255)                                                           tcAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTAACCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTcc                                                                                           >  5:250998/2‑149 (MQ=255)                                                              gTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTcc                                                                                                    <  6:340890/137‑1 (MQ=255)                                                              gTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTcc                                                                                                    >  5:340890/1‑137 (MQ=255)                                                              gTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTgg                                                                                        <  3:115610/149‑1 (MQ=255)                                                                gCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTggg                                                                                       >  8:252544/1‑148 (MQ=255)                                                                              ttGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAAc                                                                        >  5:193670/1‑149 (MQ=255)                                                                                        ttCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATTGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTc                                                              >  3:162020/1‑149 (MQ=255)                                                                                          ccAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCTTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTc                                                              <  7:250047/147‑1 (MQ=255)                                                                                          ccAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCATTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTc                                                              >  8:250047/1‑147 (MQ=255)                                                                                            aGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCg                                                                                                       <  7:247860/104‑1 (MQ=255)                                                                                            aGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCg                                                                                                       >  8:247860/1‑104 (MQ=255)                                                                                             gCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCa                                                         >  3:41968/1‑149 (MQ=255)                                                                                                  ggAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTc                                                     >  7:66883/1‑148 (MQ=255)                                                                                                   gAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGa                                                   <  7:252544/149‑1 (MQ=255)                                                                                                    aaCAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGAt                                                  >  6:322529/1‑149 (MQ=255)                                                                                                    aaCAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGAt                                                  <  4:162020/149‑1 (MQ=255)                                                                                                     aCAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGt                                                                                                                              <  2:112298/72‑1 (MQ=255)                                                                                                         gATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGgccg                                             >  3:324413/1‑149 (MQ=255)                                                                                                              ggcggcTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCgg                                                            >  4:189149/1‑129 (MQ=255)                                                                                                              ggcggcTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCgg                                                            <  3:189149/129‑1 (MQ=255)                                                                                                                      cccGGGTCGTTGGCGATGGTGTCCAGGCACTAGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAg                                >  8:221221/1‑149 (MQ=255)                                                                                                                      cccGGGTCGTTGGCGATGGTGTCCAGGCACTAGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAg                                >  8:219340/1‑149 (MQ=255)                                                                                                                           gTCGTTGGCGATGGTGTCCAGGCACTAGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGa                           <  7:219340/149‑1 (MQ=255)                                                                                                                           gTCGTTGGCGATGGTGTCCAGGCACTAGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGa                           <  7:221221/149‑1 (MQ=255)                                                                                                                               ttGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATTCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATgcg                       <  6:250998/149‑1 (MQ=255)                                                                                                                                 ggCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATgcgc                      <  5:62866/148‑1 (MQ=255)                                                                                                                                  gCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCtt                    >  2:34664/1‑149 (MQ=255)                                                                                                                                  gCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCtt                    <  4:41968/149‑1 (MQ=255)                                                                                                                                   cGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTg                   <  5:322529/149‑1 (MQ=255)                                                                                                                                   cGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTg                   <  4:324413/149‑1 (MQ=255)                                                                                                                                   cGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTg                   >  3:47980/1‑149 (MQ=255)                                                                                                                                              aGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCa        >  1:179782/1‑149 (MQ=255)                                                                                                                                                 cACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTg     <  8:66883/149‑1 (MQ=255)                                                                                                                                                  aCTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTgg    >  8:35347/1‑149 (MQ=255)                                                                                                                                                   cTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGgc   <  4:301482/149‑1 (MQ=255)                                                                                                                                                    tGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGgcg  >  1:299765/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AGGGAAATTTCCCTGGTCTCGACATTCGGGTGGAACTGGTACAGCAAGTGACACGGTGGCAGTGCCATTTCGTCCAGTTGCGCGCAGTTCCAGCCGTGGAACAGGATGCGGCGGCTACCCGGGTCGTTGGCGATGGTGTCCAGGCACTGGCGTACCTGGTCGATCGCCTTGTACAGGATGACGAACGCGACGCCGTCCTCTTCGTCCTGGGCAATGCGGCGGAAACCGGCTTGCTCGGCCATTTCGATGGCCGCCGGGTTGCTCAGCGGGATGCGCTTGTAGCCCGGCCACTGGCG  >  NC_002947/5867080‑5867375

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