Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,506,923	A→G	9.1%	intergenic (+44/‑73)	petC → / → sspA	ubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,506,923	0	A	G	9.1%	180.8 / 7.4	66	intergenic (+44/‑73)	petC/sspA	ubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base G (3/3);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

AAACTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTC  >  NC_002947/1506775‑1507068                                                                                                                                                      |                                                                                                                                                  aaaCTGGAAAGCCAGCGCATCGGCACGTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCGGGGGTGGCCTa                                                                                                                                                   <  8:82101/149‑1 (MQ=255)     tGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTGGCaa                                                                                                                                               >  1:339744/1‑148 (MQ=255)         aaGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTAtt                                                                                                                                           <  4:344418/149‑1 (MQ=255)         aaGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGGGGCCTAGCAGTAtt                                                                                                                                           <  3:1187/149‑1 (MQ=255)          aGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTg                                                                                                                                          >  4:214792/1‑149 (MQ=255)           gCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGc                                                                                                                                         <  2:280026/149‑1 (MQ=255)                 gcATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGc                                                                                                                                    <  2:339744/148‑1 (MQ=255)                  cATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGCTAGGCCa                                                                                                                                  <  8:113400/149‑1 (MQ=255)                  cATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCa                                                                                                                                  <  6:370849/149‑1 (MQ=255)                      ggCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCAcccc                                                                                                                              >  8:255752/1‑149 (MQ=255)                           cTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGgcg                                                                                                                         >  7:16252/1‑149 (MQ=255)                           cTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGgcg                                                                                                                         >  2:97840/1‑149 (MQ=255)                           cTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGgcg                                                                                                                         >  2:93020/1‑149 (MQ=255)                                 gtTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGttt                                                                                                                    >  4:200179/1‑148 (MQ=255)                                  ttGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTc                                                                                                                  <  5:224485/149‑1 (MQ=255)                                  tggCTGTAGCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGACGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTCAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTc                                                                                                                  <  5:55749/147‑1 (MQ=255)                                    gCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCa                                                                                                                 <  2:104123/148‑1 (MQ=255)                                          tCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTaaa                                                                                                           <  4:177315/148‑1 (MQ=255)                                            tGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGtt                                                                                                        >  6:9929/1‑149 (MQ=255)                                              gCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGtt                                                                                                        <  3:214792/147‑1 (MQ=255)                                              gCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTg                                                                                                       <  5:352785/148‑1 (MQ=255)                                                    ttcttcGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAg                                                                                                 >  6:368342/1‑148 (MQ=255)                                                      cttcGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGc                                                                                                <  7:255752/147‑1 (MQ=255)                                                         cGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCg                                                                                           >  6:71270/1‑149 (MQ=255)                                                           tATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTg                                                                                         >  5:143715/1‑149 (MQ=255)                                                             ttCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGc                                                                                        <  3:200179/148‑1 (MQ=255)                                                                gCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCgag                                                                                     >  6:359320/1‑148 (MQ=255)                                                                 ccTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCgagga                                                                                   <  5:359320/149‑1 (MQ=255)                                                                        gCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGc                                                                                                                                         <  1:79889/88‑1 (MQ=38)                                                                        gCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGc                                                                                                                                         >  2:79889/1‑88 (MQ=38)                                                                        gCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTc                                                                                                                  <  2:304028/111‑1 (MQ=255)                                                                         cTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACTCAGTAGTTTCT‑GCGTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTAc                                                                           <  3:32756/149‑1 (MQ=255)                                                                           cAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCTGGGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACa                                                                          >  1:279539/1‑149 (MQ=255)                                                                            aaGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATg                                                                        >  7:190703/1‑149 (MQ=255)                                                                            aaGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATg                                                                        >  7:23285/1‑149 (MQ=255)                                                                              gCGTGAATACTGGAAGGACGTGCACTGCTCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGCGCCGCTACATggg                                                                      >  5:92414/1‑149 (MQ=255)                                                                                 tGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCgag                                                                                     <  5:145737/131‑1 (MQ=255)                                                                                 tGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCgag                                                                                     >  6:145737/1‑131 (MQ=255)                                                                                    aTACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaac                                                                >  3:255979/1‑149 (MQ=255)                                                                                    aTACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaac                                                                <  1:97840/149‑1 (MQ=255)                                                                                    aTACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaac                                                                <  1:93020/149‑1 (MQ=255)                                                                                         ggAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGtt                                                                                                        >  2:39254/1‑104 (MQ=255)                                                                                         ggAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGtt                                                                                                        <  1:39254/104‑1 (MQ=255)                                                                                         ggAAGGACGTGCACTGATCACGCAGTAGGTTCT‑GGTGGTTGACCGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACa                                                           <  8:255901/149‑1 (MQ=255)                                                                                          gAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGc                                                                                                      <  1:135611/105‑1 (MQ=255)                                                                                          gAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGc                                                                                                      >  2:135611/1‑105 (MQ=255)                                                                                          gAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaaccaac                                                            <  6:143715/147‑1 (MQ=255)                                                                                          gAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACa                                                           <  8:190703/148‑1 (MQ=255)                                                                                             ggACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCAAAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAgga                                                        >  3:294485/1‑147 (MQ=255)                                                                                             ggACGGGCACTGTTCAGGCAGTTGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCTTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCTTCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAg                                                          <  4:294485/146‑1 (MQ=39)                                                                                                   gCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCt                                                  <  1:107340/148‑1 (MQ=255)                                                                                                   gCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCt                                                  >  2:107340/1‑148 (MQ=255)                                                                                                    cACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATACGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGc                                                >  3:354166/1‑149 (MQ=255)                                                                                                         aTCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTc                                           <  2:279539/149‑1 (MQ=255)                                                                                                         aTCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTc                                           >  5:149766/1‑149 (MQ=255)                                                                                                          tCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTcc                                          >  7:313541/1‑149 (MQ=255)                                                                                                          tCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTcc                                          <  8:16252/149‑1 (MQ=255)                                                                                                                               gttgAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCgagga                                                                                   >  4:138178/1‑88 (MQ=255)                                                                                                                               gttgAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCgagga                                                                                   <  3:138178/88‑1 (MQ=255)                                                                                                                                ttgAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATtct                     <  8:23285/149‑1 (MQ=255)                                                                                                                                      gcgtgcGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGTCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATtc                      <  6:376612/137‑1 (MQ=255)                                                                                                                                      gcgcCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGTCGCTACATGGGCGAAACCAACAGGTTAGCCTGCTATTACGATACCGCTGATCACTATtctagatct               >  5:376612/1‑143 (MQ=39)                                                                                                                                                tGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTgg     <  8:313541/149‑1 (MQ=255)                                                                                                                                                  gCCTGGCAATACTGCTAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGAATCCGAGCGAGGATCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCAGCTGATCACTATTCTCATCGGGTTCGCCTGGt    >  1:102642/1‑148 (MQ=38)                                                                                                                                                   ccTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTAATTTTAAGTTTCACAAGCATCCCTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTACGATCCCGCTGATCACTATTATCATCGGGTTCGCCTGGttg  >  2:172154/1‑148 (MQ=255)                                                                                                                                                    cTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGtt                                                       <  3:256783/95‑1 (MQ=255)                                                                                                                                                    cTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGtt                                                       >  4:256783/1‑95 (MQ=255)                                                                                                                                                      |                                                                                                                                                  AAACTGGAAAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTC  >  NC_002947/1506775‑1507068

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑