Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,424,070	T→C	6.0%	N101N (AAT→AAC)	PP_3915 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,424,070	0	T	C	6.0%	194.2 / 2.3	67	N101N (AAT→AAC)	PP_3915	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base C (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GATCGTCGGCCACTACATGGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGAGTACGGCAACCCACACCGAGTGTTC  >  NC_002947/4423929‑4424216                                                                                                                                               |                                                                                                                                                   gaTCGTCGGCCACTACATTGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCcggagc                                                                                                                                             <  2:375350/149‑1 (MQ=255) gaTCGTCGGCCACTACATGGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGCCCTGTATGGC‑AATCcggagc                                                                                                                                             >  6:322704/1‑149 (MQ=255)     gtcgGCCACTACATGGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCcggagcggag                                                                                                                                         <  3:376758/149‑1 (MQ=255)             cTACATGGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCAt                                                                                                                                 <  3:305777/149‑1 (MQ=255)             cTACATGGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCAt                                                                                                                                 <  7:48939/149‑1 (MQ=255)              tACATGGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCAt                                                                                                                                 >  2:363325/1‑148 (MQ=255)                  tGGGCAAAACATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGAGGAGGCGCGCATgacc                                                                                                                             >  8:96299/1‑148 (MQ=255)                     gCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCcggagc                                                                                                                                             <  1:224981/129‑1 (MQ=255)                     gCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCcggagc                                                                                                                                             >  2:224981/1‑129 (MQ=255)                        aaCCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGgcgc                                                                                                                                     <  1:214888/134‑1 (MQ=255)                        aaCCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGgcgc                                                                                                                                     >  2:214888/1‑134 (MQ=255)                         acaaTGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGAt                                                                                                                     >  6:341503/4‑149 (MQ=255)                             tgtgCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGaa                                                                                                                 <  4:76597/149‑1 (MQ=255)                                   cgGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTg                                                                                                              >  2:379893/1‑146 (MQ=255)                                   cgGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAg                                                                                                           <  2:113161/149‑1 (MQ=255)                                      gTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGATTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAATTGAAGACg                                                                                                        <  5:202345/149‑1 (MQ=255)                                        gCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGAc                                                                                                         <  2:347911/146‑1 (MQ=255)                                        gCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGAc                                                                                                         >  1:347911/1‑146 (MQ=255)                                           cGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGa                                                                                                   >  6:254149/1‑149 (MQ=255)                                             gCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGAcc                                                                                                 <  5:254149/149‑1 (MQ=255)                                                  aaCTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGa                                                                                                          <  4:158248/135‑1 (MQ=255)                                                  aaCTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGa                                                                                                          >  3:158248/1‑135 (MQ=255)                                                  aaCTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAg                                                                                                    >  4:200902/1‑141 (MQ=255)                                                  aaCTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAg                                                                                                    <  3:200902/141‑1 (MQ=255)                                                     tACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGc                                                                                         >  8:363269/1‑149 (MQ=255)                                                       cAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGc                                                                                         >  2:130192/1‑147 (MQ=255)                                                        aacgccaaCCTGCGCGCCGCCATCGCCCAGGCGCTCACCGCTGAGGCAAATCCGCCATACCGCCCGGATTACGACCTGTATGGCAAATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGa                                                                                       <  1:328332/143‑1 (MQ=255)                                                        aaCGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCtg                                                                                               <  2:207576/140‑1 (MQ=255)                                                        aaCGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCtg                                                                                               >  1:207576/1‑140 (MQ=255)                                                            agaACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGgcgc                                                                                  <  2:165328/149‑1 (MQ=255)                                                              aaCCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCtt                                                                                <  1:379893/149‑1 (MQ=255)                                                                    cgcgcCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTtgggtg                                                                           <  7:164226/148‑1 (MQ=255)                                                                      cgccgcCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGc                                                                        >  1:107950/1‑149 (MQ=255)                                                                                gACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCAcg                                                              >  4:400964/1‑149 (MQ=255)                                                                                 aCCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCAcgg                                                             >  2:310915/1‑149 (MQ=255)                                                                                  ccAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCAcggc                                                            <  4:333100/149‑1 (MQ=255)                                                                                         cTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGc                                                     <  3:400964/149‑1 (MQ=255)                                                                                            aCCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AACCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTg                                                  >  2:99047/1‑149 (MQ=255)                                                                                            aCCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AACCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTg                                                  >  4:291345/1‑149 (MQ=255)                                                                                             ccGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTg                                                  >  2:137418/1‑148 (MQ=255)                                                                                              cGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGa                                                >  1:228114/1‑149 (MQ=255)                                                                                              cGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGa                                                <  1:363325/149‑1 (MQ=255)                                                                                               gCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGAc                                               <  7:363269/149‑1 (MQ=255)                                                                                               gCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTCGGTGGGCCGTGGGCACGGCGGAAGGCCTGGAc                                               <  7:96299/149‑1 (MQ=255)                                                                                                cTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGt                                                                                                               >  6:167496/1‑84 (MQ=255)                                                                                                cTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGt                                                                                                               <  5:167496/84‑1 (MQ=255)                                                                                                   aGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AACCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGt                                           <  3:291345/149‑1 (MQ=255)                                                                                                   aGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AACCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGt                                           <  1:99047/149‑1 (MQ=255)                                                                                                          tCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCAcgg                                                             >  1:153233/1‑124 (MQ=255)                                                                                                          tCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCAcgg                                                             <  2:153233/124‑1 (MQ=255)                                                                                                           ccGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGACGGCCTGGACATGTTCATGgcg                                   >  2:354178/1‑149 (MQ=255)                                                                                                                 gagCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGccgccg                             <  5:341503/149‑1 (MQ=255)                                                                                                                   gcgcCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGa                           <  1:310915/149‑1 (MQ=255)                                                                                                                       ccGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGAGTAc                       >  2:34427/1‑149 (MQ=255)                                                                                                                          gAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGAGTACgg                     >  7:75309/1‑148 (MQ=255)                                                                                                                           aGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGGGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGAGTACGGCa                   <  5:322704/149‑1 (MQ=255)                                                                                                                                gACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGAGTACGGCaa                  <  6:402447/145‑1 (MQ=255)                                                                                                                                gACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGAGTACGGCaa                  >  5:402447/1‑145 (MQ=255)                                                                                                                                        tGGC‑AATCCGGAGCGGAGGCGCGCATGACCGAGCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGAGTACGGCAACCCACACCGAgt      >  6:248099/1‑149 (MQ=255)                                                                                                                                        tGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGAGTACGGCAACCCACACCGAgt      >  1:58851/1‑149 (MQ=255)                                                                                                                                        tGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACAGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGAGTACGGCAACCCACACCGAgt      >  7:279968/1‑149 (MQ=255)                                                                                                                                          gC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGt                                           <  6:284264/110‑1 (MQ=255)                                                                                                                                          gC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGt                                           >  5:284264/1‑110 (MQ=255)                                                                                                                                          gC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGcgccgc                               <  3:172727/122‑1 (MQ=255)                                                                                                                                          gC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGcgccgc                               >  4:172727/1‑122 (MQ=255)                                                                                                                                          gC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGAGTACGGCAACCCACACCGAgtgt    >  2:311599/1‑149 (MQ=255)                                                                                                                                             aaTCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGAGTACGGCAACCCACACCGAGTGTTc  >  4:21118/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                   GATCGTCGGCCACTACATGGGCAAACCATGTGCGCGGGTGCTCGGCGAGAACTACAACGAGAACCTGCGCGCCGCCATCGACCAGGCGCTCACCGCTGAGGCATATCCGCCAGAGCGCCCGGAGTACGACCTGTATGGC‑AATCCGGAGCGGAGGCGCGCATGACCGACCTGATCGAAGTGAAGACGGCAGACCTGGTGGGCGAGGCGCTTGGGTGGGCCGTGGGCACGGCGGAAGGCCTGGACCTGTTCATGGCGCCGCCGGAGTACGGCAACCCACACCGAGTGTTC  >  NC_002947/4423929‑4424216

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑