Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,143,481	C→G	5.8%	T87T (ACC→ACG)	PP_1899 →	ExbD/TolR family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,143,481	0	C	G	5.8%	213.8 / 2.3	69	T87T (ACC→ACG)	PP_1899	ExbD/TolR family protein
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base G (2/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

TTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACC  >  NC_002947/2143345‑2143628                                                                                                                                          |                                                                                                                                                      ttCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAATGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGc                                                                                                                                            >  1:106078/1‑149 (MQ=255)    gcgTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATa                                                                                                                                        <  3:118240/149‑1 (MQ=255)    gcgTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑CCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCa                                                                                                                                          >  4:18651/1‑147 (MQ=255)              gCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCTCGCTCAGGGTGGCCaggt                                                                                                                                       >  1:131157/1‑137 (MQ=255)              gCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCTCGCTCAGGGTGGCCaggt                                                                                                                                       <  2:131157/140‑4 (MQ=255)                  gAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCggg                                                                                                                          >  2:81308/1‑149 (MQ=255)                   caGCCCCCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAACGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGggc                                                                                                                         <  4:293165/148‑1 (MQ=255)                      ccagcaGCGCCGAGCCGTCGCCGGCCGACGCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCGGATGGCGTGTATTCGGTGAAGAACCACCTGCTGCCCAAAAGGGACCTGGCCACCCTGAGCGAGGGCATAGAGCGCGAGTCGggcggc                                                                                                                      <  2:106078/149‑1 (MQ=255)                        accaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGa                                                                                                                    <  4:211260/149‑1 (MQ=255)                         ccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGCGGCCATAGAGCGCGAGTCGGgcggcgt                                                                                                                    >  4:68818/1‑147 (MQ=255)                          caGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATa                                                                                                                                        <  4:300493/127‑1 (MQ=255)                          caGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATa                                                                                                                                        >  3:300493/1‑127 (MQ=255)                           aGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGacaa                                                                                                                 <  3:107515/149‑1 (MQ=255)                            gcgcCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGacaac                                                                                                                <  1:81308/149‑1 (MQ=255)                              gcCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGAcaacaa                                                                                                              <  4:232637/149‑1 (MQ=255)                               ccGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAACGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAg                                                                                                             <  3:18651/149‑1 (MQ=255)                                 gAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCTCGAGTCGGGCGGCGACAACAAGtt                                                                                                           >  8:27361/1‑149 (MQ=255)                                          ccggccgACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAATGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGc                                                                                                  >  1:262262/1‑149 (MQ=255)                                               cgACCCGGGCAA‑CCTGGTGGAAATCACCATCAGCGCCGAAGGCGAGTATTCGGTGAACAACAAACTGATGCCCAAAAGCGAACTCGCCACACTGAGCGAGCCCATAGAGCGAGAGTCGGGAGGCGACAACAAGTTGCCGCTTGCCATCa                                                                                             >  8:185051/1‑149 (MQ=255)                                                        cAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACg                                                                                  <  7:185051/149‑1 (MQ=255)                                                             cTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCaa                                                                              <  3:276171/149‑1 (MQ=255)                                                                       tCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCAt                                                                    <  3:298201/149‑1 (MQ=255)                                                                       tCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCAt                                                                    >  5:153149/1‑149 (MQ=255)                                                                         aCCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATc                                                                   <  2:186232/148‑1 (MQ=255)                                                                           cATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCAAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑CCGGCAAGACACCGCCTAAgg                                                                >  5:97396/1‑149 (MQ=255)                                                                            aTCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCAGCATCAGGc                                                               >  6:258104/1‑149 (MQ=255)                                                                            aTCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAACTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGGACCCGCATCAGGc                                                               >  8:163331/1‑149 (MQ=255)                                                                               aGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgt                                                            >  1:190808/1‑149 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCTACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCTCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgtg                                                           <  5:17112/149‑1 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGtt                                                                                                           >  3:252848/1‑103 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGtt                                                                                                           <  4:252848/103‑1 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑Ac                                                                                   >  3:194713/1‑125 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑Ac                                                                                   <  4:194713/125‑1 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgtg                                                           <  2:120923/149‑1 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgtg                                                           <  3:208673/149‑1 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgtg                                                           <  5:378321/149‑1 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgtg                                                           >  5:240158/1‑149 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgtg                                                           >  3:174238/1‑149 (MQ=255)                                                                                    cGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCt                                                              >  3:182777/1‑142 (MQ=255)                                                                                    cGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCt                                                              <  4:182777/142‑1 (MQ=255)                                                                                    cGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCa                                                       >  4:163521/1‑149 (MQ=255)                                                                                          cGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCaa                                                 >  8:38825/1‑149 (MQ=255)                                                                                             gtATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATgg                                              <  6:153149/149‑1 (MQ=255)                                                                                              tATTCGTTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGa                                             <  7:27361/149‑1 (MQ=255)                                                                                                ttCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCTCTGTCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGAt                                            <  5:258104/148‑1 (MQ=255)                                                                                                ttCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGAtg                                           >  3:167108/1‑149 (MQ=255)                                                                                                  cGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGAtgc                                          <  2:262262/148‑1 (MQ=255)                                                                                                   ggTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGAtgctg                                        >  2:36288/1‑149 (MQ=255)                                                                                                               ccTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgtg                                                           >  8:178107/1‑118 (MQ=255)                                                                                                               ccTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgtg                                                           <  7:178107/118‑1 (MQ=255)                                                                                                                   ctgcCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGttt                        <  1:36288/149‑1 (MQ=255)                                                                                                                     gcCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGttt                        >  7:220068/1‑147 (MQ=255)                                                                                                                     gcCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGttt                        >  7:214406/1‑147 (MQ=255)                                                                                                                     gcCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGttt                        >  5:265849/1‑147 (MQ=255)                                                                                                                      cccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAg                     >  2:217816/1‑149 (MQ=255)                                                                                                                       ccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCACGCATCAGGCTGAGGTCACCGCAATGGATGCAGCCGGCAAGCTCGGATTCAGc                    >  4:265861/1‑149 (MQ=255)                                                                                                                        cAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCGTCACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAg                     <  1:110634/149‑1 (MQ=255)                                                                                                                        cAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGc                    >  6:120487/1‑148 (MQ=255)                                                                                                                         aaaaGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGc                    >  5:246143/1‑147 (MQ=255)                                                                                                                             gctctatGGCCTCGCTCAGGGTGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                      >  4:265690/7‑141 (MQ=255)                                                                                                                             gctctatGGCCTCGCTCAGGGTGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                      <  3:265690/135‑1 (MQ=255)                                                                                                                                  ctcgCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCAAAAGGCTGTGGTCACCGCAATCCATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCAt         >  8:25878/4‑149 (MQ=255)                                                                                                                                     gCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGAtg                                           <  6:42812/112‑1 (MQ=255)                                                                                                                                     gCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGAtg                                           >  5:42812/1‑112 (MQ=255)                                                                                                                                     cccACCCTGACCGACGCCATAGAGCGCCAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATg        <  4:245136/146‑1 (MQ=255)                                                                                                                                       cACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                      <  2:82230/131‑1 (MQ=255)                                                                                                                                       cACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                      >  1:82230/1‑131 (MQ=255)                                                                                                                                       cACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGacca    >  7:274126/1‑149 (MQ=255)                                                                                                                                        aCCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGaccac   >  5:69181/1‑149 (MQ=255)                                                                                                                                         cccTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGaccacc  >  8:20275/1‑149 (MQ=255)                                                                                                                                         cccTGACCGAGGACATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGaccacc  <  7:163331/149‑1 (MQ=255)                                                                                                                                          ccTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCgg                           <  6:171658/123‑1 (MQ=255)                                                                                                                                          ccTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCgg                           >  5:171658/1‑123 (MQ=255)                                                                                                                                          |                                                                                                                                                      TTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACC  >  NC_002947/2143345‑2143628

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