Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,785,970	C→G	5.6%	P325R (CCG→CGG)	PP_5068 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,785,970	0	C	G	5.6%	289.2 / 6.0	89	P325R (CCG→CGG)	PP_5068	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (44/40);  new base G (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 = 8.75e-01

TTCATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGG  >  NC_002947/5785822‑5786116                                                                                                                                                      |                                                                                                                                                   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                                                                                                                                       >  7:255573/1‑149 (MQ=255)              aCTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGcc                                                                                                                                       >  3:266482/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                                                                                                                         >  6:17427/1‑148 (MQ=255)                             gCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCt                                                                                                                         >  8:398798/1‑148 (MQ=255)                               cGCTACAGCTATGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGAt                                                                                                                      >  7:228988/1‑149 (MQ=255)                               cGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGAt                                                                                                                      >  1:283839/1‑149 (MQ=255)                               cGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGAt                                                                                                                      <  5:17427/149‑1 (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                                                  <  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: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)                                                                                                                                                      |                                                                                                                                                   TTCATCTGCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGC‑ACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAGATGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGG  >  NC_002947/5785822‑5786116

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