Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,732,263	A→G	6.3%	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.3%	237.5 / 4.2	79	I98T (ATC→ACC)	hutI	imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base A (33/41);  new base G (3/2);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GCGCTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGC  >  NC_002947/5732118‑5732410                                                                                                                                                   |                                                                                                                                                      tcgcTCGTTGCCCAGGTCCAGGCCGTAGCCTGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGAt                                                                                                                                                      <  6:262459/146‑1 (MQ=255) gcgcTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGAt                                                                                                                                                      <  6:262455/147‑1 (MQ=255)  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                                                                                                                                      <  2:73131/149‑1 (MQ=255)              aGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTgc                                                                                                                                      <  8:143493/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:351505/1‑147 (MQ=255)                                           gAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTg                                                                                                           >  5:346882/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                                                                                       <  2:336395/148‑1 (MQ=255)                                                              gcgcATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTcg                                                                                       <  6:90037/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                                                                        >  3:240912/1‑149 (MQ=255)                                                                            gAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTgccg                                                                        >  2:289439/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                                                   <  8:65753/142‑1 (MQ=255)                                                                                                        cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTgg                                                   >  7:65753/1‑142 (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                              >  5:116451/1‑136 (MQ=38)                                                                                                                                   ccGCACAG‑GGTGGCGGTATCGCCAGCACTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTg                              <  6:116451/136‑1 (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)                                                                                                                                                   |                                                                                                                                                      GCGCTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGC  >  NC_002947/5732118‑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: ‑