Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,732,266	C→T	6.6%	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.6%	230.1 / 4.2	76	G97D (GGT→GAT)	hutI	imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base C (33/38);  new base T (3/2);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

CGCTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGC  >  NC_002947/5732119‑5732410                                                                                                                                                     |                                                                                                                                                   cgcTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATAc                                                                                                                                                    <  1:371680/148‑1 (MQ=255)  gcTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCg                                                                                                                                                  >  4:354226/1‑149 (MQ=255)    tCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGcc                                                                                                                                                >  8:57137/1‑149 (MQ=255)       ttGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGTGTGGTCACGGCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑c                                                                                                                                             >  1:269931/1‑149 (MQ=255)             aGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTgc                                                                                                                                      <  8:143493/149‑1 (MQ=255)             aGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTgc                                                                                                                                      <  2:73131/149‑1 (MQ=255)               ggccAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTgcgg                                                                                                                                    <  5:365424/147‑1 (MQ=255)               gTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGAGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTgcgg                                                                                                                                    >  8:273341/1‑149 (MQ=255)                   aGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTAACCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGAt                                                                                                                                >  1:120433/1‑149 (MQ=255)                     gCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGGTGGCGGTATCGCCagcact                                                                                                                                          >  5:43855/1‑133 (MQ=255)                     gCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGGTGGCGGTATCGCCagcact                                                                                                                                          <  6:43855/138‑6 (MQ=255)                           gCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCa                                                                                                                         >  1:336395/1‑148 (MQ=255)                             cGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAg                                                                                                                      >  7:254049/1‑149 (MQ=255)                              ggACTTGATCTCGAGGGTTGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCa                                                                                                                         <  8:284106/145‑1 (MQ=255)                              ggAATTGATCTCGAGGGTGGTCAC‑GCAATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGACGAACGGGTAGACCGCACATTGCTGGCGATACAGCCA‑CCCTTTGATGCGATTTCGGCa                                                                                                                         >  7:284106/1‑145 (MQ=255)                                     aTCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGcc                                                                                                              >  4:94717/1‑149 (MQ=255)                                       ctcGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCt                                                                                                            >  4:79398/1‑149 (MQ=255)                                        tcGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTg                                                                                                           >  1:332174/1‑149 (MQ=255)                                          gAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTg                                                                                                           >  5:346882/1‑147 (MQ=255)                                          gAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTg                                                                                                           >  5:351505/1‑147 (MQ=255)                                            gggTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACTGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAgg                                                                                                       <  8:312568/149‑1 (MQ=255)                                                  tCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑Gc                                                                                                >  5:209020/1‑149 (MQ=255)                                                  tCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑Gc                                                                                                <  3:94717/149‑1 (MQ=255)                                                  tCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACG‑CCTGCAgggcgtgag                                                                                                <  6:142003/149‑8 (MQ=255)                                                     c‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCg                                                                                             <  7:9446/149‑1 (MQ=255)                                                     c‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCg                                                                                             <  5:371585/149‑1 (MQ=255)                                                       gCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACGGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCg                                                                                             >  1:304900/1‑148 (MQ=255)                                                        ccATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGaa                                                                                           <  6:187597/149‑1 (MQ=255)                                                        ccATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGaa                                                                                           <  8:254049/149‑1 (MQ=255)                                                         cATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAAc                                                                                          >  3:155245/1‑149 (MQ=255)                                                             gcgcATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTcg                                                                                       <  6:90037/148‑1 (MQ=255)                                                             gcgcATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTcg                                                                                       <  2:336395/148‑1 (MQ=255)                                                              cgcATTAGTGCCGGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCTATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTcgcc                                                                                     <  6:197184/149‑1 (MQ=255)                                                               gcATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTcgccg                                                                                    >  4:181929/1‑149 (MQ=255)                                                                     gTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCgg                                                                              >  6:9606/1‑149 (MQ=255)                                                                     gTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCgg                                                                              <  2:304900/149‑1 (MQ=255)                                                                     gTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACAGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCgg                                                                              >  3:59071/1‑149 (MQ=255)                                                                        ccTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTg                                                                           <  6:351505/149‑1 (MQ=255)                                                                           gAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTgccg                                                                        >  2:289439/1‑149 (MQ=255)                                                                           gAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTgccg                                                                        >  3:240912/1‑149 (MQ=255)                                                                            aaCGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTgccgc                                                                       <  1:63877/149‑1 (MQ=255)                                                                             aCGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTgccgcc                                                                      <  4:240912/149‑1 (MQ=255)                                                                                               gCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGgtgt                                                     <  1:133411/148‑1 (MQ=255)                                                                                               gCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGtgtg                                                    >  5:301244/1‑149 (MQ=255)                                                                                                cGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTgg                                                   <  3:79398/149‑1 (MQ=255)                                                                                                  aacaacTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTgg                                                   >  5:247970/1‑147 (MQ=255)                                                                                                      acTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTgccgc                                                                       >  2:48340/1‑123 (MQ=255)                                                                                                      acTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTgccgc                                                                       <  1:48340/123‑1 (MQ=255)                                                                                                       cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTgg                                                   >  7:65753/1‑142 (MQ=255)                                                                                                       cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTgg                                                   <  8:65753/142‑1 (MQ=255)                                                                                                       cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGa                                            <  6:343272/149‑1 (MQ=255)                                                                                                        tCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATGGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGAAGTCGAt                                           >  4:191620/1‑149 (MQ=255)                                                                                                          gTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTCGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGAt                                           <  1:289439/148‑1 (MQ=255)                                                                                                          gTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGAt                                           <  3:367930/147‑1 (MQ=255)                                                                                                          gTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCa                                         <  5:9606/149‑1 (MQ=255)                                                                                                          gTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCa                                         <  3:181929/149‑1 (MQ=255)                                                                                                          gTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCa                                         <  7:57137/149‑1 (MQ=255)                                                                                                                  tGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCgg                                   >  5:122143/1‑147 (MQ=255)                                                                                                                    gCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGa                                 <  6:209020/147‑1 (MQ=255)                                                                                                                    gCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGa                                 <  2:269931/147‑1 (MQ=255)                                                                                                                      cGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑g                                                                                                 <  5:70085/81‑1 (MQ=255)                                                                                                                      cGCACGGGTAGCCCGCACAGTGCTGGCGATACAGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑g                                                                                                 >  6:70085/1‑81 (MQ=255)                                                                                                                         aCGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGAccc                          <  6:247970/149‑1 (MQ=255)                                                                                                                          cGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGAGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCa                         <  3:74639/149‑1 (MQ=255)                                                                                                                             gTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGc                      >  7:267950/1‑149 (MQ=255)                                                                                                                                gCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCGGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCAt                                                           >  2:105223/1‑109 (MQ=255)                                                                                                                                gCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCGGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCAt                                                           <  1:105223/109‑1 (MQ=255)                                                                                                                                gCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGcac                    <  7:166476/148‑1 (MQ=255)                                                                                                                                 cccGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTcgccg                                                                                    >  7:341792/1‑83 (MQ=255)                                                                                                                                 cccGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTcgccg                                                                                    <  8:341792/83‑1 (MQ=255)                                                                                                                                  ccGCACAG‑GGTGGCGGTATCGCCAGCACTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTg                              <  6:116451/136‑1 (MQ=38)                                                                                                                                  ccGCACAG‑GGTGGCGGTATCGCCAGCACTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTg                              >  5:116451/1‑136 (MQ=38)                                                                                                                                       acagggtGGCGGTATCGCCAGCACTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAg             >  1:18323/7‑149 (MQ=255)                                                                                                                                          tGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAgg                       <  4:109438/135‑1 (MQ=255)                                                                                                                                          tGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAgg                       >  3:109438/1‑135 (MQ=255)                                                                                                                                               gCGATACCGCCA‑GCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGt    <  4:59071/149‑1 (MQ=255)                                                                                                                                                 gATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGc  >  3:344334/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CGCTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGC  >  NC_002947/5732119‑5732410

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