Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,506,927	G→A	7.4%	intergenic (+48/‑69)	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,927	0	G	A	7.4%	194.8 / 4.1	67	intergenic (+48/‑69)	petC/sspA	ubiquinol‑‑cytochrome c reductase cytochrome c1 subunit/stringent starvation protein A
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base A (2/3);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-01

AAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTC  >  NC_002947/1506783‑1507068                                                                                                                                                  |                                                                                                                                              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                                                                                                                                         >  2:79889/1‑88 (MQ=38)                                                                gCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTGGCAATACTGc                                                                                                                                         <  1:79889/88‑1 (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                                                                                     >  6:145737/1‑131 (MQ=255)                                                                         tGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCgag                                                                                     <  5:145737/131‑1 (MQ=255)                                                                            aTACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaac                                                                <  1:97840/149‑1 (MQ=255)                                                                            aTACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaac                                                                <  1:93020/149‑1 (MQ=255)                                                                            aTACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGcaac                                                                >  3:255979/1‑149 (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                                                                                                      >  2:135611/1‑105 (MQ=255)                                                                                  gAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGc                                                                                                      <  1:135611/105‑1 (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                                           >  5:149766/1‑149 (MQ=255)                                                                                                 aTCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTc                                           <  2:279539/149‑1 (MQ=255)                                                                                                  tCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTcc                                          <  8:16252/149‑1 (MQ=255)                                                                                                  tCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTcc                                          >  7:313541/1‑149 (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                                                       >  4:256783/1‑95 (MQ=255)                                                                                                                                            cTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGtt                                                       <  3:256783/95‑1 (MQ=255)                                                                                                                                                 aGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCa                            >  6:150200/1‑117 (MQ=255)                                                                                                                                                 aGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCa                            <  5:150200/117‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                              AAGCCAGCGCATCGGCACCTATGTGTTGCTGTATCTGGCTTTCTTCTTCGTATTCGCCTACTTGCTCAAGCGTGAATACTGGAAGGACGTGCACTGATCACGCAGTAGTTTCT‑GGTTGTTGAACGCGCCCGGGGTGGCCTAGCAGTATTGCCAGGCCACCCCGGGCGCGTTTTCATTTAAAGTTGCACAAGCATCCGTGCGAGGAGCCGCTACATGGGCGCAACCAACAGGTTAGCCTGCTATTCCGATCCCGCTGATCACTATTCTCATCGGGTTCGCCTGGTTC  >  NC_002947/1506783‑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: ‑