Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,182,059	G→T	8.6%	intergenic (‑117/+37)	PP_4560 ← / ← csbD	ribonuclease BN/stress response protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,182,059	0	G	T	8.6%	155.8 / 4.4	57	intergenic (‑117/+37)	PP_4560/csbD	ribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base G (29/23);  new base T (3/2);  total (32/25)
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.23e-01

GGGGCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTG  >  NC_002947/5181912‑5182204                                                                                                                                                     |                                                                                                                                                   ggggCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTcgc                                                                                                                                                   >  7:36567/1‑149 (MQ=255) ggggCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTcgc                                                                                                                                                   <  6:107684/149‑1 (MQ=255)  gggCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTcgcg                                                                                                                                                  >  5:221522/1‑149 (MQ=255)   ggCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCgg                                                                                                                                                 <  1:11955/149‑1 (MQ=255)   ggCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCgg                                                                                                                                                 <  8:36567/149‑1 (MQ=255)             gcAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTATTCGCGGGTTCAccc                                                                                                                                         >  3:263180/1‑147 (MQ=255)                 gccGGGTAAAAGCATAGCGTCTCCTTTCGCCGCTTCTTCGTTACACTGCGTATATCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaa                                                                                                                                    <  1:284208/147‑1 (MQ=255)                 gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaag                                                                                                                                   <  5:252171/149‑1 (MQ=255)                 gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCATCCGagaag                                                                                                                                   <  7:167652/149‑1 (MQ=255)                 gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCgag                                                                                                                                      >  7:201628/1‑146 (MQ=255)                  tCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaaga                                                                                                                                  <  2:83769/149‑1 (MQ=255)                       taaaaTCATAGCGTCTCCTTTCGCCGCTTTATCGTTACACAG‑GCACGTCTAATAGACCGGGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACTAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAATAGGCCg                                                                                                                             <  6:45277/148‑1 (MQ=255)                        aaaaTCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGcctg                                                                                                                            >  2:90615/1‑147 (MQ=255)                        aaaaTCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCgg                                                                                                                            >  8:66622/1‑149 (MQ=255)                         aaaTCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCgg                                                                                                                            <  8:201628/148‑1 (MQ=255)                                   gTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAagag                                                                                                                                 <  8:156078/133‑1 (MQ=255)                                   gTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAagag                                                                                                                                 >  7:156078/1‑133 (MQ=255)                                   atctcCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACCAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaa                                                                                                                 <  3:204202/148‑1 (MQ=255)                                    tctcCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaag                                                                                                                <  5:171077/149‑1 (MQ=255)                                       ccTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAACCGAGAAGAGGCCGGCACCGCCtcaagtc                                                                                                              >  3:134171/1‑148 (MQ=255)                                        cTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaa                                                                                                                 >  4:204202/1‑144 (MQ=255)                                              ccgcTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGcttc                                                                                                      >  2:93970/1‑149 (MQ=255)                                                     aTCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTcgc                                                                                                                                                   >  3:322466/1‑97 (MQ=255)                                                     aTCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTcgc                                                                                                                                                   <  4:322466/97‑1 (MQ=255)                                                     aTCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGACACCGCCTCAAGTCAAGGCTTCTTCACCg                                                                                               >  3:13486/1‑149 (MQ=255)                                                        gTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCg                                                                                                                             <  8:199124/116‑1 (MQ=255)                                                        gTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCg                                                                                                                             >  7:199124/1‑116 (MQ=255)                                                        gTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCAt                                                                                            >  4:145776/1‑149 (MQ=255)                                                         ttACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTTAAGAATTTGCTTAGTAATAAGTCCCGGCCTCTTCGCGTGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATc                                                                                           <  7:66622/149‑1 (MQ=255)                                                         ttACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCAt                                                                                            >  6:186199/1‑148 (MQ=255)                                                          tACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATcc                                                                                          >  7:93888/1‑149 (MQ=255)                                                             acaG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCtt                                                                                        >  2:70989/1‑148 (MQ=255)                                                               aG‑GCAATTCTAATAGACCTTGCCCGAGCTGTGCTGCTCGCGGCGGTGGTTTAAGAATTTGCTTACGAATAAGTGCCGGGCTCTTCGCGGTTTCACCCTATAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTcac                                                                                     <  3:285972/149‑1 (MQ=255)                                                                        atAATAGACCGTGCCCGAGCTTTGCTGCTCCCGGCCGTGGTTGAATAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGATAAGAGGCCGGCACCGCCTCAAGTCAAGTCTTCTTCACCGCATCCTTCACATTTccc                                                                              <  4:263180/147‑1 (MQ=255)                                                                          aaTAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCt‑t                                                                           >  6:301844/1‑148 (MQ=255)                                                                           aTAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGa                                                                         <  1:93970/149‑1 (MQ=255)                                                                            tAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGAc                                                                        >  3:234500/1‑149 (MQ=255)                                                                               aCCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGc                                                                                                                               >  3:36902/1‑92 (MQ=255)                                                                                      ccGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCGTCAACGCATGCTTGACATTTCGCA‑TGACCatctgctcg                                                              >  8:280956/1‑146 (MQ=255)                                                                                           cTGTGCTGCTCCGGGCCGTGGTTTAAGAATTTGCTTACTAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTgcgc                                                               <  6:348224/143‑1 (MQ=255)                                                                                           cTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAACAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTgcgc                                                               >  5:348224/1‑143 (MQ=255)                                                                                                 tgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGcgcctc                                                            >  3:139057/1‑140 (MQ=255)                                                                                                 tcctCCGGGCCGTTGTTGAAGAATTTGCTTACTAATATGGGCCGGCCTCTTGGCGGTTTCACCCTTTATGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGcgcctc                                                            <  4:139057/138‑1 (MQ=255)                                                                                                      ccGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCt                                              >  7:287955/1‑149 (MQ=255)                                                                                                           cGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgcc                                         >  6:140005/1‑149 (MQ=255)                                                                                                             tGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGa                                                    >  2:251403/1‑136 (MQ=255)                                                                                                             tGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGa                                                    <  1:251403/136‑1 (MQ=255)                                                                                                             tGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgcctt                                       <  5:301844/149‑1 (MQ=255)                                                                                                                   aagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAACCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgccttgcctt                                  <  4:134171/148‑1 (MQ=255)                                                                                                                    agaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCg                                >  5:110102/1‑149 (MQ=255)                                                                                                                     gaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTCTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCg                                <  1:70989/149‑1 (MQ=255)                                                                                                                        tttGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGccc                            >  5:117519/1‑149 (MQ=255)                                                                                                                                cGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCttt                    <  3:145776/149‑1 (MQ=255)                                                                                                                                 gAATAAGTGCCGGCCTCTTCGCGGGTTCGCCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCtttt                   >  2:188982/1‑149 (MQ=255)                                                                                                                                   aTAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCg                 <  6:117519/149‑1 (MQ=255)                                                                                                                                            ggCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGt         >  3:97289/1‑148 (MQ=255)                                                                                                                                              cgtctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCAc      <  5:186199/147‑1 (MQ=255)                                                                                                                                               ctctTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCgtt                                                      >  4:333647/1‑100 (MQ=37)                                                                                                                                               ctctTCTCGGGTGAACCCGCGAAGAGACCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCAcc     >  1:10668/1‑149 (MQ=255)                                                                                                                                               ctctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTGCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTGTGAGCTCATGCGCCTTGCCTTCGGCCCGCCGCTGTGCGTTGTCGGTCTcc     >  8:288179/1‑149 (MQ=255)                                                                                                                                                tctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCt    >  2:158901/1‑149 (MQ=255)                                                                                                                                                  ttCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTg  >  1:18237/1‑149 (MQ=255)                                                                                                                                                    cgcgGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCtt   <  6:110102/146‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGGGCAAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTG  >  NC_002947/5181912‑5182204

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑