Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,622,932	A→C	11.9%	M269R (ATG→AGG) ‡	PP_4941 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,622,932	0	A	C	11.9%	153.3 / 13.6	59	M269R (ATG→AGG) ‡	PP_4941	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/21);  new base C (3/4);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

CAAGGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTG  >  NC_002947/5622784‑5623077                                                                                                                                                     |                                                                                                                                                  cAAGGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCa                                                                                                                                                   >  7:23411/1‑149 (MQ=255)   aGGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATg                                                                                                                                                 <  6:43330/149‑1 (MQ=255)     gCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGACCGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATgg                                                                                                                                                >  2:208522/1‑148 (MQ=255)          cAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCCAtgcaggca                                                                                                                                          <  8:103686/149‑1 (MQ=255)           aGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCCATgcaggcag                                                                                                                                         <  3:244146/149‑1 (MQ=255)           aGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGCTCGCCTTCACGTGAGCCGGCCGGCTGGTCGGGTGCTGCCTGCCATgcaggcag                                                                                                                                         <  6:14710/149‑1 (MQ=255)                   cagcagGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCAt                                                                                                                                                  >  2:35458/1‑132 (MQ=255)                   cagcagGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCAt                                                                                                                                                  <  1:35458/132‑1 (MQ=255)                      cagGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCaa                                                                                                                               >  2:121706/1‑148 (MQ=255)                        ggTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAg                                                                                                                            <  5:126690/149‑1 (MQ=255)                         gTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcag                                                                                                                                             >  8:110393/1‑131 (MQ=255)                         gTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcag                                                                                                                                             <  7:110393/131‑1 (MQ=255)                          ttGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTa                                                                                                                          <  5:181867/149‑1 (MQ=255)                             tgCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGcc                                                                                                                       <  1:121706/149‑1 (MQ=255)                              gCCAGGTGGATGCGCCATCGCGATCGCCCGATAACAAAAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGAGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGcc                                                                                                                       >  8:48746/1‑148 (MQ=255)                               ccAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCaaaa                                                                                                                             <  2:163382/141‑1 (MQ=255)                               ccAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAgagatcgg                                                                                                                     >  1:163382/1‑142 (MQ=255)                                        tgCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGc                                                                                                            <  7:125200/149‑1 (MQ=255)                                                  cgGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCAGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAAcgc                                                                                                       <  6:297952/144‑1 (MQ=255)                                                  cgGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCAGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAAcgc                                                                                                       >  5:297952/1‑144 (MQ=255)                                                         ccGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACg                                                                                           >  7:147004/1‑149 (MQ=255)                                                                 gggaacgggTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGccc                                                                                    >  2:226128/1‑148 (MQ=255)                                                                             aaCTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCCGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCacgac                                                                       >  5:76274/1‑149 (MQ=255)                                                                                  gCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCAGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGc                                                                  >  1:195321/1‑149 (MQ=255)                                                                                    ggAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTg                                                                >  1:133177/1‑149 (MQ=255)                                                                                     gAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGt                                                               >  7:202481/1‑149 (MQ=255)                                                                                     gAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGt                                                               >  4:240061/1‑149 (MQ=255)                                                                                        cacGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTcc                                                             >  2:282994/1‑148 (MQ=255)                                                                                        cacGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCCGCATGGCAGGCAGCGGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTccc                                                            <  6:76274/149‑1 (MQ=255)                                                                                         acGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCa                                                           <  8:23411/149‑1 (MQ=255)                                                                                           gcgGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAgg                                                         >  6:177654/1‑149 (MQ=255)                                                                                           gcgGAGATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGCCCCGCCGCTGCTGTCCCAgg                                                         >  7:225555/1‑149 (MQ=255)                                                                                             ggATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGATGATGCTGTCCCAGGCa                                                       >  6:255632/1‑149 (MQ=255)                                                                                                    ggCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGgc                                                >  7:225414/1‑149 (MQ=255)                                                                                                    ggCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGgc                                                >  8:267652/1‑149 (MQ=255)                                                                                                        aGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGt                                                                                                                >  1:303441/1‑81 (MQ=255)                                                                                                        aGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGt                                                                                                                <  2:303441/81‑1 (MQ=255)                                                                                                             ggCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCa                                       <  1:313119/149‑1 (MQ=255)                                                                                                                ttCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgc                                    >  3:223969/1‑149 (MQ=255)                                                                                                                   aCGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgcgcg                                 >  6:254853/1‑149 (MQ=255)                                                                                                                   aCGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgcgcg                                 >  8:172294/1‑149 (MQ=255)                                                                                                                     gTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGTATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGcgcgcg                               >  7:165869/1‑149 (MQ=255)                                                                                                                     gTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGcgcgcg                               <  6:102011/149‑1 (MQ=255)                                                                                                                     gTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGcgcgcg                               <  8:147004/149‑1 (MQ=255)                                                                                                                             gccggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctcggc                       >  7:296776/1‑149 (MQ=255)                                                                                                                                 gcgggTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCttt                   <  8:24919/146‑1 (MQ=255)                                                                                                                                 gcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAATGCGCCTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCttt                   >  7:287269/1‑149 (MQ=255)                                                                                                                                  cTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTg                  <  7:48746/149‑1 (MQ=255)                                                                                                                                     gTCGGGTGCTGCCTGCATGGCAGGCAGCTGAAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCTGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGcc                <  4:223969/148‑1 (MQ=255)                                                                                                                                      tCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAg              >  5:273595/1‑149 (MQ=255)                                                                                                                                       ctgcagctgcCTGCCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAgg             >  4:208806/6‑149 (MQ=255)                                                                                                                                       ctgcagctgcCTGCCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAgg             >  5:196985/6‑149 (MQ=255)                                                                                                                                         ggTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGcc           >  5:315659/1‑149 (MQ=255)                                                                                                                                             ctgcctgcCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTTCGGCACGCCGGCCCAACGCGGCAcgacga                                                                     <  1:191696/87‑1 (MQ=255)                                                                                                                                             ctgcctgcCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCAcgacga                                                                     >  2:191696/1‑87 (MQ=255)                                                                                                                                             ctgcctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATc       <  3:240061/149‑1 (MQ=255)                                                                                                                                               gcctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGc            >  1:14101/1‑142 (MQ=255)                                                                                                                                               gcctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGc            <  2:14101/142‑1 (MQ=255)                                                                                                                                                cctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctgc    <  1:282994/149‑1 (MQ=255)                                                                                                                                                  tgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctgctg  >  6:314297/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CAAGGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTG  >  NC_002947/5622784‑5623077

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