Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,732,263	A→G	6.2%	I98T (ATC→ACC)	hutI ←	imidazolonepropionase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,732,263	0	A	G	6.2%	191.9 / 2.5	64	I98T (ATC→ACC)	hutI	imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base G (2/2);  total (31/33)
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.81e-01

TGCGCTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCG  >  NC_002947/5732117‑5732411                                                                                                                                                   |                                                                                                                                                       tGCGCTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATa                                                                                                                                                      >  3:63112/1‑149 (MQ=255)        gTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑c                                                                                                                                             >  1:198208/1‑149 (MQ=255)             ccAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAG‑GGTGGCGGTATCG‑CCagcactgt                                                                                                                                        <  8:280362/149‑8 (MQ=255)              cAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGAt                                                                                                                                                       >  4:54148/1‑135 (MQ=255)              cAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGAt                                                                                                                                                       <  3:54148/135‑1 (MQ=255)               aGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTgc                                                                                                                                      <  5:252766/149‑1 (MQ=255)                   ccAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTgcggcg                                                                                                                                  <  7:237509/149‑1 (MQ=255)                    cAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGa                                                                                                                                 >  5:100296/1‑149 (MQ=255)                         cGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCg                                                                                                                            >  2:182073/1‑149 (MQ=255)                              ccGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCa                                                                                                                         >  2:62289/1‑147 (MQ=255)                                      gATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCAc                                                                                                                 <  1:182073/147‑1 (MQ=255)                                        tctcGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGccc                                                                                                             >  3:173312/1‑149 (MQ=255)                                        tctcGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGccc                                                                                                             >  5:229966/1‑149 (MQ=255)                                           cGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTg                                                                                                           >  5:207869/1‑148 (MQ=255)                                            gAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCt                                                                                                            <  6:100296/146‑1 (MQ=255)                                              gggTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAgg                                                                                                       >  2:204178/1‑149 (MQ=255)                                              gggTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAgg                                                                                                       <  1:32791/149‑1 (MQ=255)                                               ggtggtCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAgg                                                                                                       <  1:62289/149‑1 (MQ=255)                                                   gtCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGtt                                                                                                   <  4:127396/148‑1 (MQ=255)                                                    tCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAggagatc                                                                                                  >  7:110934/1‑143 (MQ=255)                                                    tCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAg                                                                                                        <  8:110934/142‑1 (MQ=255)                                                      aCGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCt                                                                                                <  6:207869/148‑1 (MQ=255)                                                        gCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGTAGGCGTTGCTCg                                                                                              >  6:92390/1‑148 (MQ=255)                                                                      gTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCg                                                                                <  2:52417/148‑1 (MQ=255)                                                                       tGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGt                                                                              >  6:122476/1‑149 (MQ=255)                                                                         ccTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTg                                                                            >  6:205176/1‑149 (MQ=255)                                                                             aaCGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTgccgc                                                                        <  2:198208/149‑1 (MQ=255)                                                                             aaCGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTgcc                                                                          <  7:241256/147‑1 (MQ=255)                                                                               cgcgCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGccgccg                                                                      <  1:204178/149‑1 (MQ=255)                                                                                gcgcTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCgg                                                                                                                           >  6:312668/1‑95 (MQ=255)                                                                                gcgcTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCgg                                                                                                                           <  5:312668/95‑1 (MQ=255)                                                                                   cgtgTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAaca                                                                  <  4:38614/146‑1 (MQ=255)                                                                                          gCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGGGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCTCTGCGGTTGCCGCCGAACACCGCATg                                                           <  7:201817/149‑1 (MQ=255)                                                                                           cACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATgg                                                          <  4:37352/149‑1 (MQ=255)                                                                                                 cGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCa                                                             >  4:134982/1‑140 (MQ=255)                                                                                                 cGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCa                                                             <  3:134982/140‑1 (MQ=255)                                                                                                  gAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCGGCAGGCGTGGCTCGAACTCGCCGCTGCGGGTGCCGCCGACCACCGCATGGGGGTgg                                                    >  4:185810/1‑148 (MQ=255)                                                                                                   aacaacTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCa                                                  <  5:205176/149‑1 (MQ=255)                                                                                                        cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATg                                                           >  5:175856/1‑135 (MQ=255)                                                                                                        cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATg                                                           >  5:171368/1‑135 (MQ=255)                                                                                                        cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATg                                                           <  6:171368/135‑1 (MQ=255)                                                                                                        cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATg                                                           <  6:175856/135‑1 (MQ=255)                                                                                                             cTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGc                                        >  5:24887/1‑149 (MQ=255)                                                                                                              ttCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGcc                                       >  4:43500/1‑149 (MQ=255)                                                                                                                 gCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATgg                                                          >  1:57092/1‑127 (MQ=255)                                                                                                                 gCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATgg                                                          <  2:57092/127‑1 (MQ=255)                                                                                                                     gCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGa                                  <  6:229966/147‑1 (MQ=255)                                                                                                                     gCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGa                                  <  8:206247/147‑1 (MQ=255)                                                                                                                        gCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCgg                                    <  2:141139/142‑1 (MQ=255)                                                                                                                        gCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCgg                                    >  1:141139/1‑142 (MQ=255)                                                                                                                        gCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGAc                             >  6:127196/1‑149 (MQ=255)                                                                                                                         cACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGAc                             <  3:21947/148‑1 (MQ=255)                                                                                                                         cACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGAc                             >  4:21947/1‑148 (MQ=255)                                                                                                                          aCGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGAccc                           <  6:24887/149‑1 (MQ=255)                                                                                                                              gTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGTAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGc                       <  5:92390/149‑1 (MQ=255)                                                                                                                                  cccGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGCTCAGCCCGGGAGTGACCCAGGcacca                   >  6:281719/1‑149 (MQ=255)                                                                                                                                   ccGCACAG‑GGTGGCGGTATCG‑CCAGCACTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGcaccac                  >  1:307420/1‑149 (MQ=38)                                                                                                                                         aGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTTGCAGTCGATCAGCCCGGGATTGACCCAGGCACCACCCAGGt            >  7:219618/1‑149 (MQ=255)                                                                                                                                         aGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGt            >  3:178588/1‑149 (MQ=255)                                                                                                                                          agggtGGCGGTATCG‑CCAGCACTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGt                                                >  6:58359/5‑113 (MQ=38)                                                                                                                                          agggtGGCGGTATCG‑CCAGCACTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGt                                                <  5:58359/109‑1 (MQ=38)                                                                                                                                                gCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGt     >  4:220711/1‑149 (MQ=255)                                                                                                                                                 cGATACCG‑CCA‑CCCTGTGCGGCGATTTCTGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTg    <  5:127196/149‑1 (MQ=255)                                                                                                                                                  gATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGc   <  4:173312/149‑1 (MQ=255)                                                                                                                                                   aTACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCg  <  3:43500/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                       TGCGCTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCG  >  NC_002947/5732117‑5732411

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