Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,645,307	C→T	6.7%	A220A (GCC→GCT)	PP_4954 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,645,307	0	C	T	6.7%	217.6 / 4.2	74	A220A (GCC→GCT)	PP_4954	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base T (3/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CCGGTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGA  >  NC_002947/5645164‑5645455                                                                                                                                                |                                                                                                                                                      ccGGTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCt                                                                                                                                                  <  4:359891/149‑1 (MQ=255)   ggTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGa                                                                                                                                                <  6:211911/149‑1 (MQ=255)         gcagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAg                                                                                                                                          >  6:235824/1‑149 (MQ=255)          cagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGc                                                                                                                                         >  1:15698/1‑149 (MQ=255)           agagCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGc                                                                                                                                         >  3:89293/1‑148 (MQ=255)                       gACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGat‑catc                                                                                                                           <  1:368755/149‑1 (MQ=255)                                 cGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCggg                                                                                                                 >  4:24917/1‑149 (MQ=255)                                   aGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCTAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTg                                                                                                               <  2:231393/149‑1 (MQ=255)                                        ggCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGCT‑CATCCTGCACCGGGTGGGCaa                                                                                                          >  7:228837/1‑149 (MQ=255)                                        ggCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGc                                                                                                            <  1:42501/147‑1 (MQ=255)                                         gCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTgg                                                                                                              <  7:297353/144‑1 (MQ=255)                                         gCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTgg                                                                                                              >  8:297353/1‑144 (MQ=255)                                           aGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAg                                                                                                       <  3:24917/149‑1 (MQ=255)                                                ccTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGccaa                                                                                                   >  7:182973/1‑148 (MQ=255)                                                ccTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGccaa                                                                                                   <  8:182973/148‑1 (MQ=255)                                                   gAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTgg                                                                                                              <  4:226748/134‑1 (MQ=255)                                                   gAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTgg                                                                                                              >  3:226748/1‑134 (MQ=255)                                                      gCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACt                                                                                            >  5:91198/1‑149 (MQ=255)                                                        tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTc                                                                                           >  1:260435/1‑148 (MQ=255)                                                            ggCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAAc                                                                                                         <  1:318095/130‑1 (MQ=255)                                                            ggCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAAc                                                                                                         >  2:318095/1‑130 (MQ=255)                                                            ggCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGGAgg                                                                                      >  3:116677/1‑149 (MQ=255)                                                               cAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAg                                                                                                                                          >  2:174604/1‑95 (MQ=255)                                                               cAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAg                                                                                                                                          <  1:174604/95‑1 (MQ=255)                                                               cAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGGAGGTaa                                                                                   <  4:116677/149‑1 (MQ=255)                                                               cAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGt                                                                                     >  3:107633/1‑147 (MQ=255)                                                               cAGCCTGTTTGAACCGGGCCGCAAGCTCGACGTCGAGGCACTGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAATGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTaa                                                                                   <  5:256106/149‑1 (MQ=255)                                                                 gCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGc                                                                                 <  6:91198/149‑1 (MQ=255)                                                                  ccTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGcc                                                                                <  4:169908/149‑1 (MQ=255)                                                                  ccTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGcc                                                                                <  2:260435/149‑1 (MQ=255)                                                                      tttGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTa                                                                            >  3:216163/1‑149 (MQ=255)                                                                       ttGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTAc                                                                           >  7:198562/1‑149 (MQ=255)                                                                        tGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTa                                                                            >  3:78791/1‑147 (MQ=255)                                                                        tGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTa                                                                            <  4:78791/147‑1 (MQ=255)                                                                        tGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑AATCCTGCACAGGGTGCGCAACAGCCCACCAAACTCGAAGGTAAGCCGCTAc                                                                           >  2:363851/1‑148 (MQ=255)                                                                         gAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTAcc                                                                          <  5:235824/148‑1 (MQ=255)                                                                             cGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAg                                                                     >  2:232227/1‑149 (MQ=255)                                                                             cGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAg                                                                     <  4:216163/149‑1 (MQ=255)                                                                                  cGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                 <  1:379564/148‑1 (MQ=255)                                                                                         tCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGt                                                                                                                <  5:332020/94‑1 (MQ=255)                                                                                         tCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGt                                                                                                                >  6:332020/1‑94 (MQ=255)                                                                                               tCGAGGCACGGCGCCTGGAACTGGCCCAGCATAACGGCAAGGTGGTGCCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTTGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAg                                                   <  8:261315/149‑1 (MQ=255)                                                                                                cGAGGCACTGCGCCTGTAACTTGCCCAGCAGAACGGCAAGGTGGTGTCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAg                                                   <  8:266410/148‑1 (MQ=255)                                                                                                 gAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAg                                                                                                       >  5:167623/1‑95 (MQ=255)                                                                                                 gAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAgg                                                  <  3:167595/148‑1 (MQ=255)                                                                                                 gAGGCACGGCGACTGGAACTGGCCCAGCAGAACTGCAAGGTGGTGGACGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAg                                                                                                       <  6:167623/95‑1 (MQ=255)                                                                                                  aGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGat                                                <  2:15698/149‑1 (MQ=255)                                                                                                      aCGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATcc                                            <  1:38201/149‑1 (MQ=255)                                                                                                      aCGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATc                                             >  3:93842/1‑148 (MQ=255)                                                                                                           gcCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAg                                       >  3:65129/1‑149 (MQ=255)                                                                                                               ggAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAg                                                                     <  1:220502/115‑1 (MQ=255)                                                                                                                gAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCgt                                    <  4:93842/147‑1 (MQ=255)                                                                                                                gAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCgt                                    <  3:137150/147‑1 (MQ=255)                                                                                                                      gCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAATGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTg                             <  1:363851/148‑1 (MQ=255)                                                                                                                        ccAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGc                                      <  6:7573/137‑1 (MQ=255)                                                                                                                        ccAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGc                                      >  5:7573/1‑137 (MQ=255)                                                                                                                        ccAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCt                          >  2:312623/1‑149 (MQ=255)                                                                                                                          agcagAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCg                        >  2:333596/1‑149 (MQ=255)                                                                                                                            cagAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCa                      <  4:107633/149‑1 (MQ=255)                                                                                                                             agAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCAc                     >  5:284525/1‑149 (MQ=255)                                                                                                                                  ggCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAccc                <  8:228837/149‑1 (MQ=255)                                                                                                                                    cAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGc                                      >  3:155740/1‑125 (MQ=255)                                                                                                                                    cAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGc                                      <  4:155740/125‑1 (MQ=255)                                                                                                                                    cAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAcccc               >  4:253779/1‑148 (MQ=255)                                                                                                                                        gtggtgGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACaa          >  5:256054/1‑149 (MQ=255)                                                                                                                                         tggtggCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAAc         <  8:198562/149‑1 (MQ=255)                                                                                                                                         ccttggCTGCGTTCAGCGCGGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAAc         >  1:39673/4‑149 (MQ=255)                                                                                                                                           ttggCTGCGTTCAGCGCGGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAccc                <  8:99996/139‑1 (MQ=255)                                                                                                                                           ttggCTGCGTTCAGCGCGGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAccc                >  7:99996/2‑140 (MQ=255)                                                                                                                                           ttggCTGCGTTCAGCGCGGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAAc         >  7:54876/2‑147 (MQ=255)                                                                                                                                           ttggCTGCGTTCAGCGCGGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTCTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAAc         <  8:54876/146‑1 (MQ=255)                                                                                                                                           gtggCCGCGCTGAACGCAGCCAAGGAGATGCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg        >  1:28958/1‑149 (MQ=255)                                                                                                                                           gtggCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAc                  <  2:81777/138‑1 (MQ=255)                                                                                                                                           gtggCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAc                  >  1:81777/1‑138 (MQ=255)                                                                                                                                                cgcgCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCa                                           <  5:104679/108‑1 (MQ=255)                                                                                                                                                cgcgCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCa                                           >  6:104679/1‑108 (MQ=255)                                                                                                                                                cgcgCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGa  >  8:253375/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                      CCGGTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCAAGGAGAT‑CATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGA  >  NC_002947/5645164‑5645455

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