Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,972,973	A→G	6.4%	C301R (TGT→CGT)	PP_2600 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,972,973	0	A	G	6.4%	178.0 / 2.6	62	C301R (TGT→CGT)	PP_2600	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base G (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

ATAATGTGCGGGAAGTGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAGCTGCTCGCCCTGTGGCGTCGCGACC  >  NC_002947/2972837‑2973121                                                                                                                                         |                                                                                                                                                     aTAATGTGCGGGAAGTGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGggcgg                                                                                                                                          >  7:502/1‑149 (MQ=255)  tAATGTGCGGGAAGTGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGggcgg                                                                                                                                          >  4:207958/1‑148 (MQ=255)   aaTGTGCGGGAAGTGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGggcgg                                                                                                                                          >  4:214484/1‑147 (MQ=255)        gCGGGAAGTGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAgcgc                                                                                                                                    >  7:243103/1‑148 (MQ=255)        gCGGGAAGTGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAgcgc                                                                                                                                    >  3:299178/1‑148 (MQ=255)          gggAAGTGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAg                                                                                                                                 >  6:182836/1‑149 (MQ=255)           ggAAGTGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAgg                                                                                                                                >  1:317122/1‑149 (MQ=255)           ggAAGTGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAgg                                                                                                                                <  8:243103/149‑1 (MQ=255)           ggAAGTGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAgg                                                                                                                                <  2:317122/149‑1 (MQ=255)            gAAGTGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGt                                                                                                                               >  5:27768/1‑149 (MQ=255)                tGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTc                                                                                                                           <  8:133125/149‑1 (MQ=255)                  gCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGt                                                                                                                         <  6:27768/149‑1 (MQ=255)                  gCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGt                                                                                                                         >  4:228688/1‑149 (MQ=255)                  gCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGt                                                                                                                         <  5:182836/149‑1 (MQ=255)                  gCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGt                                                                                                                         <  6:122271/149‑1 (MQ=255)                   ccGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTg                                                                                                                        <  8:502/149‑1 (MQ=255)                    cGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTgg                                                                                                                       <  2:332090/149‑1 (MQ=255)                    cGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTgg                                                                                                                       >  1:332090/1‑149 (MQ=255)                            ccGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTg                                                                                                                 <  4:299178/147‑1 (MQ=255)                                       cGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGa                                                                                                    <  4:160538/149‑1 (MQ=255)                                        ggCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGaa                                                                                                   <  3:207958/149‑1 (MQ=255)                                           aCCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAg                                                                                                >  4:39103/1‑149 (MQ=255)                                           aCCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAg                                                                                                >  6:302705/1‑149 (MQ=255)                                           aCCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAg                                                                                                >  5:72500/1‑149 (MQ=255)                                             cTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATcgcg                                                                                                         >  2:156219/1‑138 (MQ=255)                                             cTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATcgcg                                                                                                         <  1:156219/138‑1 (MQ=255)                                             cTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGtt                                                                                              <  3:214484/149‑1 (MQ=255)                                             cTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGtt                                                                                              >  6:120403/1‑149 (MQ=255)                                               gtgAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGc                                                                                            >  2:285869/1‑149 (MQ=255)                                                tgAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCt                                                                                           >  1:74077/1‑149 (MQ=255)                                                           cccAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGc                                                                                 <  2:74077/148‑1 (MQ=255)                                                                  cACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCAt                                                                         >  3:52973/1‑149 (MQ=255)                                                                   aCTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTg                                                                                          >  3:284118/1‑131 (MQ=255)                                                                   aCTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTg                                                                                          <  4:284118/131‑1 (MQ=255)                                                                       gCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCgg                                                                                                      >  5:270915/1‑115 (MQ=255)                                                                       gCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCgg                                                                                                      <  6:270915/115‑1 (MQ=255)                                                                       gCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCt                                                                                           >  6:141940/1‑126 (MQ=255)                                                                       gCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCt                                                                                           <  5:141940/126‑1 (MQ=255)                                                                         gACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGa                                                                   >  6:275182/1‑148 (MQ=255)                                                                                ggTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCt                                                           <  5:275182/149‑1 (MQ=255)                                                                                ggTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCt                                                           <  5:120403/149‑1 (MQ=255)                                                                                 gTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTg                                                          >  5:267587/1‑149 (MQ=255)                                                                                  tCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTgc                                                         >  5:317496/1‑149 (MQ=255)                                                                                             cACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCgg                                                                                   <  4:223520/112‑1 (MQ=255)                                                                                             cACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCgg                                                                                   >  3:223520/1‑112 (MQ=255)                                                                                                  cgccgGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGcc                                         <  3:39103/149‑1 (MQ=255)                                                                                                  cgccgGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGcc                                         <  5:302705/149‑1 (MQ=255)                                                                                                            cTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGa                               <  6:72500/149‑1 (MQ=255)                                                                                                               gCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAgg                                                                                                                                >  7:312660/1‑49 (MQ=255)                                                                                                               gCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAgg                                                                                                                                <  8:312660/49‑1 (MQ=255)                                                                                                               gCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCgg                                >  5:64475/1‑145 (MQ=255)                                                                                                               gCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCgg                                <  6:64475/145‑1 (MQ=255)                                                                                                                 ccAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGa                               >  4:79946/1‑144 (MQ=255)                                                                                                                 ccAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGa                               <  3:79946/144‑1 (MQ=255)                                                                                                                 ccAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAgc                          >  6:158842/1‑149 (MQ=255)                                                                                                                     ccGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAgctgct                      <  6:267587/149‑1 (MQ=255)                                                                                                                              aTCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAGCTGCTCGCCCtgtg             <  3:312070/149‑1 (MQ=255)                                                                                                                                   cAGAACGCCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGCAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAGCTGCTCGCCCtgt              <  7:114593/143‑1 (MQ=255)                                                                                                                                   cAGAACGCCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGCAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAGCTGCTCGCCCtgt              <  7:119515/143‑1 (MQ=255)                                                                                                                                   cAGAACGCCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGCAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAGCTGCTCGCCCtgt              >  8:114593/1‑143 (MQ=255)                                                                                                                                   cAGAACGCCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGCAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAGCTGCTCGCCCtgt              >  8:119515/1‑143 (MQ=255)                                                                                                                                    aGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAGCTGCTCGCCCTGTGGCGTcg       >  8:75469/1‑149 (MQ=255)                                                                                                                                       acacCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAGCTGCTCGCCCTGTGGCGTCGCGa    <  4:52973/149‑1 (MQ=255)                                                                                                                                        cacCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAGCTGCTCGCCCTGTGGCGTCGCGAc   <  6:317496/149‑1 (MQ=255)                                                                                                                                         acCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAGCTGCTCGCCCTGTGGCGTCGCGAcc  >  6:192681/1‑149 (MQ=255)                                                                                                                                         |                                                                                                                                                     ATAATGTGCGGGAAGTGGCCGACGCTGCCGTCGCGATACGGCACCTGTGAAACCGGTGCCCAGGCCACTGGCGACTTCAGGTCGGTGGCAACCACGCCGCCGGCGGACTCGCCCAGCCGCAGGCCATCGCCAGAACACCCCTTGGGCGGCAGCGCCAGGTTGTCGTGGCCGCTGGCATCGCGCGGGAACAGTTGCTGGCGCCGGGCGTTGTCATTGGGGAAACCGCCTGCCGCCAGCACCACGGCCTTGGCGCGGATCAGCTGCTCGCCCTGTGGCGTCGCGACC  >  NC_002947/2972837‑2973121

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