Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,732,271:1	+G	5.3%	coding (285/1206 nt)	hutI ←	imidazolonepropionase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,732,271	1	.	G	5.3%	312.1 / 8.2	75	coding (285/1206 nt)	hutI	imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base . (32/39);  new base G (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

GTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGT  >  NC_002947/5732124‑5732413                                                                                                                                                      |                                                                                                                                               gTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑c                                                                                                                                               >  5:217671/1‑149 (MQ=255) gTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑c                                                                                                                                               >  6:289513/1‑149 (MQ=255)  ttGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑cc                                                                                                                                              <  1:21600/149‑1 (MQ=255)     gCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCAGCACAGTGCTGGCGATACCGCCA‑CCCtg                                                                                                                                           >  1:293805/1‑149 (MQ=255)      ccAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCtgt                                                                                                                                          >  3:314097/1‑149 (MQ=255)        aGGTCCAGGACGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTgc                                                                                                                                        >  5:93057/1‑149 (MQ=255)          gTCCAGGACGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTgcgg                                                                                                                                      <  6:93057/149‑1 (MQ=255)             cAGACCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGa                                                                                                                                   >  2:47407/1‑149 (MQ=255)              aGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAG‑GGTGGCGGTATCGCCAGCACTGTgcgg                                                                                                                                      >  5:90453/1‑145 (MQ=38)              aGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAG‑GGTGGCGGTATCGCCAGCACTGTgcgg                                                                                                                                      <  6:90453/145‑1 (MQ=38)                gCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGAttt                                                                                                                                <  4:85314/149‑1 (MQ=255)                   gTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCgg                                                                                                                             <  1:47407/149‑1 (MQ=255)                        cGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTGTAACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAG‑GGTGGCGGTATCGCCAGCACTGTgcgg                                                                                                                                      >  1:74067/1‑136 (MQ=37)                        cGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAG‑GGTGGCGGTATCGCCAGCACTGTgcgg                                                                                                                                      <  2:74067/135‑1 (MQ=38)                           aCTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCt                                                                                                                      >  7:119191/1‑148 (MQ=255)                           aCTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTc                                                                                                                     <  1:8202/149‑1 (MQ=255)                           aCTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTc                                                                                                                     <  3:62484/149‑1 (MQ=255)                            cTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCa                                                                                                                    <  5:289513/149‑1 (MQ=255)                             ttGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCAc                                                                                                                   >  7:129382/1‑149 (MQ=255)                             ttGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCAc                                                                                                                   >  5:77342/1‑149 (MQ=255)                               gATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCAc                                                                                                                   <  6:77342/147‑1 (MQ=255)                               gATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACg                                                                                                                  <  4:269654/148‑1 (MQ=255)                                      aGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAg                                                                                                          >  6:72357/1‑149 (MQ=255)                                      aGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAg                                                                                                          >  2:78038/1‑149 (MQ=255)                                       gggTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAgg                                                                                                         >  7:159965/1‑149 (MQ=255)                                          tggtCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGt                                                                                                      <  5:72357/149‑1 (MQ=255)                                             tCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGc                                                                                                   <  5:298376/149‑1 (MQ=255)                                                cGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTc                                                                                                 >  8:226369/1‑148 (MQ=255)                                                 gCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGc                                                                                                   <  1:184811/145‑1 (MQ=255)                                                 gCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGc                                                                                                   >  2:184811/1‑145 (MQ=255)                                                    aTCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACt                                                                                            >  3:194421/1‑149 (MQ=255)                                                    aTCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACt                                                                                            >  1:324078/1‑149 (MQ=255)                                                     tCGCGCGTTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTc                                                                                           >  7:247087/1‑149 (MQ=255)                                                         gcATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTcgccg                                                                                       <  8:159965/149‑1 (MQ=255)                                                                   cTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGTGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCCCCGCTGCGGTTgc                                                                             >  7:254539/1‑149 (MQ=255)                                                                       aaCGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTgccg                                                                           <  3:129899/148‑1 (MQ=255)                                                                        aCGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTgccgcc                                                                         <  8:46504/149‑1 (MQ=255)                                                                           cgcTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGaa                                                                      <  8:119191/149‑1 (MQ=255)                                                                            gcTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAac                                                                     >  5:146641/1‑149 (MQ=255)                                                                                gTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCTCTCAGTGCTGGCGATACCGCCA‑CCCTGTGAGGCGATTTCGGCATAGCTCACGGCCTGCTGGCGTTGTTCGAACTCGCCGCTGCGGTTGCCGCATAACACCg                                                                 <  2:293805/149‑1 (MQ=255)                                                                                gTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAGCTCGCCGCTGCGGTTGCCGCCGAACAcc                                                                  <  3:258215/148‑1 (MQ=255)                                                                                   aGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCAt                                                              >  6:52761/1‑149 (MQ=255)                                                                                       cTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGgt                                                          <  8:75122/149‑1 (MQ=255)                                                                                        tGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑ACCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGgtg                                                         >  2:19689/1‑149 (MQ=255)                                                                                           cGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTgg                                                      <  6:146641/149‑1 (MQ=255)                                                                                           cGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTgg                                                      <  8:129382/149‑1 (MQ=255)                                                                                               caacTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGt                                                  <  8:32764/149‑1 (MQ=255)                                                                                               caacTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGt                                                  <  6:30617/149‑1 (MQ=255)                                                                                                 acTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGgtgt                                                        <  8:183810/141‑1 (MQ=255)                                                                                                 acTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGgtgt                                                        >  7:183810/1‑141 (MQ=255)                                                                                                 acTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCg                                                <  6:217671/149‑1 (MQ=255)                                                                                                  cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTgg                                                      >  5:90198/1‑142 (MQ=255)                                                                                                  cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTgg                                                      <  6:90198/142‑1 (MQ=255)                                                                                                     gTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGc                                                                <  5:320765/129‑1 (MQ=255)                                                                                                     gTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGc                                                                >  6:320765/1‑129 (MQ=255)                                                                                                          cGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCg                                       >  1:165578/1‑149 (MQ=255)                                                                                                              ggCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCgg                                      >  6:9357/1‑146 (MQ=255)                                                                                                              ggCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGa                                    <  5:9357/148‑1 (MQ=255)                                                                                                               gCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGa                                    <  4:194421/147‑1 (MQ=255)                                                                                                                 cGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGa                                >  2:200843/1‑149 (MQ=255)                                                                                                                 cGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGa                                <  4:281885/149‑1 (MQ=255)                                                                                                                    aCGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGAcc                              >  5:286382/1‑148 (MQ=255)                                                                                                                       ggTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGAGTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAgg                          <  8:273628/149‑1 (MQ=255)                                                                                                                        gTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCCCGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAgg                          >  5:94463/1‑148 (MQ=255)                                                                                                                        gTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGc                         <  8:247087/149‑1 (MQ=255)                                                                                                                             ccGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACAcc                                                                  >  3:94783/1‑103 (MQ=255)                                                                                                                             ccGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACAcc                                                                  <  4:94783/103‑1 (MQ=255)                                                                                                                                        tGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAAcc         >  3:218240/1‑149 (MQ=255)                                                                                                                                           cGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTg      <  7:226369/149‑1 (MQ=255)                                                                                                                                           cGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTg      <  2:165578/149‑1 (MQ=255)                                                                                                                                            gATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCaa           >  5:166039/1‑143 (MQ=255)                                                                                                                                            gATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCaa           <  6:166039/143‑1 (MQ=255)                                                                                                                                             aTACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCg    <  6:94463/149‑1 (MQ=255)                                                                                                                                             aTACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCg    <  4:132470/149‑1 (MQ=255)                                                                                                                                             aTACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGAAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCg    <  1:78038/149‑1 (MQ=255)                                                                                                                                               aCCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGt  >  2:166812/1‑149 (MQ=255)                                                                                                                                                ccGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGa                                                                       <  3:73380/79‑1 (MQ=255)                                                                                                                                                ccGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGa                                                                       >  4:73380/1‑79 (MQ=255)                                                                                                                                                      |                                                                                                                                               GTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑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: ‑