Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	495,966	A→T	7.2%	L23Q (CTG→CAG)	PP_0408 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	495,966	0	A	T	7.2%	204.0 / 5.4	69	L23Q (CTG→CAG)	PP_0408	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base T (2/3);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CGGCCTGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTC  >  NC_002947/495820‑496107                                                                                                                                                   |                                                                                                                                               cGGCCTGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAgg                                                                                                                                              >  1:229638/1‑149 (MQ=255)        cAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGa                                                                                                                                      <  6:286181/149‑1 (MQ=255)           ggTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGCCCTGGCCGC‑CGatgc                                                                                                                                   <  2:308497/149‑4 (MQ=255)            gTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGCCGTGGCCGC‑CGatgcc                                                                                                                                  <  3:237484/149‑5 (MQ=255)             tCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTgg                                                                                                                                  >  3:169015/1‑148 (MQ=255)                gctgctCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCCCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGAGGAAg                                                                                                                               >  1:306373/1‑148 (MQ=255)                  tgctCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGAc                                                                                                                            >  1:160057/1‑149 (MQ=255)                    ctCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGAcc                                                                                                                           >  2:43152/1‑148 (MQ=255)                       aGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAg                                                                                                                        <  1:43152/148‑1 (MQ=255)                           ttGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAg                                                                                                                   >  8:165278/1‑149 (MQ=255)                               gcagcaGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGgcag                                                                                                               <  2:154270/149‑1 (MQ=255)                               gcagcaGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGgcag                                                                                                               >  2:19546/1‑149 (MQ=255)                                     gcCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACAtt                                                                                                         >  6:280643/1‑149 (MQ=255)                                      ccGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTg                                                                                                        <  2:291537/149‑1 (MQ=255)                                        gTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGcc                                                                                                      <  4:357045/149‑1 (MQ=255)                                         ttCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCa                                                                                                     <  3:31525/149‑1 (MQ=255)                                                  cAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGc                                                                                            <  5:280643/149‑1 (MQ=255)                                                  cAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGc                                                                                            <  2:278108/149‑1 (MQ=255)                                                   aGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGCCCTGGCCGC‑CGatg                                                                                                                                    >  7:280358/1‑106 (MQ=255)                                                   aGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGCCCTGGCCGC‑CGatg                                                                                                                                    <  8:280358/108‑3 (MQ=255)                                                        aTCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCCCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATaaa                                                                                      <  2:306373/149‑1 (MQ=255)                                                           ctcctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAg                                                                                                                        <  4:28337/112‑1 (MQ=255)                                                           ctcctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAg                                                                                                                        >  3:28337/1‑112 (MQ=255)                                                           ctcctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAgtg                                                                                   <  6:335462/149‑1 (MQ=255)                                                            tcctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGCCCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAgtgt                                                                                  >  6:74021/1‑149 (MQ=255)                                                            tcctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGcc                                                                                                                     >  3:20796/1‑114 (MQ=255)                                                            tcctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGcc                                                                                                                     <  4:20796/114‑1 (MQ=255)                                                              ctGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTgg                                                                                <  2:160057/149‑1 (MQ=255)                                                                 gCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACa                                                                             <  4:169015/149‑1 (MQ=255)                                                                     cGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGa                                                                          <  1:19546/148‑1 (MQ=255)                                                                       ttCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGc                                                                                                       >  6:213023/1‑117 (MQ=255)                                                                       ttCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGc                                                                                                       <  5:213023/117‑1 (MQ=255)                                                                        tCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTc                                                                       >  3:53675/1‑148 (MQ=255)                                                                          cAGCAAGGGTGCGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCACGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAAc                                                                    >  2:235493/1‑149 (MQ=255)                                                                           aGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTgg                                                                                >  6:92236/1‑136 (MQ=255)                                                                           aGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTgg                                                                                <  5:92236/136‑1 (MQ=255)                                                                             cAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAg                                                                                                                                     >  3:224606/1‑81 (MQ=255)                                                                             cAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAg                                                                                                                                     <  4:224606/81‑1 (MQ=255)                                                                              aaGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCaa                                                                <  2:229638/149‑1 (MQ=255)                                                                                  gTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGc                                                             >  8:163548/1‑148 (MQ=255)                                                                                       cGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAAt                                                       <  1:235493/149‑1 (MQ=255)                                                                                                tCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCgg                                               >  5:315442/1‑148 (MQ=255)                                                                                                  ggCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGc                                                             >  6:99185/1‑132 (MQ=255)                                                                                                  ggCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGc                                                             <  5:99185/132‑1 (MQ=255)                                                                                                         ggTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTaa                                       >  4:258970/1‑147 (MQ=255)                                                                                                         ggTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATa                                     <  7:165278/149‑1 (MQ=255)                                                                                                         ggTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATa                                     <  5:373718/149‑1 (MQ=255)                                                                                                             ccATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTAt                                 >  8:12353/1‑149 (MQ=255)                                                                                                             ccATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTAt                                 >  1:345042/1‑149 (MQ=255)                                                                                                               attgCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGAGCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTAttt                               <  2:162599/146‑1 (MQ=255)                                                                                                               aTGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTAttt                               >  3:318655/1‑149 (MQ=255)                                                                                                               aTGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTAttt                               >  1:344119/1‑149 (MQ=255)                                                                                                                tGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTg                              >  5:198798/1‑149 (MQ=255)                                                                                                                   cATTTGCAGGGGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCTTCGTGCGATAAAGTGTGGACAGGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGccc                           <  6:324269/149‑1 (MQ=255)                                                                                                                    aTTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTAtacaccgg                          >  1:332622/1‑142 (MQ=255)                                                                                                                     tttGCAGGTGGTTTTTCTGCATCGGCGGCCAGGGAGC‑CGAGTGGAAGGACAGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAAACCAAGGCCAGAATGAAACAGGAAGCATAATAGTATTTGCCCGt                         >  6:313797/1‑149 (MQ=255)                                                                                                                     tttGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGt                         >  8:258891/1‑149 (MQ=255)                                                                                                                           ggtggtTTTTCGGCATCGGCGGCCAGGGCGCGCGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCgg                                               >  1:232077/1‑122 (MQ=255)                                                                                                                           ggtggtTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCgg                                               <  2:232077/121‑1 (MQ=255)                                                                                                                           ggtggtTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGa                       >  7:47750/1‑145 (MQ=255)                                                                                                                           ggtggtTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGa                       <  8:47750/145‑1 (MQ=255)                                                                                                                            gtggtTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGaaa                                                   <  4:60944/116‑1 (MQ=255)                                                                                                                            gtggtTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGaaa                                                   >  3:60944/1‑116 (MQ=255)                                                                                                                               gtTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAt                <  4:53675/148‑1 (MQ=255)                                                                                                                                   ttCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAtgctg            <  2:345042/148‑1 (MQ=255)                                                                                                                                    tCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAtgctgc           <  3:258970/148‑1 (MQ=255)                                                                                                                                       gCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATgg       >  8:274057/1‑149 (MQ=255)                                                                                                                                           cggcggCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGt   >  8:260039/1‑149 (MQ=255)                                                                                                                                            ggcgccCTGGCCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTc  >  8:77586/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                               CGGCCTGCAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTC  >  NC_002947/495820‑496107

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