Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,002,833	A→C	9.8%	I31L (ATT→CTT)	PP_3531 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,002,833	0	A	C	9.8%	180.7 / 9.0	60	I31L (ATT→CTT)	PP_3531	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base C (3/3);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCG  >  NC_002947/4002692‑4002981                                                                                                                                              |                                                                                                                                                     cTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAATGCCTTcgcgg                                                                                                                                               <  7:312912/149‑2 (MQ=255)  tttCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGtt                                                                                                                                              >  3:314701/1‑149 (MQ=255)        tttGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTAc                                                                                                                                           <  7:181964/146‑1 (MQ=255)        tttGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGt                                                                                                                                         >  2:118124/1‑148 (MQ=255)           gCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGag                                                                                                                                      <  3:43461/148‑1 (MQ=255)           gCCGCGCATACTGAAACACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGag                                                                                                                                      >  2:309208/1‑148 (MQ=255)                 cATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCa                                                                                                                               <  4:314701/149‑1 (MQ=255)                 cATACTGAAATACCCACCCCTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCa                                                                                                                               <  5:175584/149‑1 (MQ=255)                     cTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGgc                                                                                                                            <  1:66955/148‑1 (MQ=255)                      tGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGgcgc                                                                                                                          >  8:10732/1‑149 (MQ=255)                      tGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGgcgc                                                                                                                          >  3:129486/1‑149 (MQ=255)                        aaaTACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGgcgc                                                                                                                          <  2:254759/147‑1 (MQ=255)                        aaaTACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGgcgc                                                                                                                          >  1:254759/1‑147 (MQ=255)                         aaTACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTc                                                                                                                       >  3:168853/1‑149 (MQ=255)                              ccacccTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGcc                                                                                                                         >  8:10099/1‑142 (MQ=255)                              ccacccTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGcc                                                                                                                         <  7:10099/142‑1 (MQ=255)                               cacccTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTa                                                                                                                 >  5:260041/1‑149 (MQ=255)                               cacccTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTa                                                                                                                 <  1:309208/149‑1 (MQ=255)                               cacccTTGCAGGGTACCGGCCATGAACCGCTGCCTGCGGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTa                                                                                                                 >  6:25131/1‑149 (MQ=255)                                              ccGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCgtgg                                                                                                  <  4:168853/149‑1 (MQ=255)                                              ccGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCgtg                                                                                                   <  2:304852/148‑1 (MQ=255)                                                 gCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGgtgt                                                                                               >  3:180177/1‑149 (MQ=255)                                                       aaCCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCAc                                                                                         <  4:129486/149‑1 (MQ=255)                                                         ccGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCa                                                                                       <  4:180177/149‑1 (MQ=255)                                                         ccGCTACCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCa                                                                                       >  2:118578/1‑149 (MQ=255)                                                          cGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCAt                                                                                      >  7:53309/1‑149 (MQ=255)                                                          cGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCAt                                                                                      <  4:198450/149‑1 (MQ=255)                                                              gcctgcTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGa                                                                                    >  8:178366/1‑147 (MQ=255)                                                               cctgcTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGaaa                                                                                  <  7:178366/148‑1 (MQ=255)                                                                   cagctCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCa                                                                              <  5:25131/146‑1 (MQ=255)                                                                      ctCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGt                                                                          >  7:195055/1‑149 (MQ=255)                                                                                 gctgctGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCa                                                                                             >  2:317090/1‑119 (MQ=255)                                                                                 gctgctGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCa                                                                                             <  1:317090/119‑1 (MQ=255)                                                                                        tgCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGAccc                                                        >  6:41841/1‑148 (MQ=255)                                                                                             ctctGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCgtg                                                                                                   >  6:303462/1‑101 (MQ=255)                                                                                             ctctGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCgtg                                                                                                   <  5:303462/101‑1 (MQ=255)                                                                                             ctctGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAg                                                    >  3:264123/1‑148 (MQ=255)                                                                                                  cccTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGAt                                                 <  6:260041/146‑1 (MQ=255)                                                                                                  cccTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAAc                                              <  8:53309/149‑1 (MQ=255)                                                                                                  cccTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAAc                                              >  5:228593/1‑149 (MQ=255)                                                                                                       gccgccGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTg                                           <  1:200787/147‑1 (MQ=255)                                                                                                          gccgATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGt                                      >  3:261157/1‑149 (MQ=255)                                                                                                          gccgATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCCGCAGGTTATCCTGCAGGt                                      <  5:41841/149‑1 (MQ=255)                                                                                                             gATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGc                                   >  4:305350/1‑149 (MQ=255)                                                                                                                   tGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCCGCACCAGGCAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATCACCTGCGGGTGGCCTGc                               >  8:151921/1‑147 (MQ=255)                                                                                                                              cTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACg                  >  7:103111/1‑149 (MQ=255)                                                                                                                               cgtaaacgcgAATGCCTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGc                 >  4:235650/7‑149 (MQ=255)                                                                                                                                gtaaacgcgAATGCCTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGAc                                                          >  3:87838/6‑107 (MQ=255)                                                                                                                                gtaaacgcgAATGCCTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGAc                                                          >  3:92514/6‑107 (MQ=255)                                                                                                                                gtaaacgcgAATGCCTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGAc                                                          <  4:87838/102‑1 (MQ=255)                                                                                                                                gtaaacgcgAATGCCTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGAc                                                          <  4:92514/102‑1 (MQ=255)                                                                                                                                aCGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGTGTGGCCTGCAAGTGGCTGTACGcc                >  8:147699/1‑149 (MQ=255)                                                                                                                                 cGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCt               <  7:160188/149‑1 (MQ=255)                                                                                                                                 cGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCt               <  8:195055/149‑1 (MQ=255)                                                                                                                                  gACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCg                                         <  5:166789/122‑1 (MQ=255)                                                                                                                                  gACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCg                                         >  6:166789/1‑122 (MQ=255)                                                                                                                                  gACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCt                                 >  2:235340/1‑130 (MQ=255)                                                                                                                                  gACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCt                                 <  1:235340/130‑1 (MQ=255)                                                                                                                                   aCCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCtgt             <  3:305350/149‑1 (MQ=255)                                                                                                                                         aaGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTg                                     <  2:268767/119‑1 (MQ=255)                                                                                                                                              aTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCg  >  8:267825/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                     CTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCG  >  NC_002947/4002692‑4002981

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