Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,477,813	T→G	5.9%	Q60P (CAG→CCG) ‡	PP_3968 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,477,813	0	T	G	5.9%	246.3 / 4.3	81	Q60P (CAG→CCG) ‡	PP_3968	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base G (2/3);  total (40/41)
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.56e-01

GCGATAGAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCA  >  NC_002947/4477667‑4477960                                                                                                                                                   |                                                                                                                                                    gCGATAGAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGa                                                                                                                                                   >  4:44269/1‑149 (MQ=255)  cGATAGAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGAt                                                                                                                                                  <  5:305436/149‑1 (MQ=255)    aTAGAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATgc                                                                                                                                                <  7:110085/149‑1 (MQ=255)       gAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCa                                                                                                                                             <  4:328401/149‑1 (MQ=255)         agtaccGAGCGCTGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGt                                                                                                                                           >  3:343753/6‑149 (MQ=255)           acTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCt                                                                                                                                         <  4:71948/149‑1 (MQ=255)              cGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGATGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAg                                                                                                                                      >  8:358154/1‑149 (MQ=255)              cGAGCGCAGCGCTCGGTACTCACCGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATAAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCGCATCAGCAGTCTGAg                                                                                                                                      <  3:39868/149‑1 (MQ=255)               gAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCGCATCAGCAGTCTGa                                                                                                                                       <  2:242306/147‑1 (MQ=255)               gAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCGCATCAGCAGTCTGAg                                                                                                                                      <  8:160769/148‑1 (MQ=255)                   gctgcgcTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTa                                                                                                                                 <  6:336075/146‑1 (MQ=255)                          tCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCTCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCatat                                                                                                                          <  8:275984/149‑1 (MQ=255)                            ggTACTCATAGACCAATTGATCCCTCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCatat                                                                                                                          <  8:60286/147‑1 (MQ=255)                                 tCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCgg                                                                                                                    >  7:338034/1‑148 (MQ=255)                                 tCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCgg                                                                                                                    >  2:234177/1‑148 (MQ=255)                                   agagACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCaa                                                                                                                 >  2:39217/1‑149 (MQ=255)                                       aCCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGt                                                                                                             <  4:293107/149‑1 (MQ=255)                                         cAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTTCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTgc                                                                                                           <  2:13687/149‑1 (MQ=255)                                         cAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTgc                                                                                                           >  7:153777/1‑149 (MQ=255)                                           aTTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTgcg                                                                                                          >  6:302714/1‑148 (MQ=255)                                                 ccATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAAtc                                                                                                   >  3:203638/1‑149 (MQ=255)                                                   aTCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAAtctt                                                                                                 >  3:87468/1‑149 (MQ=255)                                                   aTCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTCATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGAAGAAgctt                                                                                                 >  6:229104/1‑149 (MQ=255)                                                        aaaGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAAtcttct                                                                                               >  4:287009/1‑146 (MQ=255)                                                         aaGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCaga                                                                                           <  2:145859/149‑1 (MQ=255)                                                          aGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCagag                                                                                          <  2:77840/149‑1 (MQ=255)                                                            cccGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGt                                                                                         >  6:30227/1‑148 (MQ=255)                                                                cTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTc                                                                                    <  4:87468/149‑1 (MQ=255)                                                                     aaaCGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCTGAGTGATTCCTGa                                                                                <  6:124989/148‑1 (MQ=255)                                                                     aaaCGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGa                                                                                <  1:234177/148‑1 (MQ=255)                                                                     aaaCGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAg                                                                               <  5:154001/149‑1 (MQ=255)                                                                     aaaCGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAg                                                                               >  3:92830/1‑149 (MQ=255)                                                                     aaaCGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAg                                                                               <  5:120737/149‑1 (MQ=255)                                                                     aaaCGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAg                                                                               >  6:154001/1‑149 (MQ=255)                                                                       aCGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTATCAGAGTGATTCCTGAGGc                                                                             >  8:79278/1‑149 (MQ=255)                                                                          aGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCgg                                                                          <  4:343753/149‑1 (MQ=255)                                                                               tGGCATGGGGCTGAGCCGCCGTATCGCGTCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGTGTGATTCCTGAGGCCGGCTGaa                                                                     <  3:259353/149‑1 (MQ=255)                                                                                 gCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATc                                                                   <  3:287009/149‑1 (MQ=255)                                                                                       gTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGt                                                             <  1:64243/149‑1 (MQ=255)                                                                                       gTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATATTCACAGAGAGATTCCTGAGGCCGGCTGAAACCTGAGt                                                             >  8:56824/1‑149 (MQ=255)                                                                                        tCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTg                                                            >  7:100829/1‑149 (MQ=255)                                                                                        tCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTg                                                            <  8:338034/149‑1 (MQ=255)                                                                                         cTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGc                                                                         <  4:33431/135‑1 (MQ=255)                                                                                         cTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGc                                                                         >  3:33431/1‑135 (MQ=255)                                                                                             gccgccGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGt                                                         >  7:142127/1‑147 (MQ=255)                                                                                              ccgccgTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAg                                                      >  1:161535/1‑149 (MQ=255)                                                                                                gccgTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTc                                                    >  6:79175/1‑149 (MQ=255)                                                                                                    tATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTc                                                <  8:153777/149‑1 (MQ=255)                                                                                                        gctgctCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACTGTGCGCAGAATCTTCTCAGAGTGATTc                                                                                    <  5:178221/109‑1 (MQ=255)                                                                                                        gctgctCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACAGTGCGCAGAATCTTCTCAGAGTGATTc                                                                                    >  6:178221/1‑109 (MQ=255)                                                                                                          tgctCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGAcc                                          >  5:145806/1‑149 (MQ=255)                                                                                                            ctCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGt                                                                                                             >  5:27700/1‑80 (MQ=255)                                                                                                            ctCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGt                                                                                                             <  6:27700/80‑1 (MQ=255)                                                                                                             tCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGaa                                       >  1:56721/1‑149 (MQ=255)                                                                                                              cATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAAc                                      >  7:73252/1‑149 (MQ=255)                                                                                                               aTTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACg                                     >  7:121617/1‑149 (MQ=255)                                                                                                                   gagGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGc                                  <  5:79175/148‑1 (MQ=255)                                                                                                                         ctctGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGc                                  <  6:139703/142‑1 (MQ=255)                                                                                                                         ctctGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGc                                  >  5:139703/1‑142 (MQ=255)                                                                                                                         ctctGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAAtcc                           >  3:318778/1‑149 (MQ=255)                                                                                                                         ctctGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAAtcc                           >  4:231826/1‑149 (MQ=255)                                                                                                                         ctctGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAAtcc                           >  1:337111/1‑149 (MQ=255)                                                                                                                          tctGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAAtcct                          >  5:79373/1‑149 (MQ=255)                                                                                                                            tGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAAtcct                          <  5:302714/147‑1 (MQ=255)                                                                                                                                  gcgcCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGc                   <  1:39217/148‑1 (MQ=255)                                                                                                                                   gtactgctcAGACTGCGCATCAGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCt                 >  1:171479/7‑149 (MQ=255)                                                                                                                                   cgcCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCt                 <  7:56824/149‑1 (MQ=255)                                                                                                                                   cagcacatcAGACTGCTCATCAGCAGTCTGAGCAGTACGCATATCAGCGTCAACGTTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGTGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCt                 <  3:382320/142‑1 (MQ=255)                                                                                                                                    gcCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTg                <  5:30227/149‑1 (MQ=255)                                                                                                                                     ccTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTcc                                                   <  8:189301/113‑1 (MQ=255)                                                                                                                                     ccTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTcc                                                   >  7:189301/1‑113 (MQ=255)                                                                                                                                     ccTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTg                >  4:370247/1‑148 (MQ=255)                                                                                                                                      cTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGga              <  2:56721/149‑1 (MQ=255)                                                                                                                                      cTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGga              <  3:370247/149‑1 (MQ=255)                                                                                                                                           cAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATg         >  5:377060/1‑149 (MQ=255)                                                                                                                                            aGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAAt          <  3:231826/147‑1 (MQ=255)                                                                                                                                            aGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGc        >  7:177546/1‑149 (MQ=255)                                                                                                                                            aGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGc        >  3:291952/1‑149 (MQ=255)                                                                                                                                            aGACTGCGCATCAGCAGTCTGAGAAG‑GCGCATATCAGCCGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAACTCCGTCTTGACCGAACGTGCGAAATACTCCCAAGCCTGGCGAATGCt       >  2:141816/1‑149 (MQ=255)                                                                                                                                              aCTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCtt      <  8:142127/149‑1 (MQ=255)                                                                                                                                              aCTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCtt      <  8:73252/149‑1 (MQ=255)                                                                                                                                               ctgctgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTc     <  6:79373/149‑1 (MQ=255)                                                                                                                                                 gctgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTccc   <  4:53610/149‑1 (MQ=255)                                                                                                                                                  ctgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTccc   <  7:358154/148‑1 (MQ=255)                                                                                                                                                  ctgATGCGCAGTATGAGGAGTACGCATATCAGCGTCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTcc    <  5:229104/147‑1 (MQ=255)                                                                                                                                                   tgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGc        <  3:193966/142‑1 (MQ=255)                                                                                                                                                   tgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGc        >  4:193966/1‑142 (MQ=255)                                                                                                                                                   tgATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCa  >  6:311598/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GCGATAGAACACTCGAGCGCAGCGCTCGGTACTCAGAGACCAATTGATCCATCGAAAAGCCCGCTATCAAACGAGTAATGGCATGGGTCTGAGCCGCCGTATCGCTGCTCATTAGAGGGCCTCTGCCTTGCGCCTTCTCAGACTGCTGATGCGCAGTCTGAGCAGTACGCATATCAGCGGCAACGGTGCGCAGAATCTTCTCAGAGTGATTCCTGAGGCCGGCTGAATCCTGAGTGGGTAAGTCCGTCTTGACCGAACGTGCGAAATCCTCCCAAGCCTGGAGAATGCTTCCCA  >  NC_002947/4477667‑4477960

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