Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,732,266	C→T	6.4%	G97D (GGT→GAT)	hutI ←	imidazolonepropionase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,732,266	0	C	T	6.4%	184.7 / 2.1	62	G97D (GGT→GAT)	hutI	imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base T (2/2);  total (29/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 = 4.48e-01

GTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGT  >  NC_002947/5732124‑5732413                                                                                                                                               |                                                                                                                                                      gTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑c                                                                                                                                               >  1:198208/1‑149 (MQ=255)      ccAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAG‑GGTGGCGGTATCG‑CCagcactgt                                                                                                                                          <  8:280362/149‑8 (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                                                                                                               >  5:229966/1‑149 (MQ=255)                                 tctcGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGccc                                                                                                               >  3:173312/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                                                                                                                             <  5:312668/95‑1 (MQ=255)                                                                         gcgcTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCgg                                                                                                                             >  6:312668/1‑95 (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:171368/1‑135 (MQ=255)                                                                                                 cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATg                                                             >  5:175856/1‑135 (MQ=255)                                                                                                 cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATg                                                             <  6:175856/135‑1 (MQ=255)                                                                                                 cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATg                                                             <  6:171368/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                                    <  8:206247/147‑1 (MQ=255)                                                                                                              gCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGa                                    <  6:229966/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)                                                                                                                                               ccG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGt  <  4:178588/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                      GTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCG‑CCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGT  >  NC_002947/5732124‑5732413

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