Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,247,267	G→T	6.3%	A465E (GCG→GAG)	PP_4623 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,247,267	0	G	T	6.3%	200.0 / 2.9	63	A465E (GCG→GAG)	PP_4623	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base T (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CAGATCCTTCAAGGGCGTTTGCGCACCATTGAGGCAAACCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATAGAGGTACTCCCCAACAAGGAACGGTG  >  NC_002947/5247119‑5247412                                                                                                                                                     |                                                                                                                                                  cagATCCTTCAAGGGCGTTTGCGCACCATTGAGGCAAACCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGgccagcg                                                                                                                                                   >  7:352188/1‑149 (MQ=255)  agATCCTTCAAGGGCGTTTGCGCACCATAGAGGCAAACCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGCTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGgccagcgc                                                                                                                                                  <  3:17320/149‑1 (MQ=255)   gctCCTTCAAGGGCGTTTGCGCACCATTGAGGCAAACCAGACCGGCGCAGATGACATTGGCCTGCATTTTCCCCACGCGTCCGTGGTGCAACAGATCGAACGCCGCCGCGTGAAAACCCTGATGACGAGCGGCCTGATTGgccagcgac                                                                                                                                                 >  8:321409/3‑147 (MQ=255)   gATCCTTCAAGGGCGTTTGCGCACCATTGAGGCAAACCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGgccagcgc                                                                                                                                                  >  1:185507/1‑148 (MQ=255)       cTTCAAGGGCGTTTGCGCACCATTGAGGCAAACCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCCGCCTGATTGGCCAGCGCCAGcc                                                                                                                                             <  8:79679/149‑1 (MQ=255)          cAAGGGCGTTTGCGCACCATTGAGGCAAACCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGt                                                                                                                                           <  8:69463/148‑1 (MQ=255)           aaGGGCGTTTGCGCACCATTGAGGCAAACCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGc                                                                                                                                         <  2:185507/149‑1 (MQ=255)           aaGGGCGTTTGCGCACCATTGAGGC‑AACCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGc                                                                                                                                         <  7:155181/148‑1 (MQ=255)                           cATTGAGGCAAACCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTcaca                                                                                                                          >  4:303091/1‑148 (MQ=255)                                aGGCAAACCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGc                                                                                                                    >  5:3013/1‑149 (MQ=255)                                 ggCAAACCACACCGGCGCTGATGACATTGGACTGCATTTTCCCCAGGCGTGCGTGGTGAAACAGATAGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCg                                                                                                                   >  1:405085/1‑149 (MQ=255)                                  gCAAACCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCgg                                                                                                                  >  2:27714/1‑149 (MQ=255)                                    aaaCCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGcc                                                                                                                <  8:352188/149‑1 (MQ=255)                                    aaaCCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGc                                                                                                                 >  7:391529/1‑148 (MQ=255)                                                  cTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTgtc                                                                                                  <  3:303091/149‑1 (MQ=255)                                                  cTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTgtc                                                                                                  <  6:3013/149‑1 (MQ=255)                                                   tgatgaCATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTgtcg                                                                                                 >  4:202052/1‑149 (MQ=255)                                                       gaCATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTTTCGTCa                                                                                              >  4:341612/1‑148 (MQ=255)                                                         cATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGgccagcgc                                                                                                                                                  <  6:385112/94‑1 (MQ=255)                                                         cATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGgccagcgc                                                                                                                                                  >  5:385112/1‑94 (MQ=255)                                                         cATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGt                                                                                                           <  5:230293/133‑1 (MQ=255)                                                         cATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGt                                                                                                           >  6:230293/1‑133 (MQ=255)                                                            tGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGc                                                                                        >  6:214337/1‑149 (MQ=255)                                                               ccTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTa                                                                                     <  5:181716/149‑1 (MQ=255)                                                                cTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTAc                                                                                    >  7:369632/1‑149 (MQ=255)                                                                 tGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACg                                                                                   >  4:358323/1‑149 (MQ=255)                                                                  gCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCCGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACg                                                                                   >  7:265106/1‑148 (MQ=255)                                                                        tCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTg                                                                             <  3:202052/148‑1 (MQ=255)                                                                        tCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGc                                                                            >  5:380187/1‑149 (MQ=255)                                                                          cccAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCtt                                                                          <  6:380187/149‑1 (MQ=255)                                                                           ccAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATaga                                                                                          >  1:290066/1‑132 (MQ=255)                                                                           ccAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATaga                                                                                          <  2:290066/132‑1 (MQ=255)                                                                           ccAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCtt                                                                          >  6:399244/1‑148 (MQ=255)                                                                           ccAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTg                                                                         <  7:164921/149‑1 (MQ=255)                                                                             aGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGAc                                                                       <  8:391529/149‑1 (MQ=255)                                                                               gcgtgcgtGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCt                                                                              >  1:206824/1‑140 (MQ=255)                                                                               gcgtgcgtGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCt                                                                              <  2:206824/140‑1 (MQ=255)                                                                                cgtgcgtgGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGCCGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTa                                                                    <  6:120454/149‑1 (MQ=255)                                                                                 gtgcgt‑gTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGc                                                                                                                                  <  5:218476/85‑1 (MQ=255)                                                                                 gtgcgt‑gTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGc                                                                                                                                  >  6:218476/1‑85 (MQ=255)                                                                                   gcgtgGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCaaa                                                                                                       <  7:180466/111‑1 (MQ=255)                                                                                   gcgtgGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCaaa                                                                                                       >  8:180466/1‑111 (MQ=255)                                                                                     gtggtgCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGGTTGGCCAGCTCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCa                                                               >  6:80877/1‑149 (MQ=255)                                                                                     gtggtgCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCTCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCa                                                               >  2:228542/1‑149 (MQ=255)                                                                                           cAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGa                                                                                 <  8:280127/125‑1 (MQ=255)                                                                                           cAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGa                                                                                 >  7:280127/1‑125 (MQ=255)                                                                                                 aTCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATc                                                     >  6:377229/1‑147 (MQ=255)                                                                                                     aaCGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGttt                                               >  3:368189/1‑149 (MQ=255)                                                                                                      aCGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGtttt                                              >  2:78293/1‑149 (MQ=255)                                                                                                       cgccgcCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGttttt                                             <  5:214337/149‑1 (MQ=255)                                                                                                          cgccgGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGttttt                                             >  3:205441/1‑146 (MQ=255)                                                                                                          cgccgGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGttttt                                             <  4:205441/146‑1 (MQ=255)                                                                                                            ccgGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATc                                        <  2:131970/149‑1 (MQ=255)                                                                                                             cgGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCa                                       <  3:341612/149‑1 (MQ=255)                                                                                                               ggTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACg                                     >  5:422674/1‑149 (MQ=255)                                                                                                                  gAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGctc                                  >  3:334051/1‑149 (MQ=255)                                                                                                                    aaaCCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTg                                >  4:206864/1‑149 (MQ=255)                                                                                                                      acgcTGTTGACGAGCGGCCTGATTGGCGAGCGGCAGCCGTGCATCAAGCACCTGACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACTGATCTGCTTGAGGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGa                              <  7:321409/146‑1 (MQ=255)                                                                                                                      aCCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGa                              <  1:27714/149‑1 (MQ=255)                                                                                                                        ccTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGa                              >  4:295646/1‑147 (MQ=255)                                                                                                                           gatgaCGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCg                                           >  5:209783/1‑131 (MQ=255)                                                                                                                           gatgaCGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCg                                           <  6:209783/131‑1 (MQ=255)                                                                                                                            atgaCGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGTCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCTAATCGAGTTTTTCGATCACGCTCTGGATagag                         <  2:405085/148‑1 (MQ=255)                                                                                                                            atgaCGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATagag                         <  7:77145/148‑1 (MQ=255)                                                                                                                                           gttGGCCAGCTCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATAGAGGTACTCCCCAACAAg          <  5:80877/147‑1 (MQ=255)                                                                                                                                           aTTGGCCAGCTCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATAGAGGTACTCCCCAACAAg          <  1:228542/148‑1 (MQ=255)                                                                                                                                           aTTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATAGAGGTACTCCCCAACAAg          >  6:36603/1‑148 (MQ=255)                                                                                                                                             tGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCTGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACTCTCTGGATAGATGTACTCCCCAACAAGGa        <  8:369632/148‑1 (MQ=255)                                                                                                                                              ggCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATAGAGGTACTCCCCAACAAGGAAc      >  2:395168/1‑149 (MQ=255)                                                                                                                                                 cagcgccagcCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATagag                         <  5:75867/127‑1 (MQ=255)                                                                                                                                                 cagcgccagcCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATagag                         <  5:81739/127‑1 (MQ=255)                                                                                                                                                 cagcgccagcCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATagag                         >  6:75867/1‑127 (MQ=255)                                                                                                                                                 cagcgccagcCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATagag                         >  6:81739/1‑127 (MQ=255)                                                                                                                                                 cagcgccagcCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATAGAGGTACTCCCCAACAAGGAACGGt   >  5:424513/1‑149 (MQ=255)                                                                                                                                                  agcgccagcCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATAGAGGTACTCCCCAACAAGGAACGGTg  <  1:395168/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CAGATCCTTCAAGGGCGTTTGCGCACCATTGAGGCAAACCAGACCGGCGCTGATGACATTGGCCTGCATTTTCCCCAGGCGTGCGTGGTGCAACAGATCGAACGCCGCCGGGTGAAAACCCTGATGACGAGCGGCCTGATTGGCCAGCGCCAGCCGTGCATCAAGCACCTCACACAGTGCGGCCATGGTCAAAGTGTCGTCATAGAGCGTACGGATCTGCTTGACGTAGTTCACGCGCAAATCCAGTTTTTCGATCACGCTCTGGATAGAGGTACTCCCCAACAAGGAACGGTG  >  NC_002947/5247119‑5247412

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