Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,108,703	G→C	7.8%	A312G (GCG→GGG) ‡	PP_5359 ←	CobW/P47K family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,108,703	0	G	C	7.8%	186.4 / 2.7	60	A312G (GCG→GGG) ‡	PP_5359	CobW/P47K family protein
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base C (3/2);  total (31/29)
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

GCAGATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTGGCCCGTCG  >  NC_002947/6108557‑6108848                                                                                                                                                    |                                                                                                                                                     gCAGATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAAcgcg                                                                                                                                                    >  1:166890/1‑149 (MQ=255)      tCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGggttg                                                                                                                                               <  1:131221/149‑1 (MQ=255)              gcgcCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGggtt                                                                                                                                                <  1:7208/140‑1 (MQ=255)              gcgcCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGggtt                                                                                                                                                >  2:7208/1‑140 (MQ=255)              gcgcCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATGCAACCCGCGTTGCAAGCCGggttgaatgaat                                                                                                                                        <  3:320689/148‑8 (MQ=255)                 cccTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCt                                                                                                                                 <  6:132414/149‑1 (MQ=255)                 cccTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGCGTTGCAAGCCGggttgaatgaa                                                                                                                                         >  6:169877/1‑138 (MQ=255)                   cTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGa                                                                                                                               <  3:85213/149‑1 (MQ=255)                      cccGTTCGAAATCACGACCCCACCACGCGCCGACAACGTCATCGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATg                                                                                                                             <  8:273158/148‑1 (MQ=255)                       ccGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATgg                                                                                                                            >  3:237990/1‑148 (MQ=255)                        cGTTCGTAATCAAGACCCCACCACGCGCCGAAAACGTCATAGAGATGAAGCAATGGCGACAGGATCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCAATCTTGCAACGCggg                                                                                                                                                  <  6:317389/128‑1 (MQ=255)                        cGTTCGAAATCAAGACACCACCACGCGACGACAACGTCATAGAGATTCAGCAATGGCGCCAGGAGCATAACCGCTGGCGAAACTGATCAACT‑GACGCGACATTCATTCAACCCAGCTTGCAACGCggg                                                                                                                                                  >  5:317389/1‑128 (MQ=255)                           tCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAAt                                                                                                                       <  3:8913/149‑1 (MQ=255)                             gAATTCATGACCCCACCCGGCGCCTAAAACGTCATTGAGATGCAGCATTGGCGCCAGGAGCATAACCGCTGGCGAAAATTATCAAATATACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAAt                                                                                                                       <  1:103346/148‑1 (MQ=255)                             gAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAAt                                                                                                                       >  4:8913/1‑147 (MQ=255)                             gAAATCAAGACAACACCACGATCCGACAAAGTCATAGATATGAAGCAATGGCGCCAGGAGCATAAAAGCTGGCGCAAATGATCAACT‑GACTCGAAATTCATTCAAACCGGATTTAAAAGCGGGTTG‑‑A‑GGTTGCGTCTGATGTAAAATCa                                                                                                                     >  3:348853/1‑149 (MQ=255)                                 tCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTc                                                                                                                 <  7:278651/149‑1 (MQ=255)                                     gACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATg                                                                                                                             <  8:191502/133‑1 (MQ=255)                                       ccccACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACg                                                                                                           >  3:23521/1‑149 (MQ=255)                                           accacGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTg                                                                                                       <  4:237990/149‑1 (MQ=255)                                             cacGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTc                                                                                                     >  4:187065/1‑149 (MQ=255)                                              acGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTcc                                                                                                    <  2:166890/149‑1 (MQ=255)                                                gcgcCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCt                                                                                                  >  1:72302/1‑149 (MQ=255)                                                gcgcCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCt                                                                                                  >  2:166619/1‑149 (MQ=255)                                                  gcCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTgc                                                                                                >  8:50280/1‑149 (MQ=255)                                                        aaCGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGt                                                                                          >  1:175881/1‑149 (MQ=255)                                                              aTAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCt                                                                                                  <  5:297761/135‑1 (MQ=255)                                                              aTAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCt                                                                                                  >  6:297761/1‑135 (MQ=255)                                                              aTAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTgg                                                                                    >  4:75730/1‑149 (MQ=255)                                                                       cagcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GCCGCGACATTCATTCACCCCGGCTTGCACCGCGGGTTG‑‑C‑GGTTGCGGGTGAGGGAAAAACCGCGCGATCCGGGTGGCCCTGCGCCCGGCGCTGGGCGGccccc                                                                           >  8:54772/1‑146 (MQ=255)                                                                         gcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCg                                                                         >  4:21763/1‑149 (MQ=255)                                                                            aatgcgcCATGATCAAAAACCCTTGCGCAAATTTTCAAAT‑TACGCGAACTTTATTCAACCCGGGCTGCAACGGGGGGTGTTC‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCgg                                                                        <  1:293571/146‑1 (MQ=255)                                                                            aTGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGc                                                                      <  7:50280/149‑1 (MQ=255)                                                                                       aGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCAtt                                                           <  2:175881/149‑1 (MQ=255)                                                                                            aaCCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCa                                                             <  6:333352/142‑1 (MQ=255)                                                                                            aaCCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCa                                                             >  5:333352/1‑142 (MQ=255)                                                                                            aaCCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGcc                                                      >  3:231260/1‑149 (MQ=255)                                                                                            aaCCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGc                                                       >  7:226360/1‑148 (MQ=255)                                                                                                    gcgcAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGaa                                               <  4:23521/148‑1 (MQ=255)                                                                                                    gcgcAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAAt                                              <  3:187065/149‑1 (MQ=255)                                                                                                    gcgcAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAAt                                              <  2:72302/149‑1 (MQ=255)                                                                                                     cgcAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATg                                             <  4:231260/149‑1 (MQ=255)                                                                                                     cgcAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATg                                             <  8:226360/149‑1 (MQ=255)                                                                                                        aaCTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑CGGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGac                                           >  1:93626/1‑149 (MQ=255)                                                                                                          cTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGaccac                                        >  5:343778/1‑149 (MQ=255)                                                                                                             aTCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCggg                                                                       <  4:183210/115‑1 (MQ=255)                                                                                                             aTCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCggg                                                                       >  3:183210/1‑115 (MQ=255)                                                                                                                 aCT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACAACCAGCCCt                                 >  3:197268/1‑149 (MQ=255)                                                                                                                  cT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTc                                >  5:399/1‑149 (MQ=255)                                                                                                                      aCGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCAt                             <  6:343778/149‑1 (MQ=255)                                                                                                                          gACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCAATAGAATGACCACCAGCCCTCCATGcag                         >  3:228977/1‑147 (MQ=255)                                                                                                                              ttcattcaACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGa                     <  3:75730/149‑1 (MQ=255)                                                                                                                                 ccgcaACCCGCGTTGCAAGCCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATaa                  >  4:194095/4‑149 (MQ=255)                                                                                                                                 ccgcaACCCGCGTTGCAAGCCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATaa                  >  3:83436/4‑149 (MQ=255)                                                                                                                                 attcaACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATaa                  >  2:276503/1‑149 (MQ=255)                                                                                                                                  ttcaACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCtgt                        <  7:16112/142‑1 (MQ=255)                                                                                                                                  ttcaACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCtgt                        >  8:16112/1‑142 (MQ=255)                                                                                                                                       cccGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCAATGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCt            >  1:119377/1‑149 (MQ=255)                                                                                                                                         cgcgttGCAAGCCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTg          <  4:83436/145‑1 (MQ=255)                                                                                                                                         cGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTg          <  2:93626/149‑1 (MQ=255)                                                                                                                                                 aaCGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTGGCCCGTCg  >  1:362238/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GCAGATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTGGCCCGTCG  >  NC_002947/6108557‑6108848

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