Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	495,969	G→C	7.4%	A22G (GCC→GGC)	PP_0408 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	495,969	0	G	C	7.4%	214.0 / 6.9	68	A22G (GCC→GGC)	PP_0408	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base C (2/3);  total (34/34)
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.54e-01

CAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTC  >  NC_002947/495827‑496107                                                                                                                                               |                                                                                                                                            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                                                                                                                                    <  8:280358/108‑3 (MQ=255)                                            aGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGCCCTGGCCGC‑CGatg                                                                                                                                    >  7:280358/1‑106 (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                                                                                <  5:92236/136‑1 (MQ=255)                                                                    aGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTgg                                                                                >  6:92236/1‑136 (MQ=255)                                                                      cAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAg                                                                                                                                     <  4:224606/81‑1 (MQ=255)                                                                      cAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAg                                                                                                                                     >  3:224606/1‑81 (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                                     <  5:373718/149‑1 (MQ=255)                                                                                                  ggTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATa                                     <  7:165278/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                       <  8:47750/145‑1 (MQ=255)                                                                                                                    ggtggtTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGa                       >  7:47750/1‑145 (MQ=255)                                                                                                                     gtggtTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGaaa                                                   >  3:60944/1‑116 (MQ=255)                                                                                                                     gtggtTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGaaa                                                   <  4:60944/116‑1 (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)                                                                                                                                            agagcgc‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAtgct             >  5:170888/4‑131 (MQ=255)                                                                                                                                            agagcgc‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCAtgct             <  6:170888/128‑1 (MQ=255)                                                                                                                                               |                                                                                                                                            CAGGGTCTGCTGCTCAGCCTTGGGCAGCAGCCGTTCCAGCGCCAGTGCATCCTCCTGGCTGCGTTCCAGCAAGGGTGGGCGTGGGGCCTCGGCAGCCGGTGCCATGGCATTTGCAGGTGGTTTTTCGGCATCGGCGGCCAGGGCGC‑CGAGTGGAAGGACCGAGGCCAGGCAGCACATTGCCAGCGTCGTGCGATAAAGTGTGGACATGAGTCAACCCAAGGCCAGAATGAAACCGGCAGCCTAATAGTATTTGCCCGTGACGGCCATGCTGCATGGCCGTC  >  NC_002947/495827‑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: ‑