Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,497,329	Δ1 bp	8.4%	coding (268/765 nt)	PP_5562 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,497,329	0	C	.	8.4%	220.3 / 13.6	71	coding (268/765 nt)	PP_5562	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base . (2/4);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

CCAGGTCGAAGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGC  >  NC_002947/3497185‑3497471                                                                                                                                                 |                                                                                                                                                 ccAGGTCGAAGGCGTGGAAGCCCAGCACCCAGAGGGTCGGGCTGGTGATGCGGATGACCTGGGCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGTGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTg                                                                                                                                              <  4:320218/149‑1 (MQ=255)       cGAAGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGc                                                                                                                                       >  6:234177/1‑149 (MQ=255)        gAAGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGc                                                                                                                                       <  5:234177/148‑1 (MQ=255)          aGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCT‑GGTCACCGAGCTgg                                                                                                                                    <  1:20852/149‑1 (MQ=255)          aGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCT‑GGTCACCGAGCTgg                                                                                                                                    <  5:76967/149‑1 (MQ=255)               tGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACgg                                                                                                                               >  8:304810/1‑149 (MQ=255)                ggAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGgc                                                                                                                              >  3:334566/1‑149 (MQ=255)                  aaGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGgcgc                                                                                                                            <  7:304810/149‑1 (MQ=255)                       caccacCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgcc                                                                                                                        >  2:153171/1‑148 (MQ=255)                        accaccCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgccga                                                                                                                      <  6:106164/149‑1 (MQ=255)                        accaccCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgcc                                                                                                                        >  7:293204/1‑147 (MQ=255)                        accaccCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgcc                                                                                                                        >  7:33712/1‑147 (MQ=255)                          caccCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgccgagc                                                                                                                    >  5:241177/1‑149 (MQ=255)                            cccAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgccgagccg                                                                                                                  <  5:1325/149‑1 (MQ=255)                            cccAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgccgagcc                                                                                                                   >  5:38081/1‑148 (MQ=255)                              cAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgccgagccga                                                                                                                 >  6:315338/1‑148 (MQ=255)                               aGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATa                                                                                                               >  4:278367/1‑149 (MQ=255)                                       gggCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTa                                                                                                       <  6:235173/149‑1 (MQ=255)                                       gggCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTa                                                                                                       <  5:150878/149‑1 (MQ=255)                                       gggCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTa                                                                                                       <  3:98434/149‑1 (MQ=255)                                             gTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTc                                                                                                 <  6:241177/149‑1 (MQ=255)                                                aTGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGgc                                                                                                                              >  6:25071/1‑117 (MQ=255)                                                aTGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGgc                                                                                                                              <  5:25071/117‑1 (MQ=255)                                                aTGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGa                                                                                               <  8:33712/148‑1 (MQ=255)                                                    ggATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCTGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAgg                                                                                          <  5:315338/149‑1 (MQ=255)                                                    ggATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAgg                                                                                          <  3:278367/149‑1 (MQ=255)                                                         aCCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTg                                                                                     <  4:334566/149‑1 (MQ=255)                                                              gCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTa                                                                                                       <  3:38567/126‑1 (MQ=255)                                                              gCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTa                                                                                                       >  4:38567/1‑126 (MQ=255)                                                              gCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGt                                                                                >  3:285389/1‑149 (MQ=255)                                                                cGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTgg                                                                              >  3:8680/1‑149 (MQ=255)                                                                  cTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTGGAGTTCGACCAGGGCCTGGAAGTGGcc                                                                            >  4:347175/1‑149 (MQ=255)                                                                       ttCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTggg                                                                       >  4:249400/1‑149 (MQ=255)                                                                       ttCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTgg                                                                        <  1:153171/148‑1 (MQ=255)                                                                       ttCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTgg                                                                        <  4:285389/148‑1 (MQ=255)                                                                            gggTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGAcc                                                                                             >  7:116581/1‑122 (MQ=255)                                                                            gggTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGAcc                                                                                             <  8:116581/122‑1 (MQ=255)                                                                               tCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑c                                                                  <  1:39534/145‑1 (MQ=255)                                                                               tCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑c                                                                  >  2:39534/1‑145 (MQ=255)                                                                                    tCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGaa                                                                                                                                  <  2:235326/77‑1 (MQ=255)                                                                                    tCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGaa                                                                                                                                  >  1:235326/1‑77 (MQ=255)                                                                                    tCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCgct                                                         >  6:193911/1‑149 (MQ=255)                                                                                         ttCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgccgagc                                                                                                                    >  3:284636/1‑86 (MQ=255)                                                                                         ttCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgccgagc                                                                                                                    <  4:284636/86‑1 (MQ=255)                                                                                         ttCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGaa                                                                                  <  3:334681/120‑1 (MQ=255)                                                                                         ttCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGaa                                                                                  >  4:334681/1‑120 (MQ=255)                                                                                           ccAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTgcgc                                                  >  6:134074/1‑149 (MQ=255)                                                                                            cAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCt                                                 >  2:215136/1‑149 (MQ=255)                                                                                                   aaCAGTTCGGGGGTGGGGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGcc                                          <  3:113502/149‑1 (MQ=255)                                                                                                       gTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGc                                      <  3:249400/149‑1 (MQ=255)                                                                                                          cGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAg                                    <  8:293204/148‑1 (MQ=255)                                                                                                          cGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAg                                    <  6:38081/148‑1 (MQ=255)                                                                                                           gggggTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAg                                    <  1:215136/148‑1 (MQ=255)                                                                                                             gggTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTCGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCCCCTTGCAGGTAt                                <  3:347175/149‑1 (MQ=255)                                                                                                                 ggtgCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGCTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGc                            <  5:350243/149‑1 (MQ=255)                                                                                                                          cAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAggg                   >  5:282949/1‑149 (MQ=255)                                                                                                                             gTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGttt                >  7:302444/1‑149 (MQ=255)                                                                                                                                  gCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcgga           >  5:223255/1‑149 (MQ=255)                                                                                                                                   cGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCt                                                 >  5:88259/1‑110 (MQ=255)                                                                                                                                   cGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCt                                                 <  6:88259/110‑1 (MQ=255)                                                                                                                                      tACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggacgga       >  8:248174/1‑149 (MQ=255)                                                                                                                                       cgttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGttt                <  2:294972/137‑1 (MQ=255)                                                                                                                                       cgttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCCCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGttt                >  1:294972/3‑139 (MQ=255)                                                                                                                                        gttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTg               >  7:117668/2‑139 (MQ=255)                                                                                                                                        gttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTg               <  8:117668/138‑1 (MQ=255)                                                                                                                                          tCCTGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCg                                                           <  4:201794/93‑1 (MQ=255)                                                                                                                                          tCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCg                                                           >  3:201794/1‑93 (MQ=255)                                                                                                                                           ccAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGc  >  5:320660/1‑149 (MQ=255)                                                                                                                                            cAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGc  <  4:8680/148‑1 (MQ=255)                                                                                                                                               cTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGc                                      >  3:238614/1‑109 (MQ=255)                                                                                                                                               cTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGc                                      <  4:238614/109‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                 CCAGGTCGAAGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTT‑CGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGC  >  NC_002947/3497185‑3497471

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