Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,200,465	G→A	7.5%	V41M (GTG→ATG) ‡	gspG →	protein secretion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,200,465	0	G	A	7.5%	194.7 / 4.0	66	V41M (GTG→ATG) ‡	gspG	protein secretion protein
Reads supporting (aligned to +/- strand):  ref base G (32/29);  new base A (2/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

ATGCAACTCAACCAGGGCCTGCAATTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCA  >  NC_002947/1200319‑1200613                                                                                                                                                   |                                                                                                                                                     aTGCAACTCAACCAGGGCCTGCAATTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTg                                                                                                                                                    >  2:299757/1‑149 (MQ=255)          aaCCAGGGCCTGCAATTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCACGCTGGGcgc                                                                                                                                            <  3:95585/148‑3 (MQ=255)               gggCCTGCAATTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGa                                                                                                                                      <  7:297873/149‑1 (MQ=255)               gggCCTGCAATTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGa                                                                                                                                      >  7:235646/1‑149 (MQ=255)                    tGCAATTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAgg                                                                                                                                 >  3:67680/1‑149 (MQ=255)                    tGCAATTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAgg                                                                                                                                 <  3:2936/149‑1 (MQ=255)                    tGCAATTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAgg                                                                                                                                 >  1:286376/1‑149 (MQ=255)                     gCAATTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGc                                                                                                                                >  8:145674/1‑149 (MQ=255)                      cAATTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAgg                                                                                                                                 >  4:2936/1‑147 (MQ=255)                      cAATTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGc                                                                                                                                >  2:3771/1‑148 (MQ=255)                        aTTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCAt                                                                                                                             <  7:156617/149‑1 (MQ=255)                             gacgacGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAg                                                                                                                         >  8:38745/1‑148 (MQ=255)                                gacgaGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAg                                                                                                                      >  1:33321/1‑148 (MQ=255)                                     ggagTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGt                                                                                                                 >  8:31899/1‑148 (MQ=255)                                       agTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGAt                                                                                                              <  4:67680/149‑1 (MQ=255)                                        gTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGaa                                                                                                                          >  1:70957/1‑136 (MQ=255)                                        gTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGaa                                                                                                                          <  2:70957/136‑1 (MQ=255)                                        gTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATg                                                                                                             <  2:286376/149‑1 (MQ=255)                                           ttGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATTGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGc                                                                                                           <  2:82839/148‑1 (MQ=255)                                               cATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGAcc                                                                                                      >  2:41227/1‑149 (MQ=255)                                                aTGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAAcc                                                                                                                                          >  3:165851/1‑112 (MQ=255)                                                aTGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAAcc                                                                                                                                          <  4:165851/112‑1 (MQ=255)                                                            cGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGa                                                                                         <  1:299757/149‑1 (MQ=255)                                                               aaCCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAg                                                                                        >  7:82493/1‑147 (MQ=255)                                                                  cGCCAACGCGGTTTCACCCTCATGTAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCACGCTGGGcgcta                                                                                                                                          <  4:297371/94‑5 (MQ=255)                                                                  cGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTgg                                                                                          >  8:251022/1‑142 (MQ=255)                                                                  cGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTgg                                                                                          <  7:251022/142‑1 (MQ=255)                                                                  cGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGgc                                                                                   <  8:235646/149‑1 (MQ=255)                                                                     cAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCt                                                                                >  4:140771/1‑149 (MQ=255)                                                                      aaCGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTg                                                                               <  2:33321/149‑1 (MQ=255)                                                                      aaCGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTg                                                                               >  6:191032/1‑149 (MQ=255)                                                                      aaCGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTg                                                                               <  7:38745/149‑1 (MQ=255)                                                                           ggTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACAt                                                                          <  7:145674/149‑1 (MQ=255)                                                                             tttCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGt                                                                        <  2:210762/149‑1 (MQ=255)                                                                                cACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGa                                                                                         <  5:145628/129‑1 (MQ=255)                                                                                cACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGa                                                                                         >  6:145628/1‑129 (MQ=255)                                                                                   ccTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAg                                                                                                                   <  8:218735/100‑1 (MQ=255)                                                                                   ccTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAg                                                                                                                   >  7:218735/1‑100 (MQ=255)                                                                                   ccTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATgg                                                                                                            <  4:77065/107‑1 (MQ=255)                                                                                     tCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCCCCCAGGAGCAGGCGCTGGACATGTATCGCCTg                                                                >  8:143636/1‑149 (MQ=255)                                                                                     tCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTg                                                                                           <  5:211780/122‑1 (MQ=255)                                                                                     tCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTg                                                                                           >  6:211780/1‑122 (MQ=255)                                                                                       aTGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTg                                                                <  1:3771/147‑1 (MQ=255)                                                                                         ggAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGAc                                                             >  8:125328/1‑148 (MQ=255)                                                                                                      gtgATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTccc                                               >  4:300604/1‑149 (MQ=255)                                                                                                      gtgATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTccc                                               <  8:82493/149‑1 (MQ=255)                                                                                                            tttATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACAt                                                                          <  1:157646/116‑1 (MQ=255)                                                                                                            tttATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACAt                                                                          >  2:157646/1‑116 (MQ=255)                                                                                                              tATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAAc                                        >  7:113023/1‑148 (MQ=255)                                                                                                              tATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAAc                                        >  1:164529/1‑148 (MQ=255)                                                                                                               atcatcGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACg                                       >  2:36825/1‑148 (MQ=255)                                                                                                                      gCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAAc                                        >  4:285841/1‑140 (MQ=255)                                                                                                                      gCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAAc                                        <  3:285841/140‑1 (MQ=255)                                                                                                                           ctgATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTgg                           <  2:164529/148‑1 (MQ=255)                                                                                                                           ctgATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGc                          >  1:250232/1‑149 (MQ=255)                                                                                                                                tGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCgcgc                     >  2:88125/1‑149 (MQ=255)                                                                                                                                tGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACAAACAGGGCCTGGCCgcgc                     >  1:283108/1‑149 (MQ=255)                                                                                                                                  ccGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGATAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTg                   >  7:135076/1‑149 (MQ=255)                                                                                                                                    gTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGt                 <  3:300604/149‑1 (MQ=255)                                                                                                                                       tggttgcCCAGCACGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGa               >  4:1488/6‑148 (MQ=255)                                                                                                                                          ttgcCCAGCACGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTg                            >  8:43304/3‑132 (MQ=255)                                                                                                                                          ttgcCCAGCACGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTg                            <  7:43304/130‑1 (MQ=255)                                                                                                                                             cccAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGcc        <  7:134055/149‑1 (MQ=255)                                                                                                                                              ccAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCg       <  7:31899/149‑1 (MQ=255)                                                                                                                                                 gCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCt    <  1:41227/149‑1 (MQ=255)                                                                                                                                                 gCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCt    <  1:88125/149‑1 (MQ=255)                                                                                                                                                   gTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCACGCTGGTGAAAAAGCCGGCTCa  <  1:144592/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     ATGCAACTCAACCAGGGCCTGCAATTGTGACGACGAGGAGTATTGCCATGCAGCATCGACGTAACCGCCAACGCGGTTTCACCCTCATGGAAATCATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCA  >  NC_002947/1200319‑1200613

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