Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,108,695	C→G	7.8%	A315P (GCC→CCC)	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,695	0	C	G	7.8%	185.5 / 5.0	61	A315P (GCC→CCC)	PP_5359	CobW/P47K family protein
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base G (3/2);  total (32/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

CACAGCGTGCAGATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTG  >  NC_002947/6108549‑6108840                                                                                                                                                    |                                                                                                                                                     cacaGCGTGCAGATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCtt                                                                                                                                                    >  5:132414/1‑149 (MQ=255)  acacCGTGCAGATCCAACAACCCCCCTACCCGTTCGAAATCCACACCCCACCACCCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTg                                                                                                                                                   <  4:310683/149‑1 (MQ=255)         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                                                                                          >  2:166619/1‑149 (MQ=255)                                                        gcgcCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCt                                                                                          >  1:72302/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                                                                                          >  6:297761/1‑135 (MQ=255)                                                                      aTAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCt                                                                                          <  5:297761/135‑1 (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                                                     >  5:333352/1‑142 (MQ=255)                                                                                                    aaCCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCa                                                     <  6:333352/142‑1 (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                                                               >  3:183210/1‑115 (MQ=255)                                                                                                                     aTCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCggg                                                               <  4:183210/115‑1 (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          >  3:83436/4‑149 (MQ=255)                                                                                                                                         ccgcaACCCGCGTTGCAAGCCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATaa          >  4:194095/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)                                                                                                                                                 cGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTg  <  2:93626/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CACAGCGTGCAGATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACT‑GACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTG‑‑C‑GGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTG  >  NC_002947/6108549‑6108840

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