Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,732,257	C→G	5.1%	S100T (AGC→ACC)	hutI ←	imidazolonepropionase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,732,257	0	C	G	5.1%	257.6 / 3.2	79	S100T (AGC→ACC)	hutI	imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base C (33/42);  new base G (2/2);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

AGCATCTTGCGCTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCAC  >  NC_002947/5732110‑5732392                                                                                                                                                     |                                                                                                                                          agcaTCTTGCGCTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCt                                                                                                                                          <  7:24289/149‑1 (MQ=255)         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)                           cAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTgg                                                                                                                                        >  5:182540/1‑125 (MQ=255)                           cAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTgg                                                                                                                                        <  6:182540/125‑1 (MQ=255)                            aGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTAACCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGAt                                                                                                              >  1:120433/1‑149 (MQ=255)                              gCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGGTGGCGGTATCGCCagcact                                                                                                                        <  6:43855/138‑6 (MQ=255)                              gCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGGTGGCGGTATCGCCagcact                                                                                                                        >  5:43855/1‑133 (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                                                                           <  5:371585/149‑1 (MQ=255)                                                              c‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCg                                                                           <  7:9446/149‑1 (MQ=255)                                                                gCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACGGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCg                                                                           >  1:304900/1‑148 (MQ=255)                                                                 ccATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGaa                                                                         <  8:254049/149‑1 (MQ=255)                                                                 ccATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGaa                                                                         <  6:187597/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                                                            <  2:304900/149‑1 (MQ=255)                                                                              gTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCgg                                                            >  6:9606/1‑149 (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                                                     <  1:48340/123‑1 (MQ=255)                                                                                                               acTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTgccgc                                                     >  2:48340/1‑123 (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               <  2:269931/147‑1 (MQ=255)                                                                                                                             gCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGa               <  6:209020/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                                         <  1:105223/109‑1 (MQ=255)                                                                                                                                         gCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCGGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCAt                                         >  2:105223/1‑109 (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)                                                                                                                                                   tGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAgg     >  3:109438/1‑135 (MQ=255)                                                                                                                                                   tGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAgg     <  4:109438/135‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                          AGCATCTTGCGCTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCAC‑GCCATCGCGCATTAGTGCCTGAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTT‑GCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCAC  >  NC_002947/5732110‑5732392

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