Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,785,975	T→G	5.6%	Y327D (TAC→GAC) ‡	PP_5068 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,785,975	0	T	G	5.6%	281.6 / 6.8	90	Y327D (TAC→GAC) ‡	PP_5068	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/43);  new base G (2/3);  total (44/46)
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.54e-01

TGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGC  >  NC_002947/5785828‑5786122                                                                                                                                                     |                                                                                                                                                    tGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTa                                                                                                                                                    >  4:146208/1‑149 (MQ=255)       cagtCCGACGACGGTGGCCGCTTCATCTACGCCAACCAGGTGCCCTTCGCCCAGTTGAACCTCAGCGCGGTGGGCCAGGC‑ACTGACCACCGTCTTCGAAGTTGAGCCTGTGAAGGATGCTTTGGGCCTGTTCCTGCCGCTGTACCAGGc                                                                                                                                              <  2:370264/146‑1 (MQ=255)        aCTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGcc                                                                                                                                             >  3:266482/1‑149 (MQ=255)        aCTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGcc                                                                                                                                             >  7:255573/1‑149 (MQ=255)        aCTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGcc                                                                                                                                             <  6:113210/149‑1 (MQ=255)         cTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGccc                                                                                                                                            <  4:11103/149‑1 (MQ=255)         cTCCGACGACCGTGGACGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCACAGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGccc                                                                                                                                            <  5:31777/149‑1 (MQ=255)          tCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCa                                                                                                                                           >  6:394393/1‑149 (MQ=255)          tCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGccc                                                                                                                                            >  2:281260/1‑148 (MQ=255)            cgacgaCCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCAtt                                                                                                                                         <  1:271894/149‑1 (MQ=255)            cgacgaCCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCAt                                                                                                                                          >  1:251074/1‑148 (MQ=255)            cgacgaCCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCAt                                                                                                                                          >  6:101018/1‑148 (MQ=255)            cgacgaCCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGccc                                                                                                                                            <  1:281260/146‑1 (MQ=255)            cgacgaCCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTGGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCAt                                                                                                                                          >  7:347857/1‑148 (MQ=255)                gacCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCt                                                                                                                                     <  7:146504/149‑1 (MQ=255)                gacCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCt                                                                                                                                     >  5:46889/1‑149 (MQ=255)                    gTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGAc                                                                                                                                 <  2:392543/149‑1 (MQ=255)                     tGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGAcc                                                                                                                                >  6:150420/1‑149 (MQ=255)                       gCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCt                                                                                                                               >  6:17427/1‑148 (MQ=255)                       gCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCt                                                                                                                               >  8:398798/1‑148 (MQ=255)                         cGCTACAGCTATGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGAt                                                                                                                            >  7:228988/1‑149 (MQ=255)                         cGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGAt                                                                                                                            <  5:17427/149‑1 (MQ=255)                         cGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGAt                                                                                                                            >  1:283839/1‑149 (MQ=255)                          gCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATg                                                                                                                           <  5:101018/149‑1 (MQ=255)                              cAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATgcg                                                                                                                         <  7:87480/147‑1 (MQ=255)                               aGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCg                                                                                                                      <  2:283839/149‑1 (MQ=255)                                gCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCGGATGACAACCGTGATCGAAGTTGCGCCACTGAAGGAACCGTTGGGCCTGTTCCTGCCGCTGTA‑CAGGCCCATTACCTGGCCCTGATGCGCCgg                                                                                                                     >  1:337521/1‑149 (MQ=255)                                    cGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGgcgc                                                                                                                  <  4:18623/148‑1 (MQ=255)                                    cGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGgcg                                                                                                                   <  2:294125/147‑1 (MQ=255)                                       cAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAg                                                                                                              <  8:327015/149‑1 (MQ=255)                                        aaCCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAgg                                                                                                             >  7:110875/1‑149 (MQ=255)                                        aaCCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAgg                                                                                                             <  3:204202/149‑1 (MQ=255)                                         aCCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAggg                                                                                                            <  8:347857/149‑1 (MQ=255)                                              gggCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCTAAGTTGAGCCACTGAAGTAAGCGTTGTTCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGac                                                                                                       <  1:349239/147‑1 (MQ=255)                                              gTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGac                                                                                                       <  7:319031/149‑1 (MQ=255)                                                 cccATCGCCAACTGGAACCTCAGCGAGCTGGACCAGGC‑ACTGACAACCGTGATCGAAGTTGATCCACTGAAGGAAGCGTTGGGCATGTTCCTGCCGCTGTACCAAGCCCATTACATGGACATGATGCGCCGGCGCCTATGGCTGaccaa                                                                                                    >  8:13361/1‑148 (MQ=255)                                                   cATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCTAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCAccg                                                                                                  <  6:362527/149‑1 (MQ=255)                                                    aTCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCAccgc                                                                                                 <  8:110875/149‑1 (MQ=255)                                                       gCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCAccgccg                                                                                               <  2:251074/148‑1 (MQ=255)                                                        cccACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGaa                                                                                             <  4:338954/149‑1 (MQ=255)                                                        cccACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTGGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCAccgccg                                                                                               >  5:415730/1‑147 (MQ=255)                                                         ccACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCAccgccg                                                                                               >  4:289388/1‑146 (MQ=255)                                                         ccACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAg                                                                                            <  7:398798/149‑1 (MQ=255)                                                           aCTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAgat                                                                                          <  4:266482/149‑1 (MQ=255)                                                            cTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAgatg                                                                                         <  5:150420/149‑1 (MQ=255)                                                            cTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAgat                                                                                          >  8:114401/1‑148 (MQ=255)                                                               gAACCGCAGCGCGCTGGCCCAGGC‑ACGTACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACTTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGatgatg                                                                                      <  8:374010/149‑1 (MQ=255)                                                                aaCCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGAtgatga                                                                                     <  8:228988/149‑1 (MQ=255)                                                                   cTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACAt                                                                                  <  3:146208/149‑1 (MQ=255)                                                                   cTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACAt                                                                                  >  3:357606/1‑149 (MQ=255)                                                                   cTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACAAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGGAGATGATGACAt                                                                                  >  2:348942/1‑149 (MQ=255)                                                                    tCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACAt                                                                                  >  4:2988/1‑148 (MQ=255)                                                                      aGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACAt                                                                                  <  3:2988/146‑1 (MQ=255)                                                                               gCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTc                                                                        >  7:372184/1‑147 (MQ=255)                                                                               gCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTc                                                                        >  6:217048/1‑147 (MQ=255)                                                                                cccAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAg                                                                     <  8:372184/149‑1 (MQ=255)                                                                                  cAGGC‑ACTGACCACCGTGATCGAAGTAGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCACATTACCTGGACCTGATGCGGCGGCGCCTAGGGCTGcccacc                                                                                                   >  7:151826/1‑117 (MQ=255)                                                                                        aCTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCctgc                                                              <  3:289388/149‑1 (MQ=255)                                                                                        aCTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCctgc                                                              <  2:28718/149‑1 (MQ=255)                                                                                              caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGt                                                        >  3:10807/1‑149 (MQ=255)                                                                                              caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGt                                                        <  2:359182/149‑1 (MQ=255)                                                                                               accGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGtg                                                       >  3:427536/1‑149 (MQ=255)                                                                                                 cGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTg                                                                           <  7:263263/127‑1 (MQ=255)                                                                                                 cGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTg                                                                           >  8:263263/1‑127 (MQ=255)                                                                                                   tGATCGAAGTTGAGCCAATGAAGGAAGCGTTGGGCATGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCAGAGGGCTGACCACCGCCGAAGATGATGACATGGTGATGGTCGATCGCATGCTGCATtgcagg                                                   >  8:271929/1‑147 (MQ=255)                                                                                                        gAAGTTGAGCGACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGaa                                                                                             <  1:73363/102‑1 (MQ=255)                                                                                                        gAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGaa                                                                                             >  2:73363/1‑102 (MQ=255)                                                                                                        gAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAAcg                                              >  4:287473/1‑149 (MQ=255)                                                                                                            ttGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGt                                          <  5:394393/149‑1 (MQ=255)                                                                                                                  cACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTg                                    <  6:415730/149‑1 (MQ=255)                                                                                                                          gaagCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAg                            >  3:142491/1‑149 (MQ=255)                                                                                                                            agCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATAACCAGGACCTGATGCGCCGGCGACTAGGGCTGACCACCGCCGATGATGATGACATGGTGCTGGTCGACAGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGcc                          >  7:396861/1‑149 (MQ=255)                                                                                                                              taatGGGCCTGGTACAGCGGCAGGAACAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCa                                                          <  7:197986/112‑1 (MQ=37)                                                                                                                              taatGGGCCTGGTACAGCGGCAGGAACAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCa                                                          >  8:197986/4‑115 (MQ=37)                                                                                                                               gTTTGGCCTGTTCCTGCCGCTTTACCATGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTTACCACCGCCGAAGATGATGACATTGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTg                        <  4:97516/148‑1 (MQ=255)                                                                                                                                   ggCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTg                        <  3:119675/144‑1 (MQ=255)                                                                                                                                   ggCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTg                        >  4:119675/1‑144 (MQ=255)                                                                                                                                    gCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGtgcttc                  >  8:183515/1‑149 (MQ=255)                                                                                                                                    cccTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCACACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTCCTGCAGTGCATGCAACGCGGTCGCGTCGACTACAGCATcttctta                  >  2:86947/2‑148 (MQ=255)                                                                                                                                     ccTGTTCCTGCCTCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTcc                 <  7:13361/149‑1 (MQ=255)                                                                                                                                     ccTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTcc                 >  3:106793/1‑149 (MQ=255)                                                                                                                                     ccTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTcc                 <  3:119214/149‑1 (MQ=255)                                                                                                                                      cTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCg                >  8:334050/1‑149 (MQ=255)                                                                                                                                        gTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGc               >  1:264059/1‑148 (MQ=255)                                                                                                                                            cTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACt          >  4:409038/1‑149 (MQ=255)                                                                                                                                              gCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTgg        >  3:343825/1‑149 (MQ=255)                                                                                                                                               ccGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTgg        <  4:10807/148‑1 (MQ=255)                                                                                                                                               ccGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTgg        <  3:143488/148‑1 (MQ=255)                                                                                                                                               ccGCTGTA‑CAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTgg        <  2:337521/147‑1 (MQ=255)                                                                                                                                                 gCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGt      >  1:99292/1‑148 (MQ=255)                                                                                                                                                     tACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGagc  <  3:409038/148‑1 (MQ=255)                                                                                                                                                     tACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGag   <  5:217048/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGAGC  >  NC_002947/5785828‑5786122

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