Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,182,072	A→C	9.3%	intergenic (‑130/+24)	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,072	0	A	C	9.3%	148.5 / 6.8	54	intergenic (‑130/+24)	PP_4560/csbD	ribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base A (29/19);  new base C (3/2);  total (32/22)
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.06e-01

GTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGT  >  NC_002947/5181928‑5182216                                                                                                                                                  |                                                                                                                                                  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‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGACACCGCCTCAAGTCAAGGCTTCTTCACCg                                                                                                           >  3:13486/1‑149 (MQ=255)                                        gTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCg                                                                                                                                         >  7:199124/1‑116 (MQ=255)                                        gTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCg                                                                                                                                         <  8:199124/116‑1 (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                                                                <  1:251403/136‑1 (MQ=255)                                                                                             tGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGa                                                                >  2:251403/1‑136 (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)                                                                                                                                         gTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAacac       >  5:90581/1‑149 (MQ=255)                                                                                                                                             aCCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCAGCGCATCCTTCACATTTGCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTg   <  4:234500/149‑1 (MQ=255)                                                                                                                                             aCCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTg   >  3:227277/1‑149 (MQ=255)                                                                                                                                              cccGCGAAGAGACCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGt  <  2:10668/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAG‑GCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCT‑TGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGT  >  NC_002947/5181928‑5182216

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