Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,486,821	G→A	9.8%	intergenic (+54/‑130)	spuC‑I → / → PP_2181	polyamine:pyruvate transaminase/Cro/CI family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,486,821	0	G	A	9.8%	163.8 / 7.3	62	intergenic (+54/‑130)	spuC‑I/PP_2181	polyamine:pyruvate transaminase/Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base A (2/4);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

CCACCGCTGGTGATCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAA  >  NC_002947/2486678‑2486956                                                                                                                                                  |                                                                                                                                         ccACCGCTGGTGATCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCgg                                                                                                                                     >  8:222380/1‑149 (MQ=255)  cACCGCTGGTGATCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGt                                                                                                                                    <  4:215317/149‑1 (MQ=255)     cGCTGGTGATCAGTATCGAAGAGATCGACGAACTTGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCGCAGGCGGTGC‑t                                                                                                                                <  5:325392/149‑1 (MQ=255)     cGCTGGTGATCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGAAGCCTAATAGGCTGTTAGGTAAAAGCAGAGGGCCTGGTTAGCTGTTCAAAGCACCTGTGAAATAAAACACCGCCTGCGCGGTg                                                                                                                                   >  1:117118/1‑149 (MQ=37)     cGCTGGTGATCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑t                                                                                                                                <  4:160371/149‑1 (MQ=255)      gCTGGTGATCAGTATCGAAGAGATGGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCGAGCACCGCGCAGGCGGTGC‑tt                                                                                                                               <  4:157523/149‑2 (MQ=39)      gCTGGTGATCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑Tc                                                                                                                               <  5:214718/149‑1 (MQ=255)           tGATCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCg                                                                                                                          <  6:134469/149‑1 (MQ=255)              tCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTg                                                                                                                       >  2:205185/1‑149 (MQ=255)                       aaGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCg                                                                                                              >  2:146031/1‑149 (MQ=255)                            aTCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACAt                                                                                                         <  1:370280/149‑1 (MQ=255)                               gacgaACTGGTGGATAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACATCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACAttt                                                                                                       <  3:366043/148‑1 (MQ=255)                                acgaACTGGTGGAGAAATCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCa                                                                                                     <  6:309882/149‑1 (MQ=255)                                acgaACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGTCGCCATCGTACATTTCa                                                                                                     >  8:21540/1‑149 (MQ=255)                                acgaACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCa                                                                                                     >  8:255505/1‑149 (MQ=255)                                 cgaACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCAc                                                                                                    <  5:143408/149‑1 (MQ=255)                                     cTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCa                                                                                                <  3:77751/149‑1 (MQ=255)                                        gtggAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGCGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAgg                                                                                              >  7:140267/1‑148 (MQ=255)                                        gtggAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAgg                                                                                              >  8:226732/1‑148 (MQ=255)                                        gtggAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGc                                                                                             >  2:104796/1‑149 (MQ=255)                                         tggAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCa                                                                                            >  7:13047/1‑149 (MQ=255)                                           gagaAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCAt                                                                                                                >  3:349803/1‑127 (MQ=255)                                           gagaAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCAt                                                                                                                <  4:349803/127‑1 (MQ=255)                                           gagaAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCAcc                                                                                          >  3:360272/1‑149 (MQ=255)                                            agaAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCGAGCACCGCGCAGGCGGTGC‑TTGATctcaca                                                                                                                      <  2:353823/120‑6 (MQ=37)                                            agaAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCGAGCACCGCGCAGGCGGTGC‑TTGATctcaca                                                                                                                      >  1:353823/1‑115 (MQ=37)                                                         gTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCAtt                                                                            >  8:67359/1‑149 (MQ=255)                                                         gTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCAtt                                                                            >  3:284835/1‑149 (MQ=255)                                                            ccTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTg                                                                         >  8:255716/1‑149 (MQ=255)                                                                gacctga‑‑cctATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCt                                                                       <  6:55850/147‑1 (MQ=255)                                                                 acctga‑‑cgtATGAGGCTGTTCGGTAAGAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCTAGCACGGCCTGCGCGGTGC‑TCGATCGGTGCCGGCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTtg                                                                    <  7:43136/149‑1 (MQ=255)                                                                    tga‑‑cctATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTtgtg                                                                  >  7:41534/1‑148 (MQ=255)                                                                     ga‑‑cctATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCGCAGGCGGTGC‑ttg                                                                                                                              <  1:7019/87‑3 (MQ=38)                                                                     ga‑‑cctATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCGCAGGCGGTGC‑ttg                                                                                                                              >  2:7019/1‑85 (MQ=38)                                                                      a‑‑cctATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCAAAATGCTTtgcgg                                                                 <  4:360272/147‑4 (MQ=255)                                                                        ccctATGAGGCTGTTCTGTAACAGCAGATGGCCTGGTTATCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCTGGCACCTTCCGCAAAGCATTTTGCTTTGTGa                                                                 <  4:21458/146‑1 (MQ=255)                                                                           tATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCttt                                                                     <  8:129655/140‑1 (MQ=255)                                                                           tATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCttt                                                                     >  7:129655/1‑140 (MQ=255)                                                                            aTGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGaa                                                           >  6:350994/1‑149 (MQ=255)                                                                             tGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTtgt                                                                   >  3:146094/1‑140 (MQ=255)                                                                             tGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTtgt                                                                   <  4:146094/140‑1 (MQ=255)                                                                               cggCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCGGTTCGCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATc                                                         <  8:177672/147‑1 (MQ=255)                                                                                ggCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCt                                                       <  7:255505/149‑1 (MQ=255)                                                                                 gCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTg                                                      >  2:230120/1‑149 (MQ=255)                                                                                          tAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCt                                             <  8:13047/149‑1 (MQ=255)                                                                                                 agagGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATc                                      >  7:179931/1‑149 (MQ=255)                                                                                                 agagGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATc                                      >  5:237290/1‑149 (MQ=255)                                                                                                       ccTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTtgt                                                                   >  3:74735/1‑114 (MQ=255)                                                                                                       ccTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTtgt                                                                   <  4:74735/114‑1 (MQ=255)                                                                                                        cTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGt                                                     <  6:258938/127‑1 (MQ=255)                                                                                                        cTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGt                                                     >  5:258938/1‑127 (MQ=255)                                                                                                           gTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGTCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGt                            <  7:21540/149‑1 (MQ=255)                                                                                                            ttAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTc                           <  3:157658/149‑1 (MQ=255)                                                                                                            ttAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTc                           <  1:344059/149‑1 (MQ=255)                                                                                                              aGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCt                          <  7:255716/148‑1 (MQ=255)                                                                                                                cTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGa                        >  8:265489/1‑148 (MQ=255)                                                                                                                cTGTTCTCAGCACCTGTGAGATAAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTATGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGa                        >  3:164822/1‑148 (MQ=255)                                                                                                                        aGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGCTTCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTc                   >  1:260086/1‑146 (MQ=255)                                                                                                                        aGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTcgcg                <  7:222380/148‑1 (MQ=255)                                                                                                                                tgAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCa       <  5:350994/149‑1 (MQ=255)                                                                                                                                     tCAAGCACCGCCTTCGCGGTTC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGaaa  <  7:28036/149‑1 (MQ=255)                                                                                                                                     tCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGaaa  <  1:205185/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                         CCACCGCTGGTGATCAGTATCGAAGAGATCGACGAACTGGTGGAGAAAGCCCGCAAGTGCCTGGACCTGA‑‑CCTATGAGGCTGTTCGGTAACAGCAGAGGGCCTGGTTAGCTGTTCTCAGCACCTGTGAGATCAAGCACCGCCTGCGCGGTGC‑TCGATCGGTGCCGCCATCGTACATTTCACGTCAGGCACCTTCCGCAAAGCATTTTGCTTTGTGAAATGAATCCTGTACTTTTCTGGGCATCCATTTCCCGTCTGAAAGTCGCGATTCAACCATGAAA  >  NC_002947/2486678‑2486956

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