Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,553,210	A→C	5.6%	E84A (GAG→GCG)	PP_1363 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,553,210	0	A	C	5.6%	282.1 / 5.7	89	E84A (GAG→GCG)	PP_1363	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/47);  new base C (2/3);  total (39/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

TGGACCACCATCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACGGCCGTGGC  >  NC_002947/1553063‑1553357                                                                                                                                                    |                                                                                                                                                    tGGACCACCATCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAg                                                                                                                                                    >  4:167519/1‑149 (MQ=255) tGGACCACCATCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAg                                                                                                                                                    <  2:151337/149‑1 (MQ=255)     ccaccaTCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAgcgc                                                                                                                                                 >  3:178136/1‑148 (MQ=255)      caccaTCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGCGCTcc                                                                                                                                                >  4:392312/1‑148 (MQ=255)      caccaTCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGCGCTcc                                                                                                                                                <  3:392312/148‑1 (MQ=255)      caccaTCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAgcgc                                                                                                                                                 >  6:271540/1‑147 (MQ=255)      caccaTCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAgcgc                                                                                                                                                 >  6:87482/1‑147 (MQ=255)         caTCTGAACACTGCCTACCGCGCCTACGACCGCGGCGATTGCGCAGGGGGCATGCGCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGc                                                                                                                                            <  4:290576/149‑1 (MQ=255)         caTCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGc                                                                                                                                            <  3:5364/149‑1 (MQ=255)            cTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGtt                                                                                                                                         >  1:2835/1‑149 (MQ=255)               aaCAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGCGCTCCAGGcac                                                                                                                                          <  8:119442/145‑3 (MQ=255)               aaCAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGtt                                                                                                                                         <  3:179448/146‑1 (MQ=255)               aaCAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGtt                                                                                                                                         <  8:209846/146‑1 (MQ=255)               aaCAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGt                                                                                                                                          <  2:73556/145‑1 (MQ=255)               aaCAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTa                                                                                                                                      <  6:310420/149‑1 (MQ=255)                aCAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTAc                                                                                                                                     <  2:29375/149‑1 (MQ=255)                aCAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTAc                                                                                                                                     >  3:210749/1‑149 (MQ=255)                   aTGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTg                                                                                                                                  >  5:13653/1‑149 (MQ=255)                   aTGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTg                                                                                                                                  <  1:198457/149‑1 (MQ=255)                    tGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTg                                                                                                                                  >  3:51826/1‑148 (MQ=255)                      ccTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGa                                                                                                                                <  3:237454/148‑1 (MQ=255)                      ccTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGAc                                                                                                                               >  6:191123/1‑149 (MQ=255)                        tACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGAcgc                                                                                                                             >  5:59660/1‑149 (MQ=255)                        tACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGAcgc                                                                                                                             >  2:81641/1‑149 (MQ=255)                            gcgcCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTg                                                                                                                                  <  5:390374/140‑1 (MQ=255)                            gcgcCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTg                                                                                                                                  >  6:390374/1‑140 (MQ=255)                                cTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCaaa                                                                                                                     >  2:46290/1‑149 (MQ=255)                                   cGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAAc                                                                                                                    <  6:404220/147‑1 (MQ=255)                                   cGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTGGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAAc                                                                                                                    <  3:288982/147‑1 (MQ=255)                                     aCCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTAt                                                                                                                <  5:191123/149‑1 (MQ=255)                                         cggcgATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAg                                                                                                             >  5:247065/1‑148 (MQ=255)                                           gcgATTGCGCACGGGTCATGCTCGAACTGTCGCAGGGTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTAc                                                                                                          <  8:75984/149‑1 (MQ=255)                                            cgATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTAcc                                                                                                         <  1:94301/149‑1 (MQ=255)                                            cgATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTAcc                                                                                                         <  1:81641/149‑1 (MQ=255)                                             gATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCt                                                                                                        >  3:180227/1‑149 (MQ=255)                                              aTTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCtt                                                                                                       <  4:180227/149‑1 (MQ=255)                                               ttGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTa                                                                                                      <  7:140319/149‑1 (MQ=255)                                                  cgcACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATc                                                                                                    >  3:301475/1‑148 (MQ=255)                                                  cgcACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATc                                                                                                    >  8:294173/1‑148 (MQ=255)                                                       gggTCATGCTCGAACTGTCGCAAGTTGATCGCACCAGACGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGATGTACGTGGACGCGGCGCAAACCTATCAGTACCTTAACcag                                                                                                 >  5:121990/1‑146 (MQ=255)                                                             tGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCg                                                                                        <  4:301475/149‑1 (MQ=255)                                                                tCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCggg                                                                                      >  1:191053/1‑148 (MQ=255)                                                                    aCTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAAc                                                                                  >  4:205225/1‑148 (MQ=255)                                                                        tCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAata                                                                             >  5:365044/1‑149 (MQ=255)                                                                         cGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAatat                                                                            >  5:358251/1‑149 (MQ=255)                                                                          gCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATg                                                                           >  3:140394/1‑149 (MQ=255)                                                                           cAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGc                                                                          <  6:247065/149‑1 (MQ=255)                                                                           cAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGACACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGc                                                                          >  1:320683/1‑149 (MQ=255)                                                                                  aTCGCAGCAGCCGTGCCCGGCCGTTCATCCACGCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATcgc                                                                   <  6:121990/149‑1 (MQ=255)                                                                                   tCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATcgcg                                                                  >  3:230111/1‑149 (MQ=255)                                                                                   tCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATcgcg                                                                  >  2:238098/1‑149 (MQ=255)                                                                                     gCACCAGCCGTGCCCGTCCGTTCATCCACCCCGAAGTTTCACTGCTGCGCGGCCAGTGCCTTGAGCGCCAGGCGTTGTACGTGGACGCTGCTCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCgcgc                                                                 <  6:221279/148‑1 (MQ=255)                                                                                      cACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATc                                                                     <  2:431842/143‑1 (MQ=255)                                                                                      cACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATc                                                                     >  1:431842/1‑143 (MQ=255)                                                                                      cACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGccc                                                               >  8:370550/1‑149 (MQ=255)                                                                                            ccGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGgcgc                                                         <  2:23793/149‑1 (MQ=255)                                                                                              gTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGgcgcgc                                                       <  7:370550/149‑1 (MQ=255)                                                                                               tGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCACTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGCGCGCGcc                                                      >  4:251553/1‑149 (MQ=255)                                                                                                 cccGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTg                                                    <  4:210749/149‑1 (MQ=255)                                                                                                   cGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCa                                                  <  6:59660/149‑1 (MQ=255)                                                                                                        gTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGAccc                                             >  1:406366/1‑149 (MQ=255)                                                                                                             tCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGaa                                        <  3:167519/149‑1 (MQ=255)                                                                                                                aCCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGaaaaa                                     <  1:238098/149‑1 (MQ=255)                                                                                                                aCCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGaaaaa                                     >  3:263559/1‑149 (MQ=255)                                                                                                                    cGAAGTCTCACTGCTGCGCGGCCAGTGCCTGGAGCTCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGATTATTCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCg                                 <  6:43075/149‑1 (MQ=255)                                                                                                                       aGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTc                              <  1:46290/149‑1 (MQ=255)                                                                                                                           tCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCAttt                          <  5:87482/149‑1 (MQ=255)                                                                                                                             aCTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTgc                        >  3:194281/1‑149 (MQ=255)                                                                                                                             aCTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTgc                        >  2:271235/1‑149 (MQ=255)                                                                                                                               tgctgcGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTgcgc                      >  1:835/1‑149 (MQ=255)                                                                                                                               tgctgcGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTgcgc                      >  4:82145/1‑149 (MQ=255)                                                                                                                                 ctgcGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCATTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGcgcg                     <  2:320683/148‑1 (MQ=255)                                                                                                                                 ctgcGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCgg                    <  4:140394/149‑1 (MQ=255)                                                                                                                                    cgcgGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGt                   <  6:358251/147‑1 (MQ=255)                                                                                                                                    cgcgGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGCACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGc                 <  4:51826/149‑1 (MQ=255)                                                                                                                                         ccAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGaaa            >  2:208767/1‑149 (MQ=255)                                                                                                                                          caacgCCTGGCGCTCCAGGCGTTGTACGTGGATGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAAc           >  2:418479/5‑149 (MQ=255)                                                                                                                                          cAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAAc           >  3:407359/1‑149 (MQ=255)                                                                                                                                          cAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAAc           <  6:13653/149‑1 (MQ=255)                                                                                                                                            gTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACgg         <  5:271540/149‑1 (MQ=255)                                                                                                                                            gTGCCTGGAGCGCCAGGCGTTGTAAGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACgg         <  4:178136/149‑1 (MQ=255)                                                                                                                                              gCCTGGCGCTCCAGGCGTTGTACGTGGATGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACGGcc       <  1:418479/149‑1 (MQ=255)                                                                                                                                                cTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACGGCCGt     <  2:191053/149‑1 (MQ=255)                                                                                                                                                cTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACCCGGTCATTTGCGCGGTGCAGAAACGGCCGt     <  7:326493/149‑1 (MQ=255)                                                                                                                                                  ggAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACGGCCGTgg   <  4:230111/149‑1 (MQ=255)                                                                                                                                                   gAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACGGCCGTgg   <  7:294173/148‑1 (MQ=255)                                                                                                                                                   gAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACGGCCGTgg   <  3:82145/148‑1 (MQ=255)                                                                                                                                                   gAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACGGCCGTg    <  4:263559/147‑1 (MQ=255)                                                                                                                                                   gAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACGGCCGTGGc  <  6:365044/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TGGACCACCATCTGAACAATGCCTACCGCGCCTACGACCGCGGCGATTGCGCACGGGTCATGCTCGAACTGTCGCAGGTTGATCGCACCAGCCGTGCCCGGCCGTTCATCCACCCCGAAGTGTCACTGCTGCGCGGCCAGTGCCTGGAGCGCCAGGCGTTGTACGTGGACGCGGCGCAAACCTATCAGTACCTTATCCAGCAGTACCCGGGCAACGAATATGCCTATCGCGCCCGGGCGCGCCTGCAGACCCTGGAAAAACTCGGTCATTTGCGCGGTGCAGAAACGGCCGTGGC  >  NC_002947/1553063‑1553357

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