Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,732,255	Δ1 bp	5.6%	coding (301/1206 nt)	hutI ←	imidazolonepropionase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,732,255	0	T	.	5.6%	229.5 / 6.4	72	coding (301/1206 nt)	hutI	imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base . (2/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGC  >  NC_002947/5732118‑5732390                                                                                                                                           |                                                                                                                                         gcgcTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATAc                                                                                                                               >  3:269654/1‑149 (MQ=255)  cgcTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATAc                                                                                                                               <  7:108472/148‑1 (MQ=255)    cTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGc                                                                                                                            >  2:8202/1‑149 (MQ=255)       gTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑c                                                                                                                        >  6:289513/1‑149 (MQ=255)       gTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑c                                                                                                                        >  5:217671/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)                   cAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGc                                                                                                                                        <  6:223510/122‑1 (MQ=255)                   cAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGc                                                                                                                                        >  5:223510/1‑122 (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                                                                     >  1:324078/1‑149 (MQ=255)                                                          aTCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACt                                                                     >  3:194421/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                           <  6:30617/149‑1 (MQ=255)                                                                                                     caacTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGt                           <  8:32764/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                               <  6:90198/142‑1 (MQ=255)                                                                                                        cTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTgg                               >  5:90198/1‑142 (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         <  4:281885/149‑1 (MQ=255)                                                                                                                       cGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGa         >  2:200843/1‑149 (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                                           <  4:94783/103‑1 (MQ=255)                                                                                                                                   ccGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACAcc                                           >  3:94783/1‑103 (MQ=255)                                                                                                                                           |                                                                                                                                         GCGCTCGTTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGC  >  NC_002947/5732118‑5732390

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