Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,971,955	C→G	8.6%	intergenic (‑390/+27)	PP_4379 ← / ← flgL	3‑oxoacyl‑ACP synthase/flagellar hook‑associated protein FlgL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,971,955	0	C	G	8.6%	178.2 / 6.8	58	intergenic (‑390/+27)	PP_4379/flgL	3‑oxoacyl‑ACP synthase/flagellar hook‑associated protein FlgL
Reads supporting (aligned to +/- strand):  ref base C (28/25);  new base G (3/2);  total (31/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 = 8.89e-01

ACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATG  >  NC_002947/4971811‑4972102                                                                                                                                                 |                                                                                                                                                      aCTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCgcagca                                                                                                                                                   <  7:93898/149‑1 (MQ=255)      aCCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGcc‑cc                                                                                                                                             <  4:14643/149‑1 (MQ=255)        cccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCTGCGCAGcc‑ctat                                                                                                                                           >  5:26583/1‑146 (MQ=255)        cccGAATGATCAGCTTCGCAGCAGATATCGAAGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAg                                                                                                                                           >  6:316049/1‑149 (MQ=255)         ccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCTGCGCAGcc‑ctatc                                                                                                                                          <  6:26583/149‑5 (MQ=255)         ccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAg                                                                                                                                           >  3:313829/1‑148 (MQ=255)           gAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGAt                                                                                                                                        >  3:20378/1‑149 (MQ=255)                 tCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATcaggc                                                                                                                                   <  4:313829/148‑1 (MQ=255)                     cTTCGCAGCAGATATCGATGGGGCCGCATCGAGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATcaggcagg                                                                                                                                >  1:33877/1‑147 (MQ=255)                              atatatCGATGGGGCCGCATCGGGGCCCTTTCGCGACGCAATGCCGGTCCCAAGAATACGCCACCTCCAGGTGAGTGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCATCCTCCAGGATCAGGCAGGCAGATTAGa                                                                                                                       <  7:203046/146‑1 (MQ=255)                               gATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTa                                                                                                                         >  5:169962/1‑144 (MQ=255)                               gATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTa                                                                                                                         <  6:169962/144‑1 (MQ=255)                               gATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCt                                                                                                                    <  2:142748/149‑1 (MQ=255)                                  atCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGt                                                                                                                 >  5:106378/1‑149 (MQ=255)                                   gcgaaagggCCGCGATGCGGCCCTTTCGCGACGCAAGGCCCCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGCGGCCCTGGTGTCGCGATAGGGCTAAGCAGCAGCC‑CCACGATCAGGCAGGCAGATTAGATCTTGTc                                                                                                                >  1:110053/7‑149 (MQ=255)                                       tGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAAc                                                                                                            <  6:88387/149‑1 (MQ=255)                                        ggggCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGAtt                                                                                                                          <  6:222993/134‑1 (MQ=255)                                        ggggCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGAtt                                                                                                                          >  5:222993/1‑134 (MQ=255)                                         gggCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACa                                                                                                           <  5:285341/148‑1 (MQ=255)                                            ccGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATcaggc                                                                                                                                   >  7:35649/1‑121 (MQ=255)                                            ccGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATcaggc                                                                                                                                   <  8:35649/121‑1 (MQ=255)                                             cGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCtt                                                                                                      <  5:316049/149‑1 (MQ=255)                                                    cgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCt                                                                                                  >  5:213826/1‑146 (MQ=255)                                                                gACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCa                                                                                    >  5:300095/1‑148 (MQ=255)                                                                    cAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGc                                                                               <  2:33877/149‑1 (MQ=255)                                                                        gCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGAAGATTAGATCTTGTGGAACAGGCTTAGCTGAGAAATGCGGGCAAACGCCAGt                                                                           >  3:57100/1‑149 (MQ=255)                                                                         ccGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGtt                                                                          <  1:122849/149‑1 (MQ=255)                                                                            cTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCg                                                                       <  1:115118/149‑1 (MQ=255)                                                                             tCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCg                                                                       >  3:72093/1‑148 (MQ=255)                                                                                cACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAAc                                                                                 <  6:234287/135‑1 (MQ=255)                                                                                cACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCA‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAAc                                                                                 >  5:234287/1‑135 (MQ=255)                                                                                 aCGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCt                                                                  >  6:234116/1‑149 (MQ=255)                                                                                     aTACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTc                                                                >  1:131863/1‑147 (MQ=255)                                                                                     aTACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCaa                                                              <  4:72093/149‑1 (MQ=255)                                                                                         gCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCAt                                                         >  4:289788/1‑149 (MQ=255)                                                                                          ccACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATg                                                        <  3:289788/149‑1 (MQ=255)                                                                                          ccACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATg                                                        <  6:300095/149‑1 (MQ=255)                                                                                               tCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGc                                                      >  7:4989/1‑146 (MQ=255)                                                                                               tCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGc                                                      <  8:4989/146‑1 (MQ=255)                                                                                                           gacgacCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGaa                                                                     >  6:3624/1‑120 (MQ=255)                                                                                                           gacgacCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAg                                       >  6:122158/1‑149 (MQ=255)                                                                                                           gacgacCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAg                                       >  3:274916/1‑149 (MQ=255)                                                                                                           gacgacCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAg                                       >  7:27880/1‑149 (MQ=255)                                                                                                             cgacCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCt                                     <  5:144435/149‑1 (MQ=255)                                                                                                              gacCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCt                                     >  8:252986/1‑148 (MQ=255)                                                                                                                    tgGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAAtggt                              <  4:20378/149‑1 (MQ=255)                                                                                                                      ggggCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAAtggtgg                            >  6:7839/1‑149 (MQ=255)                                                                                                                               tgtCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGcc                    <  7:252986/148‑1 (MQ=255)                                                                                                                               tgtCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCa                   <  8:27880/149‑1 (MQ=255)                                                                                                                                 tCGCGATAGGGCTGCGCAGCATCC‑CCAGGATCAGGCAGTCAGTTTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGGGAGTTGCGTAGCTTCAACGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCAt                  <  1:229173/149‑1 (MQ=255)                                                                                                                                 tCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATg                 >  8:97593/1‑149 (MQ=255)                                                                                                                                 tCGCGATAGGGATGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGATTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATg                 >  8:78073/1‑149 (MQ=255)                                                                                                                                   gcgATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGt                >  8:59455/1‑148 (MQ=255)                                                                                                                                    cgATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCa              <  5:7839/149‑1 (MQ=255)                                                                                                                                     gATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAg             >  6:92216/1‑149 (MQ=255)                                                                                                                                      ctggggCTGCTGCGCAGCC‑CTAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCaa                                                              >  7:245640/4‑100 (MQ=255)                                                                                                                                      ctggggCTGCTGCGCAGCC‑CTAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCaa                                                              <  8:245640/97‑1 (MQ=255)                                                                                                                                      ctggggCTGCTGCGCAGCC‑CTAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAg             >  5:137354/4‑148 (MQ=255)                                                                                                                                       tAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTa           <  6:213826/149‑1 (MQ=255)                                                                                                                                               cgcAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAAt   >  5:22378/1‑149 (MQ=255)                                                                                                                                                gcagcaGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATg  <  2:110053/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                      ACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCC‑CCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAA‑GCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTATCGCCAATG  >  NC_002947/4971811‑4972102

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