Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,065,216	A→G	6.3%	L210P (CTG→CCG)	PP_0923 ←	acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,065,216	0	A	G	6.3%	183.4 / 2.2	63	L210P (CTG→CCG)	PP_0923	acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (34/25);  new base G (2/2);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

CCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACGCCACGGTCGATGGCTCCTGCAAGCAG  >  NC_002947/1065070‑1065355                                                                                                                                                   |                                                                                                                                            ccTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAgg                                                                                                                                           <  2:293006/149‑1 (MQ=255)     ccGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGc                                                                                                                                         <  4:175703/147‑1 (MQ=255)       gACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGaa                                                                                                                                      >  6:149564/1‑148 (MQ=255)         cGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAg                                                                                                                                   <  6:86809/149‑1 (MQ=255)          gTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAgg                                                                                                                                  <  5:59751/149‑1 (MQ=255)          gTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAgg                                                                                                                                  <  6:291634/149‑1 (MQ=255)          gTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCCGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAgg                                                                                                                                  <  3:305559/149‑1 (MQ=255)             gcgcCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGTACAGgcgc                                                                                                                               <  4:332023/149‑1 (MQ=255)                      ccACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTgg                                                                                                                      >  5:104577/1‑149 (MQ=255)                      ccACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTgg                                                                                                                      >  7:148412/1‑149 (MQ=255)                       cACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGGACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTgg                                                                                                                      <  1:319184/148‑1 (MQ=255)                       cACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAggg                                                                                                                                          >  3:240137/1‑128 (MQ=255)                       cACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAggg                                                                                                                                          <  4:240137/128‑1 (MQ=255)                          gTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCCCGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCCTAAGGTgggccc                                                                                                                  >  8:351313/1‑146 (MQ=255)                             ttCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGacac                                                                                                                 >  3:262979/1‑147 (MQ=255)                             ttCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGacac                                                                                                                 >  6:78772/1‑147 (MQ=255)                                   cccGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCCCGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGgtcgt                                                                                                         >  5:236753/1‑149 (MQ=255)                                   cccGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGgtcgt                                                                                                         >  8:3792/1‑149 (MQ=255)                                            cgcCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCaa                                                                                                <  5:149564/149‑1 (MQ=255)                                                  gTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACg                                                                                          >  5:373107/1‑149 (MQ=255)                                                    gAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGgc                                                                                        >  8:274442/1‑149 (MQ=255)                                                         gCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTAc                                                                                   >  3:85951/1‑149 (MQ=255)                                                          ccTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTAcc                                                                                  <  4:85951/149‑1 (MQ=255)                                                              ttGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCgg                                                                                <  7:3792/147‑1 (MQ=255)                                                                  gcCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGAt                                                                          >  1:161594/1‑149 (MQ=255)                                                                  gcCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGAt                                                                          >  2:317558/1‑149 (MQ=255)                                                                   ccTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATc                                                                         >  5:279063/1‑149 (MQ=255)                                                                    cTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCt                                                                        >  7:46759/1‑149 (MQ=255)                                                                               gCACCGCACGTTCCTTGCCCTCGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGc                                                             <  7:38506/149‑1 (MQ=255)                                                                                cACCGCACGTTCCTTGCACACGCTACCTATTGGCTGCAGCAGGTGGATGCACAACTCGCACGGCGGCAGGGAGAACAGGCGCATAAGGTGTGACACCACGTCGTCCTCACCAATGAACGCCGCTACCGGGTCGATCAGGCAATCaag                                                              >  8:147399/1‑145 (MQ=255)                                                                                 aCCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAg                                                           <  6:373107/149‑1 (MQ=255)                                                                                 aCCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAg                                                           <  2:161594/149‑1 (MQ=255)                                                                                   cGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAgg                                                          >  2:230872/1‑148 (MQ=255)                                                                                    gCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTa                                                        >  5:266457/1‑149 (MQ=255)                                                                                     cACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTAc                                                       <  6:266457/149‑1 (MQ=255)                                                                                        gTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTg                                                     >  7:120950/1‑148 (MQ=255)                                                                                             ttGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGc                                               >  7:287447/1‑149 (MQ=255)                                                                                               gCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGc                                               <  7:274442/147‑1 (MQ=255)                                                                                                   aCGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTgg                                         <  5:78772/149‑1 (MQ=255)                                                                                                       tACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGCCACCAGGTCGTCCTCACCAATGCACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTg                                      >  2:295835/1‑148 (MQ=255)                                                                                                             ttGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGac                               <  6:236753/149‑1 (MQ=255)                                                                                                               ggCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGaca                              <  8:148412/148‑1 (MQ=255)                                                                                                                  tGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACGcc                          >  3:13239/1‑149 (MQ=255)                                                                                                                           tGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACTGCGACTGGCTGCACTGACACGCCACGGTCGAt                 <  1:230872/149‑1 (MQ=255)                                                                                                                            ggATGCACACCTCGCCCCGCGGCAGGGCGAACATGCGCATAAGGTCGGACACCAGTTCGTCCTCACCAATGAACCGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACACCGACTGGCTGCACTGACACGCCACGCTCGATg                >  6:291305/1‑149 (MQ=255)                                                                                                                             gATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCACGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACGCCACGGTCGATgg               >  8:86157/1‑149 (MQ=255)                                                                                                                              aTGCACACCTCGCCCCACGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACGCCACGGTCGATGGc              >  1:94954/1‑149 (MQ=255)                                                                                                                               tGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACg                            <  6:247784/134‑1 (MQ=255)                                                                                                                               tGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACg                            >  5:247784/1‑134 (MQ=255)                                                                                                                                gCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGCGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACGCCACGGTCGATGGCTc            >  1:37001/1‑149 (MQ=255)                                                                                                                                  acacCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACt                                  >  6:218892/1‑125 (MQ=255)                                                                                                                                  acacCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACt                                  <  5:218892/125‑1 (MQ=255)                                                                                                                                   ctgttCGCCCTGCCGCGGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACGCCACGGTc                    >  6:76849/5‑138 (MQ=255)                                                                                                                                   ctgttCGCCCTGCCGCGGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACGCCACGGTc                    <  5:76849/134‑1 (MQ=255)                                                                                                                                   ctgttCGCCCTGCCGCGGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACAGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACGCCACGGTCGATg                <  6:326622/138‑1 (MQ=255)                                                                                                                                   ctgttCGCCCTGCCGCGGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACAGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACGCCACGGTCGATg                >  5:326622/5‑142 (MQ=255)                                                                                                                                   cacCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCCCCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGCCACGCCAAGTTCGATGGCTCCTg         >  7:339888/1‑149 (MQ=255)                                                                                                                                     ccTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCAc                                                                                                   <  4:263866/57‑1 (MQ=255)                                                                                                                                     ccTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCAc                                                                                                   >  3:263866/1‑57 (MQ=255)                                                                                                                                       tCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACGCCACGGTCGATGGCTCCTGCAAg     >  4:351474/1‑149 (MQ=255)                                                                                                                                        cGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGc                                                             >  2:294099/1‑92 (MQ=255)                                                                                                                                        cGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGc                                                             <  1:294099/92‑1 (MQ=255)                                                                                                                                          ccccGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGCCACGCCACGGTCGATGGCTCCTGCAAGccg  >  7:335993/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                            CCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACACCTCGCCCCGCGGCAGGGCGAACAGGCGCATAAGGTGGGACACCAGGTCGTCCTCACCAATGAACGGCGCTACCGGGTCGATCTGGCCATCACGCAGGTACTGAACAGCGACTGGCTGCACTGACACGCCACGGTCGATGGCTCCTGCAAGCAG  >  NC_002947/1065070‑1065355

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