Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,785,973	T→A	5.5%	L326Q (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,973	0	T	A	5.5%	285.6 / 4.4	91	L326Q (CTG→CAG)	PP_5068	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (43/43);  new base A (2/3);  total (45/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 = 5.09e-01

ATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGT  >  NC_002947/5785825‑5786118                                                                                                                                                      |                                                                                                                                                  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                                                                                                                                      >  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                                                                                                                                 >  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                                                                    >  6:217048/1‑147 (MQ=255)                                                                                  gCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTc                                                                    >  7:372184/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                                                    <  2:359182/149‑1 (MQ=255)                                                                                                 caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGt                                                    >  3:10807/1‑149 (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                                                      >  8:197986/4‑115 (MQ=37)                                                                                                                                 taatGGGCCTGGTACAGCGGCAGGAACAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCa                                                      <  7:197986/112‑1 (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)                                                                                                                                                      |                                                                                                                                                  ATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGT  >  NC_002947/5785825‑5786118

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