Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,426,090	G→C	7.5%	W209S (TGG→TCG) ‡	PP_4765 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,426,090	0	G	C	7.5%	248.8 / 9.3	80	W209S (TGG→TCG) ‡	PP_4765	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/40);  new base C (3/3);  total (37/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 = 2.04e-01

TGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCAC  >  NC_002947/5425943‑5426233                                                                                                                                                    |                                                                                                                                                tGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTg                                                                                                                                                 >  5:297079/1‑148 (MQ=255)  gAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTg                                                                                                                                                 >  5:108702/1‑147 (MQ=255)   aaaGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGa                                                                                                                                               <  7:285765/148‑1 (MQ=255)         tCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcgc                                                                                                                                         >  2:281838/1‑148 (MQ=255)         tCGGTCGTGGTGCTGGCGCTGGGTGGCGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcgc                                                                                                                                         >  5:63608/1‑148 (MQ=255)           ggTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTa                                                                                                                                       <  4:80093/148‑1 (MQ=255)            gTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACt                                                                                                                                     <  4:56460/149‑1 (MQ=255)            gTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCGGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACt                                                                                                                                     <  5:127033/149‑1 (MQ=255)             tCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTc                                                                                                                                    <  2:108069/149‑1 (MQ=255)                tggtgCGGGCGCTGGGTGGTGCCAGCGGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATAGCTCGCCTTTGCAGCCGAGCACCTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAg                                                                                                                                 <  1:198037/149‑1 (MQ=255)                   tgCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGAc                                                                                                                              <  5:254951/149‑1 (MQ=255)                      tGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGAc                                                                                                                              >  3:12777/1‑146 (MQ=255)                      tGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGAc                                                                                                                              >  3:8417/1‑146 (MQ=255)                        gggcTGGGTGGTGGCAGCTGGGGGCGGCTGGGTTCTGACGGCAGCTGGGATGCGGTGCTGGGTGATCGGGCTGTGGTTATCTCGCCTTTGCTGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGActa                                                                                                                            <  4:8417/144‑3 (MQ=255)                        gcgcTGGGTTGTGCCAGCTTGGCGCGGATGGGTTCTGACGGAAGCTGGCTACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCTGTTTTTAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACaa                                                                                                                            <  4:12777/146‑1 (MQ=255)                          gcTGGGTGGTGCCAGCTGGGGGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTcg                                                                                                                       <  7:271485/149‑1 (MQ=255)                          gcTGGGTGGTGCCAGCTGGGGGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTcg                                                                                                                       <  7:275796/149‑1 (MQ=255)                           cTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACt                                                                                                                                     <  4:214045/134‑1 (MQ=255)                           cTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACt                                                                                                                                     >  3:214045/1‑134 (MQ=255)                                     ccAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCg                                                                                                            <  5:263862/149‑1 (MQ=255)                                         cTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTg                                                                                                         >  3:202108/1‑148 (MQ=255)                                              cgcgcCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGCAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTgaagaa                                                                                                    >  1:303749/1‑148 (MQ=255)                                                 gcCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTACtt                                                                                                                                    <  3:79537/113‑2 (MQ=255)                                                 gcCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTACtt                                                                                                                                    >  4:79537/1‑112 (MQ=255)                                                   cTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTAc                                                                                                                                      <  5:182346/109‑1 (MQ=255)                                                   cTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTAc                                                                                                                                      >  6:182346/1‑109 (MQ=255)                                                    tGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATTCCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcgc                                                                                                                                         <  1:84759/105‑1 (MQ=255)                                                    tGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcgc                                                                                                                                         >  2:84759/1‑105 (MQ=255)                                                    tGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGACTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCC                                                                                             >  6:198132‑M1/1‑149 (MQ=255)                                                     gggTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCC                                                                                            >  4:116158‑M1/1‑149 (MQ=255)                                                      ggTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCC                                                                                             <  4:202108‑M1/147‑1 (MQ=255)                                                       gTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCT                                                                                           >  5:30684‑M1/1‑148 (MQ=255)                                                           tGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAg                                                                                                                           >  2:278166/1‑112 (MQ=255)                                                           tGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAg                                                                                                                           <  1:278166/112‑1 (MQ=255)                                                            gACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcg                                                                                                                                          >  2:291866/1‑96 (MQ=255)                                                            gACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcg                                                                                                                                          <  1:291866/96‑1 (MQ=255)                                                            gACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGt                                                                                     >  2:97726/1‑149 (MQ=255)                                                               ggCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTccc                                                                                  >  7:93192/1‑149 (MQ=255)                                                                gcctcTGGCAACCGCTGCTGTCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCg                                                                                 <  8:148515/145‑1 (MQ=255)                                                                     tGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGc                                                                                       >  3:224977/1‑138 (MQ=255)                                                                     tGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGc                                                                                       <  4:224977/138‑1 (MQ=255)                                                                       gCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAgcg                                                                          >  5:320688/1‑149 (MQ=255)                                                                                       gAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAagggcga                                                          >  1:175897‑M2/1‑142 (MQ=255)                                                                                       gAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGAGTACCGGGCAGCTAAACGCGCAagg                                                              <  7:67382‑M2/145‑4 (MQ=255)                                                                                          cGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGtt                                                       <  1:281838/149‑1 (MQ=255)                                                                                             gCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGCGCGCACCGCTGAAGAACATCGCCCTtgcgcggtgcgctgcccgggacgctgagggcgatgttct                                                    <  6:108702‑M1/149‑40 (MQ=255)                                                                                             gCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGCAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAt                                                    <  2:303749/149‑1 (MQ=255)                                                                                              cTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAtt                                                   <  8:93192/149‑1 (MQ=255)                                                                                                  ggATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCa                                               >  6:331302/1‑149 (MQ=255)                                                                                                     tatCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAag                                                               >  8:144250‑M2/1‑128 (MQ=255)                                                                                                     tatCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAag                                                               <  7:144250‑M2/130‑3 (MQ=255)                                                                                                     tatCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCg                                            >  3:257907/1‑149 (MQ=255)                                                                                                       tctcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAtt                                                   >  1:281624/1‑140 (MQ=255)                                                                                                       tctcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAtt                                                   <  2:281624/140‑1 (MQ=255)                                                                                                        ctcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCa                                                                                         <  5:60465/101‑1 (MQ=255)                                                                                                        ctcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCa                                                                                         >  6:60465/1‑101 (MQ=255)                                                                                                         tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTtgcgcggtgcgctgcccgggacgctgagggcgatgttc                                                     >  8:163076‑M1/1‑98 (MQ=255)                                                                                                         tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTtgcgcggtgcgctgcccgggacgctgagggcgatgttc                                                     <  7:163076‑M1/136‑39 (MQ=255)                                                                                                         tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAagg                                                              <  4:32754‑M2/127‑4 (MQ=255)                                                                                                         tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAagg                                                              >  3:32754‑M2/1‑124 (MQ=255)                                                                                                           gCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAg                                      <  3:116158/149‑1 (MQ=255)                                                                                                            ccTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAgg                                     >  2:211360/1‑149 (MQ=255)                                                                                                            ccTTTGCAGCCGAGCAACTGCGGGTTTGAAGTTGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCTCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGATCATTCTCACCGCTCAAgg                                     <  7:269119/149‑1 (MQ=255)                                                                                                                tGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGTCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAggatt                                  <  6:320688/148‑4 (MQ=255)                                                                                                                tGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGt                                  <  5:331302/148‑1 (MQ=255)                                                                                                                 gCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTg                                 <  2:220337/148‑1 (MQ=255)                                                                                                                   aGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTgg                                >  5:103131/1‑147 (MQ=255)                                                                                                                   aGCCGAGCAACTGCGGTTTTGAAGTAGCGGACTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGaa                              <  5:198132/149‑1 (MQ=255)                                                                                                                     ccGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAgg                            <  6:103131/149‑1 (MQ=255)                                                                                                                     ccGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCAACGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAgg                            >  8:130748/1‑149 (MQ=255)                                                                                                                      cGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAgg                            >  1:242804/1‑148 (MQ=255)                                                                                                                           aaCTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAACTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCt                      >  7:212005/1‑149 (MQ=255)                                                                                                                            aCTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTg                     <  6:297079/149‑1 (MQ=255)                                                                                                                            aCTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTg                     >  1:45250/1‑149 (MQ=255)                                                                                                                              tGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGgt                   >  1:271293/1‑149 (MQ=255)                                                                                                                                cGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGgt                   <  1:211360/147‑1 (MQ=255)                                                                                                                                  gTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTAt               <  6:30684/149‑1 (MQ=255)                                                                                                                                   ttttGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATg              >  3:251589/1‑149 (MQ=255)                                                                                                                                      ttgaGTAGCGCGCTCCAGCCCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCt           >  7:224235/4‑149 (MQ=255)                                                                                                                                        aaGTAGAGTTCTGGAGCGATCTACTCAAGGACAAGTTCGCTGGCGAACCGATGAAGAACATCGCCTTCATCGTCCCGGGCAGCGAACCGCGAAAGGGCTAGTTAATTATCAACGATAAATGCGTGGATGGCAGCCTTGTGTATGCCTg          <  8:67382/148‑1 (MQ=255)                                                                                                                                              cGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCa   >  1:204673/1‑149 (MQ=255)                                                                                                                                               gggCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCAc  <  2:204673/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                TGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCAC  >  NC_002947/5425943‑5426233

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: ‑