Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,182,992	A→T	8.2%	F248I (TTC→ATC)	PP_4563 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,182,992	0	A	T	8.2%	206.8 / 7.8	73	F248I (TTC→ATC)	PP_4563	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/39);  new base T (2/4);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.90e-01

GTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCT  >  NC_002947/5182851‑5183138                                                                                                                                               |                                                                                                                                                   ggggCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTg                                                                                                                                             <  6:261183/147‑1 (MQ=255) gTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTg                                                                                                                                             <  7:59675/149‑1 (MQ=255)    gCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGTACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTg                                                                                                                                          <  1:103082/149‑1 (MQ=255)     cGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                         <  4:210921/149‑1 (MQ=255)      gAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCATGCCCTGACTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTTAACAGCTTGCTGa                                                                                                                                         <  1:111284/149‑1 (MQ=255)         cagcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCaa                                                                                                                                     >  7:129630/1‑149 (MQ=255)          agcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGACGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAAc                                                                                                                                    <  8:321183/149‑1 (MQ=255)          agcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                    <  3:273087/149‑1 (MQ=255)           gcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACg                                                                                                                                   <  5:267197/149‑1 (MQ=255)               gACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAgg                                                                                                                               <  4:65661/149‑1 (MQ=255)               gACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAgg                                                                                                                               <  5:310503/149‑1 (MQ=255)                  ttCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                    <  2:237598/141‑1 (MQ=255)                  ttCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                    >  1:237598/1‑141 (MQ=255)                     cGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCAc                                                                                                                          <  2:368017/148‑1 (MQ=255)                          cTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCATGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGt                                                                                                                    <  1:262031/149‑1 (MQ=255)                             cGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACACCCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                                                                                                   <  8:354929/147‑1 (MQ=255)                             cGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                                                                                                   <  8:8878/147‑1 (MQ=255)                             cGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                                                                                                   <  3:36070/147‑1 (MQ=255)                               tACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGcac                                                                                                               >  1:55210/1‑149 (MQ=255)                               tACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGca                                                                                                                >  4:334859/1‑148 (MQ=255)                                  aGGCGTAAGGCCGTCAGGCTTTCTCGACGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACAcc                                                                                                            >  2:52266/1‑149 (MQ=255)                                     cGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACAccccc                                                                                                         >  4:258347/1‑149 (MQ=255)                                        aaGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTg                                                                                                                                          <  8:93754/113‑1 (MQ=255)                                        aaGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTg                                                                                                                                          >  7:93754/1‑113 (MQ=255)                                                      ttCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGc                                                                                        >  8:151251/1‑149 (MQ=255)                                                         tcGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTc                                                                                     >  2:297946/1‑149 (MQ=255)                                                           gATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTc                                                                                     <  5:287924/147‑1 (MQ=255)                                                            aTGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAg                                                                                  <  3:19416/149‑1 (MQ=255)                                                               cGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGgcg                                                                               <  6:146268/149‑1 (MQ=255)                                                                gATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAg                                                                                  <  4:318504/145‑1 (MQ=255)                                                                gATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAg                                                                                  >  3:318504/1‑145 (MQ=255)                                                                gATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGgcgc                                                                              >  3:165375/1‑149 (MQ=255)                                                                  tCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGgcgcgc                                                                            <  8:129630/149‑1 (MQ=255)                                                                     tCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTACAGGCGCGCCAt                                                                         >  5:242946/1‑149 (MQ=255)                                                                               cGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGctct                                                               >  5:267726/1‑149 (MQ=255)                                                                                  tCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTg                                                             >  4:126235/1‑148 (MQ=255)                                                                                  tCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTg                                                             >  6:6933/1‑148 (MQ=255)                                                                                  tCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGc                                                            <  8:95304/149‑1 (MQ=255)                                                                                  tCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGc                                                            <  1:297946/149‑1 (MQ=255)                                                                                     ggCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAgg                                                         <  3:67381/149‑1 (MQ=255)                                                                                     ggCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAgg                                                         <  5:6933/149‑1 (MQ=255)                                                                                        ccTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCg                                                                                          >  6:325490/1‑113 (MQ=255)                                                                                        ccTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCg                                                                                          <  5:325490/113‑1 (MQ=255)                                                                                        ccTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGgtg                                                      >  4:187593/1‑149 (MQ=255)                                                                                             cTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACCCCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCg                                                  >  4:37459/1‑149 (MQ=255)                                                                                             cTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                   <  4:165375/148‑1 (MQ=255)                                                                                             cTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCg                                                  >  4:42142/1‑149 (MQ=255)                                                                                             cTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCg                                                  >  4:105863/1‑149 (MQ=255)                                                                                                 cataacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAg                                              <  5:225325/146‑1 (MQ=255)                                                                                                  acaacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAgg                                             >  5:190324/1‑149 (MQ=255)                                                                                                    aacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGggc                                           >  6:368256/1‑149 (MQ=255)                                                                                                       ccGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAgg                                                         >  2:232764/1‑132 (MQ=255)                                                                                                       ccGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAgg                                                         <  1:232764/132‑1 (MQ=255)                                                                                                        cGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTT‑GGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCAAGGgcggcg                                       >  7:90404/1‑148 (MQ=255)                                                                                                             aTGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTAc                                  <  7:217949/149‑1 (MQ=255)                                                                                                                 ccTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGc                              <  3:37459/149‑1 (MQ=255)                                                                                                                 ccTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGc                              <  3:105863/149‑1 (MQ=255)                                                                                                                 ccTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGc                              <  3:42142/149‑1 (MQ=255)                                                                                                                  cttcttGAAGTTGCCCTTGCTGACGTTGAACAGCATGCTGATCAACCCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCt                             >  1:177566/1‑149 (MQ=255)                                                                                                                     cttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGc                           <  2:55210/148‑1 (MQ=255)                                                                                                                        gAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCa                          <  5:368256/146‑1 (MQ=255)                                                                                                                        gAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCa                          <  7:151251/146‑1 (MQ=255)                                                                                                                              gCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGt                  <  2:177566/148‑1 (MQ=255)                                                                                                                              gCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGt                  <  6:267726/148‑1 (MQ=255)                                                                                                                                ccTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGa                >  3:150494/1‑148 (MQ=255)                                                                                                                                     cctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCt                      >  4:254801/5‑137 (MQ=255)                                                                                                                                     cctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCt                      <  3:254801/133‑1 (MQ=255)                                                                                                                                       gACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGc        >  7:18458/1‑149 (MQ=255)                                                                                                                                        aCGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg         <  3:126235/147‑1 (MQ=255)                                                                                                                                        aCGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg         <  3:334859/147‑1 (MQ=255)                                                                                                                                        aCGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGcc       <  4:7653/149‑1 (MQ=255)                                                                                                                                          gTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGt     >  3:18097/1‑149 (MQ=255)                                                                                                                                             gAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCt  >  1:361183/1‑149 (MQ=255)                                                                                                                                              aaCAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGc                                                            <  7:137099/90‑1 (MQ=255)                                                                                                                                              aaCAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGc                                                            >  8:137099/1‑90 (MQ=255)                                                                                                                                               |                                                                                                                                                   GTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTG‑CTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCT  >  NC_002947/5182851‑5183138

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