Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,196,067	G→C	9.7%	G43A (GGC→GCC)	xcpQ →	type II secretion pathway protein XcpQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,196,067	0	G	C	9.7%	176.5 / 8.7	62	G43A (GGC→GCC)	xcpQ	type II secretion pathway protein XcpQ
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base C (4/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

AAGCGGCAAGGATGAACGATGCCTGTTAGATACTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGG  >  NC_002947/1195922‑1196214                                                                                                                                                  |                                                                                                                                                     aaGCGGCAAGGATGAACGATGCCTGTTAGATACTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGa                                                                                                                                                   <  6:137164/149‑1 (MQ=255) aaGCGGCAAGGATGAACGATGCCTGTTAGATACTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGa                                                                                                                                                   <  8:150408/149‑1 (MQ=255)      gCAAGGATGAACGATGCCTGTTAGATACTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCa                                                                                                                                              >  4:161845/1‑149 (MQ=255)           gATGAACGATGCCTGTTAGATACTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGCCGAGTTcggtatc                                                                                                                                         <  1:75962/149‑7 (MQ=255)               aaCGATGCCTGTTAGATACTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGc                                                                                                                                      <  3:237092/148‑1 (MQ=255)                  gATGCCTGTTAGATACTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTACATGGCCTATGCCCAAGATCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGAgtgt                                                                                                                                  >  2:282945/1‑149 (MQ=255)                       cTGTTAGATACTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTcgcg                                                                                                                              >  4:302020/1‑148 (MQ=255)                               tACTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCgg                                                                                                                     <  8:52121/149‑1 (MQ=255)                               tACTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCgg                                                                                                                     <  5:177989/149‑1 (MQ=255)                                aCTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCTTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGc                                                                                                                    >  7:337911/1‑149 (MQ=255)                                 cTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGca                                                                                                                   >  6:271275/1‑149 (MQ=255)                                    cATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGcacc                                                                                                                 >  3:36958/1‑148 (MQ=255)                                    cATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGcacc                                                                                                                 <  4:36958/148‑1 (MQ=255)                                               cgcTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCAtcg                                                                                                      <  2:2223/147‑1 (MQ=255)                                                gcTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCAtcgtc                                                                                                    >  7:115572/1‑148 (MQ=255)                                                 cTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGa                                                                                                  >  5:273474/1‑149 (MQ=255)                                                 cTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGa                                                                                                  >  3:27457/1‑149 (MQ=255)                                                  tGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGAc                                                                                                 >  1:260411/1‑149 (MQ=255)                                                   gAGCCTGGCCTTGTCCATTGCCTATGCCCAAGAGCCTGTGTTCGATGACAACGGTACGCCGATGTAAGAAGTGAATTTCGTAGATACCGAACTAGGCGAGTTCATCGACAGCGTGTCGCGCATAACCGGCACC‑ACCTTCATAGTCTAcc                                                                                                >  4:153043/1‑149 (MQ=255)                                                   gAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGAcc                                                                                                >  4:147294/1‑149 (MQ=255)                                                      ccTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGAccccc                                                                                             >  6:372480/1‑149 (MQ=255)                                                       cTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCg                                                                                            >  2:1710/1‑149 (MQ=255)                                                       cTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCAGCACC‑ACCTTCATCGTCGACCCCCg                                                                                            <  2:103180/149‑1 (MQ=255)                                                               gTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCaa                                                                                      <  1:376413/147‑1 (MQ=255)                                                               gTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCACgg                                                                                    >  4:6148/1‑149 (MQ=255)                                                               gTCCATGGCCGATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCaa                                                                                      <  4:35080/147‑1 (MQ=255)                                                                      gCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGt                                                                             <  3:161845/149‑1 (MQ=255)                                                                      gCCCATGCCCTAGAGCCGGTGTTCGCTGACACCGGTACTCCGATGTACGAAGTTACTTTCGTCTATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAg                                                                               <  4:338896/147‑1 (MQ=255)                                                                       ccTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTc                                                                            >  7:253863/1‑149 (MQ=255)                                                                            gCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCgg                                                                                           >  3:231805/1‑129 (MQ=255)                                                                            gCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCgg                                                                                           <  4:231805/129‑1 (MQ=255)                                                                             cccAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGt                                                                       <  6:273474/148‑1 (MQ=255)                                                                              ccAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTgc                                                                     >  3:310051/1‑149 (MQ=255)                                                                              ccAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTgc                                                                     <  1:330854/149‑1 (MQ=255)                                                                               cAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTgcg                                                                    >  2:218913/1‑149 (MQ=255)                                                                               cAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTgcg                                                                    <  6:252671/149‑1 (MQ=255)                                                                                       ggTGTTCGATGACAACGGTTCGCCGATGTAAGAAGTGAATTTCGTCGCTTCCGAACTCGGCGATTTCATCGACATCGTGGCGCGCATCACCGGCACCTACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTc                                                             <  1:282945/149‑1 (MQ=255)                                                                                        gtgtTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGa                                                           <  2:260411/149‑1 (MQ=255)                                                                                         tgtTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGAc                                                          <  2:176528/149‑1 (MQ=255)                                                                                         tgtTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGAc                                                          <  6:153452/149‑1 (MQ=255)                                                                                               aTGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCac                                                    <  8:253863/149‑1 (MQ=255)                                                                                                  aCAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCacgac                                                 >  1:170436/1‑149 (MQ=255)                                                                                                  aCAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCacga                                                  <  3:302020/148‑1 (MQ=255)                                                                                                    aaCGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGTCAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCgg                                                                                           <  8:184583/105‑1 (MQ=255)                                                                                                     aCGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCAcgacgcc                                              <  8:115572/149‑1 (MQ=255)                                                                                                     aCGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCAcgacgcc                                              >  7:138369/1‑149 (MQ=255)                                                                                                               gATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCt                                    <  1:1710/149‑1 (MQ=255)                                                                                                                     cGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGAc                               <  5:685/148‑1 (MQ=255)                                                                                                                         gTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCtt                          >  7:309853/1‑149 (MQ=255)                                                                                                                          tGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTc                         <  5:372480/149‑1 (MQ=255)                                                                                                                            aaTTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCt                       <  1:29977/149‑1 (MQ=255)                                                                                                                             aTTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCt                       <  8:138369/148‑1 (MQ=255)                                                                                                                                 cttcgATACCTAACTCGGCGAGTTCATCGACAGAGTGTCGCGCATCACCGTCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGcc                   <  8:337911/146‑1 (MQ=255)                                                                                                                                     gATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACt              <  5:271275/149‑1 (MQ=255)                                                                                                                                        gatgAACTCGCCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACt              >  7:118217/4‑146 (MQ=255)                                                                                                                                        gatgAACTCGCCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACt              <  8:118217/143‑1 (MQ=255)                                                                                                                                        gatgAACTCGCCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTgcg           >  4:70916/4‑149 (MQ=255)                                                                                                                                        gatgAACTCGCCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTgcg           >  8:78255/4‑149 (MQ=255)                                                                                                                                        gatgAACTCGCCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCATCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTgcg           >  8:242864/4‑149 (MQ=255)                                                                                                                                         ccGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTgcgc          >  3:61381/1‑149 (MQ=255)                                                                                                                                          cGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTgcgcg         >  8:349696/1‑149 (MQ=255)                                                                                                                                          cGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTgcgcg         >  6:316265/1‑149 (MQ=255)                                                                                                                                                 ggCTAGTTCATCGACAGCGTGTCTCGCATCACCGGCAGC‑ACCTTCATCTTCGACCACCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAggg  <  3:153043/149‑1 (MQ=255)                                                                                                                                                 ggCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAggg  <  3:147294/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     AAGCGGCAAGGATGAACGATGCCTGTTAGATACTGCATGGCGGCTGCGCTGAGCCTGGCCTTGTCCATGGCCTATGCCCAAGAGCCGGTGTTCGATGACAACGGTACGCCGATGTACGAAGTGAATTTCGTCGATACCGAACTCGGCGAGTTCATCGACAGCGTGTCGCGCATCACCGGCACC‑ACCTTCATCGTCGACCCCCGGGTCAAGGGCAAGGTCACCGTGCGCACCGTCGACCTGCACGACGCCGACGCCATCTACGACATCTTCCTGGCCCAACTGCGCGCCCAGGG  >  NC_002947/1195922‑1196214

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