Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,732,273	C→A	5.3%	G95C (GGT→TGT)	hutI ←	imidazolonepropionase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,732,273	0	C	A	5.3%	243.2 / 3.0	75	G95C (GGT→TGT)	hutI	imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base C (32/39);  new base A (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 = 6.48e-01

TTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTC  >  NC_002947/5732125‑5732420                                                                                                                                                       |                                                                                                                                                    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                                                                                                                                             <  6:90453/145‑1 (MQ=38)             aGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAG‑GGTGGCGGTATCGCCAGCACTGTgcgg                                                                                                                                             >  5:90453/1‑145 (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                                                                                                                            <  3:62484/149‑1 (MQ=255)                          aCTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTc                                                                                                                            <  1:8202/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                                                                                                                 >  2:78038/1‑149 (MQ=255)                                     aGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAg                                                                                                                 >  6:72357/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                                       <  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                                                                         >  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)                                                                                                                                                       ccTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTc  >  8:308215/1‑148 (MQ=255)                                                                                                                                                       |                                                                                                                                                    TTGGCCAGGTCCAGGCCGTAGCCGGACTTGATCTCGAGGGTGGTCACGCCATCGCGCATTAGTGCCTG‑AACGCGCTGGTGAGCACTGGCGAACAACTCGTCTTCGCTGGCCGCACGGGTAGCCCGCACAGTGCTGGCGATACCGCCA‑CCCTGTGCGGCGATTTCGGCATAGCTCACGCCCTGCAGGCGTTGCTCGAACTCGCCGCTGCGGTTGCCGCCGAACACCGCATGGGTGTGGCAGTCGATCAGCCCGGGAGTGACCCAGGCACCACCCAGGTCAACCGTGCGGTCAGCTTC  >  NC_002947/5732125‑5732420

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