Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,143,487	C→G	5.2%	S89R (AGC→AGG)	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,487	0	C	G	5.2%	250.9 / 2.9	77	S89R (AGC→AGG)	PP_1899	ExbD/TolR family protein
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base G (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

TTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGA  >  NC_002947/2143345‑2143633                                                                                                                                                |                                                                                                                                                     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                                                                                                                                            <  2:131157/140‑4 (MQ=255)              gCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCTCGCTCAGGGTGGCCaggt                                                                                                                                            >  1:131157/1‑137 (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                                                                >  5:240158/1‑149 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgtg                                                                <  5:378321/149‑1 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgtg                                                                <  3:208673/149‑1 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgtg                                                                <  2:120923/149‑1 (MQ=255)                                                                                gcgcCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCtgtg                                                                >  3:174238/1‑149 (MQ=255)                                                                                    cGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCt                                                                   <  4:182777/142‑1 (MQ=255)                                                                                    cGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCt                                                                   >  3:182777/1‑142 (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                             >  5:265849/1‑147 (MQ=255)                                                                                                                     gcCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGttt                             >  7:214406/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                           >  1:82230/1‑131 (MQ=255)                                                                                                                                       cACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                           <  2:82230/131‑1 (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                                >  5:171658/1‑123 (MQ=255)                                                                                                                                          ccTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCgg                                <  6:171658/123‑1 (MQ=255)                                                                                                                                             gAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                           <  2:367231/125‑1 (MQ=255)                                                                                                                                             gAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                           >  1:367231/1‑125 (MQ=255)                                                                                                                                              aGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGa  <  7:38825/149‑1 (MQ=255)                                                                                                                                              aGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGa  <  1:217816/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     TTCGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAA‑GCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCGACAACAAGTTGCCGCTGGCCATCAGCGCCG‑‑ACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGA  >  NC_002947/2143345‑2143633

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