Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,002,830	G→T	10.0%	G30C (GGC→TGC)	PP_3531 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,002,830	0	G	T	10.0%	184.3 / 9.8	60	G30C (GGC→TGC)	PP_3531	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base T (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 = 2.27e-01

ACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGT  >  NC_002947/4002682‑4002970                                                                                                                                                     |                                                                                                                                             aCCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAg                                                                                                                                              <  8:66066/149‑1 (MQ=255)           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                                               <  4:92514/102‑1 (MQ=255)                                                                                                                                          gtaaacgcgAATGCCTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGAc                                               >  3:92514/6‑107 (MQ=255)                                                                                                                                          gtaaacgcgAATGCCTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGAc                                               <  4:87838/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                              >  6:166789/1‑122 (MQ=255)                                                                                                                                            gACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCg                              <  5:166789/122‑1 (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)                                                                                                                                                     |                                                                                                                                             ACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGT  >  NC_002947/4002682‑4002970

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