Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,295,088	G→C	6.0%	Q190H (CAG→CAC)	PP_4668 →	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,295,088	0	G	C	6.0%	209.4 / 3.2	66	Q190H (CAG→CAC)	PP_4668	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (37/25);  new base C (2/2);  total (39/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 = 3.88e-01

CATCAGCCTGTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACGGCAGCGCTGC  >  NC_002947/5294941‑5295236                                                                                                                                                    |                                                                                                                                                      catcaGCCTGTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAg                                                                                                                                                       <  7:87956/148‑1 (MQ=255)     aGCCTGTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGAc                                                                                                                                                  <  1:39428/149‑1 (MQ=255)     aGCCTGTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGAc                                                                                                                                                  >  2:39428/1‑149 (MQ=255)       ccTGTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACg                                                                                                                                                 >  1:160613/1‑148 (MQ=255)       ccTGTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACg                                                                                                                                                 >  1:60160/1‑148 (MQ=255)          gTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATc                                                                                                                                             >  1:133489/1‑149 (MQ=255)           tCCAAGCGTGAGGCGGACATCGGCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGAt                                                                                                                                              <  2:241634/147‑1 (MQ=255)             cAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGTCACCTGGCGGATcggc                                                                                                                                          <  2:301656/149‑4 (MQ=255)              aaGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGc                                                                                                                                         >  2:54763/1‑149 (MQ=255)              aaGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGc                                                                                                                                         >  7:28594/1‑149 (MQ=255)                            cATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTtca                                                                                                                           >  1:541/1‑149 (MQ=255)                             aTCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTtcat                                                                                                                          <  2:160613/149‑1 (MQ=255)                                gTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAg                                                                                                                       >  2:245232/1‑149 (MQ=255)                                  cATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGTCACCTGGCGGATcgg                                                                                                                                           >  2:234609/1‑125 (MQ=255)                                  cATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGTCACCTGGCGGATcgg                                                                                                                                           <  1:234609/127‑3 (MQ=255)                                   aTCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGtt                                                                                                                     <  8:235753/148‑1 (MQ=255)                                    tCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTAc                                                                                                                   >  3:224822/1‑149 (MQ=255)                                                 gCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTgg                                                                                                     >  5:91257/1‑149 (MQ=255)                                                   ccGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCg                                                                                                   >  1:3410/1‑149 (MQ=255)                                                    cGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGt                                                                                                  <  4:261347/149‑1 (MQ=255)                                                     gAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGtt                                                                                                 <  1:245232/149‑1 (MQ=255)                                                     gAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGtt                                                                                                 <  6:91257/149‑1 (MQ=255)                                                       gcctGGGCCCTATCGGTGCGCCACGCTGTGCCACTCCCGCCTGCGGCTGTATGCCACCCCGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGtt                                                                                                 <  4:224822/144‑1 (MQ=255)                                                            ggCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTc                                                                                           <  2:60160/148‑1 (MQ=255)                                                             gCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTcc                                                                                          >  6:284427/1‑148 (MQ=255)                                                             gCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTcc                                                                                          >  2:283795/1‑148 (MQ=255)                                                                  aTGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAgc                                                                                      <  1:283795/147‑1 (MQ=255)                                                                     tgtgCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCg                                                                                         >  7:33465/1‑141 (MQ=255)                                                                     tgtgCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCg                                                                                         <  8:33465/141‑1 (MQ=255)                                                                     tgtgCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGctgctg                                                                                 >  7:40508/1‑149 (MQ=255)                                                                        gctgcAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTAc                                                                              >  2:294897/1‑149 (MQ=255)                                                                         ctgcAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGAc                                                                                                         <  1:88369/121‑1 (MQ=255)                                                                         ctgcAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGAc                                                                                                         >  2:88369/1‑121 (MQ=255)                                                                         ctgcAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTcc                                                                                          <  5:150684/136‑1 (MQ=255)                                                                         ctgcAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTcc                                                                                          >  6:150684/1‑136 (MQ=255)                                                                           gcAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTg                                                                           >  6:52828/1‑149 (MQ=255)                                                                            cAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACGCTGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCt                                                                            <  1:54763/148‑1 (MQ=255)                                                                            cAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTg                                                                           <  8:40508/148‑1 (MQ=255)                                                                                 tgtgCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAAc                                                                     >  1:72584/1‑149 (MQ=255)                                                                                   tgCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAAc                                                                     <  5:284427/147‑1 (MQ=255)                                                                                     cGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGa                                                                 >  1:299648/1‑149 (MQ=255)                                                                                         tACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATccc                                                             >  1:127844/1‑149 (MQ=255)                                                                                            cGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCaa                                                          <  1:294897/149‑1 (MQ=255)                                                                                              ccTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACg                                                        <  2:133489/149‑1 (MQ=255)                                                                                                          tGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACAtt                                            >  2:102934/1‑149 (MQ=255)                                                                                                           gCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACAttt                                           >  4:17382/1‑149 (MQ=255)                                                                                                             cACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTgc                                         >  1:106838/1‑149 (MQ=255)                                                                                                              aCCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTgcg                                        >  6:121851/1‑149 (MQ=255)                                                                                                               cccAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTgcg                                        >  1:176691/1‑148 (MQ=255)                                                                                                                ccAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTgcg                                        >  1:296245/1‑147 (MQ=255)                                                                                                                         ccTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCg                             <  2:296245/149‑1 (MQ=255)                                                                                                                          cTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGt                            <  2:114022/149‑1 (MQ=255)                                                                                                                          cTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGt                            >  1:114022/1‑149 (MQ=255)                                                                                                                           tGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGt                            >  5:62942/1‑148 (MQ=255)                                                                                                                             gACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGa                          <  1:102934/148‑1 (MQ=255)                                                                                                                                  ccACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACAttt                                           >  4:235158/1‑126 (MQ=255)                                                                                                                                  ccACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACAttt                                           <  3:235158/126‑1 (MQ=255)                                                                                                                                   cACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCa                   >  3:201018/1‑149 (MQ=255)                                                                                                                                    aCGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATcgcg                     >  6:230355/1‑146 (MQ=255)                                                                                                                                    aCGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAg                  >  1:41855/1‑149 (MQ=255)                                                                                                                                    aCGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAg                  <  5:121851/149‑1 (MQ=255)                                                                                                                                    aCGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAg                  >  3:288581/1‑149 (MQ=255)                                                                                                                                        gcCGATCGGGGAGGTGAAGGATCTGGAGAATCAGCCCTTGATCATTTTCG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTTGCTGGCCAACGTTATCGCGAACGCTTGCGCTCTTTTTCGCAGTAACAGCGTGATc                        <  4:283635/139‑1 (MQ=255)                                                                                                                                        gcCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGATGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATc                        >  3:283635/1‑139 (MQ=255)                                                                                                                                         ccGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTacac             <  6:62942/149‑1 (MQ=255)                                                                                                                                          agATCCGTCACCTGGCGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACg            >  1:113802/2‑149 (MQ=255)                                                                                                                                             tCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGTTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACgg           >  5:261236/1‑147 (MQ=255)                                                                                                                                                    ggTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACGGCAGCGCTGc  <  4:288581/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CATCAGCCTGTCCAAGCGTGAGGCGGACATCGTCATCGCGCTGGAGCGGCCCGAGCATGGGCCTTATGTGTGCTGCAAGCTGTGCGACTACCGCCTGCGGCTGTATGCCACCCAGGGCTACCTGGACAGCCACGCGCCGATCCGCCAGGTGACGGATCTGGCCAAACACCCCTTCATCAGTTACG‑TGGATGACCTGGCGTTCAGCTCCGAGCTGCTGTACCTGGCCAACCTGATCCCCAACGCTAGCGCACATTTGCGCAGTACCAGCGTGATCGCGCAGTACACGGCAGCGCTGC  >  NC_002947/5294941‑5295236

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