Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,785,967	T→A	5.6%	L324Q (CTG→CAG)	PP_5068 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,785,967	0	T	A	5.6%	283.2 / 4.3	89	L324Q (CTG→CAG)	PP_5068	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (44/40);  new base A (2/3);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

AACTTCATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACT  >  NC_002947/5785819‑5786114                                                                                                                                                      |                                                                                                                                                    aaCTTCATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCt                                                                                                                                                     <  2:53870/149‑1 (MQ=255)    ttCATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGcc                                                                                                                                                  >  1:28718/1‑149 (MQ=255)     tCATCTGCAACCACTCCGACGACCGTGGACGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCg                                                                                                                                                 >  3:97516/1‑149 (MQ=255)       aTCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCt                                                                                                                                               >  6:63451/1‑149 (MQ=255)        tCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTg                                                                                                                                              >  4:204202/1‑149 (MQ=255)          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                                                                                                                                  >  6:101018/1‑148 (MQ=255)                     cgacgaCCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCAt                                                                                                                                  >  1:251074/1‑148 (MQ=255)                     cgacgaCCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGccc                                                                                                                                    <  1:281260/146‑1 (MQ=255)                     cgacgaCCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTGGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCAt                                                                                                                                  >  7:347857/1‑148 (MQ=255)                         gacCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCt                                                                                                                             >  5:46889/1‑149 (MQ=255)                         gacCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCt                                                                                                                             <  7:146504/149‑1 (MQ=255)                             gTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGAc                                                                                                                         <  2:392543/149‑1 (MQ=255)                              tGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGAcc                                                                                                                        >  6:150420/1‑149 (MQ=255)                                gCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCt                                                                                                                       >  8:398798/1‑148 (MQ=255)                                gCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCt                                                                                                                       >  6:17427/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                                                      <  2:28718/149‑1 (MQ=255)                                                                                                 aCTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCctgc                                                      <  3:289388/149‑1 (MQ=255)                                                                                                       caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGt                                                <  2:359182/149‑1 (MQ=255)                                                                                                       caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGt                                                >  3:10807/1‑149 (MQ=255)                                                                                                        accGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGtg                                               >  3:427536/1‑149 (MQ=255)                                                                                                          cGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTg                                                                   >  8:263263/1‑127 (MQ=255)                                                                                                          cGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTg                                                                   <  7:263263/127‑1 (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                >  4:119675/1‑144 (MQ=255)                                                                                                                                            ggCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTg                <  3:119675/144‑1 (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:119214/149‑1 (MQ=255)                                                                                                                                              ccTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTcc         >  3:106793/1‑149 (MQ=255)                                                                                                                                               cTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCg        >  8:334050/1‑149 (MQ=255)                                                                                                                                                 gTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGc       >  1:264059/1‑148 (MQ=255)                                                                                                                                                     cTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACt  >  4:409038/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AACTTCATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACT  >  NC_002947/5785819‑5786114

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