Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,734,779	G→C	8.0%	T295T (ACC→ACG)	hutH ←	histidine ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,734,779	0	G	C	8.0%	250.9 / 8.0	82	T295T (ACC→ACG)	hutH	histidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base G (38/37);  new base C (5/2);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGGCCAGGTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCCGCGTTGGC  >  NC_002947/5734631‑5734927                                                                                                                                                     |                                                                                                                                                     ccagggccaggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTgg                                                                                                                                                      >  1:130018/1‑149 (MQ=255)      gccaggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTTACGCAGCTGGGTcag                                                                                                                                                 <  2:96324/149‑1 (MQ=255)      gccaggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcag                                                                                                                                                 <  2:293077/149‑1 (MQ=255)         aggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcaggca                                                                                                                                              <  6:368980/149‑1 (MQ=255)         aggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcaggca                                                                                                                                              <  5:152223/149‑1 (MQ=255)           gTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcaggcagg                                                                                                                                            <  6:27465/149‑1 (MQ=255)            ttGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCaggcaggc                                                                                                                                           <  5:128062/149‑1 (MQ=255)              accggtggCCATGGCCGCCGGTTCGGCGTGGAAGTTGCCGCCGGGGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCAc                                                                                                                                         >  4:153921/3‑149 (MQ=255)                        aTGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTg                                                                                                                               <  2:130018/149‑1 (MQ=255)                        aTGGCCACCGGTTCCGCGTGGAAGTCGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTCTCCGCCACAGCGTTGGCTTCCATGACCAGTACCTCGGCGGCCTGACGAAGCTGCGTGATACAGGCACCCaagcccgc                                                                                                                               >  6:285201/1‑142 (MQ=255)                           gCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCgg                                                                                                                            <  2:133655/149‑1 (MQ=255)                           gCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCgg                                                                                                                            <  3:307391/149‑1 (MQ=255)                            ccACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGc                                                                                                                           <  3:322759/149‑1 (MQ=255)                            ccACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGc                                                                                                                           >  2:395271/1‑149 (MQ=255)                             cACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCTAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGGCCAGTACCTCGGCGGCCTGACGCATCTGGGTCAGGCATGCACCCATGACCTGCGGCt                                                                                                                          <  4:232633/149‑1 (MQ=255)                             cACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCt                                                                                                                          >  7:129899/1‑149 (MQ=255)                               ccGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTgg                                                                                                                        <  1:395271/149‑1 (MQ=255)                               ccGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTgg                                                                                                                        <  7:9096/149‑1 (MQ=255)                                     cGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcgc                                                                                                                  >  7:31615/1‑149 (MQ=255)                                      ggCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcgc                                                                                                                  <  7:354516/148‑1 (MQ=255)                                      ggCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcg                                                                                                                   >  6:57321/1‑147 (MQ=255)                                       gCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcgcag                                                                                                                <  7:8413/149‑1 (MQ=255)                                          tGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCg                                                                                                                             >  8:303250/1‑133 (MQ=255)                                            gAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAg                                                                                                            >  4:185116/1‑148 (MQ=255)                                            gAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAg                                                                                                            <  2:324657/148‑1 (MQ=255)                                                ttGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGTCCTGCGGCTGGCAGCGCAGCGAGTACgg                                                                                                       <  5:334465/149‑1 (MQ=255)                                                    cgccgGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTcc                                                                                                   >  5:47896/1‑149 (MQ=255)                                                        ggAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGa                                                                                               >  5:229262/1‑149 (MQ=255)                                                         gagaTCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGAc                                                                                              >  6:245360/1‑149 (MQ=255)                                                          agaTCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGAcc                                                                                             >  3:128046/1‑149 (MQ=255)                                                              cacatcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGt                                                                                         >  5:406856/1‑149 (MQ=255)                                                               acatcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTc                                                                                        >  2:30253/1‑149 (MQ=255)                                                                catcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCt                                                                                                  <  5:218489/138‑1 (MQ=255)                                                                catcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCt                                                                                                  >  6:218489/1‑138 (MQ=255)                                                                 atcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCg                                                                                       >  2:118211/1‑148 (MQ=255)                                                                  tcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCg                                                                                       >  6:128751/1‑147 (MQ=255)                                                                     ccTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGAGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGAAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGt                                                                                   >  2:315464/1‑148 (MQ=255)                                                                            cAGCGAACACCGGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgc                                                                           <  2:232835/149‑1 (MQ=255)                                                                             aGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgcg                                                                          <  5:245360/149‑1 (MQ=255)                                                                                 aaCACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgcgacag                                                                      >  4:239191/1‑149 (MQ=255)                                                                                      cAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgcgacagcgaca                                                                 >  4:145263/1‑149 (MQ=255)                                                                                      cAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgcgacagcgaca                                                                 >  3:272371/1‑149 (MQ=255)                                                                                      cAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgcgacagcgac                                                                  >  7:82726/1‑148 (MQ=255)                                                                                        gCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACt                                                               <  1:30253/149‑1 (MQ=255)                                                                                             ttGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCt                                                          <  6:47896/149‑1 (MQ=255)                                                                                               gTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                                         <  1:93617/148‑1 (MQ=255)                                                                                                tCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGa                                                       <  5:128751/149‑1 (MQ=255)                                                                                                  cGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCAGCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGt                                                     >  5:117408/1‑149 (MQ=255)                                                                                                     cacaGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGc                                                  <  5:52236/149‑1 (MQ=255)                                                                                                       caGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGccc                                                <  6:229262/149‑1 (MQ=255)                                                                                                          cGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                                         >  6:172909/1‑137 (MQ=255)                                                                                                          cGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                                         <  5:172909/137‑1 (MQ=255)                                                                                                          cGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcagc                                             <  8:82726/149‑1 (MQ=255)                                                                                                          cGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcag                                              >  1:183766/1‑148 (MQ=255)                                                                                                           gTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                                         <  1:37736/136‑1 (MQ=255)                                                                                                           gTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                                         >  2:37736/1‑136 (MQ=255)                                                                                                           gTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcag                                              >  1:320588/1‑147 (MQ=255)                                                                                                            ttGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCcagcag                                           >  1:275689/1‑149 (MQ=255)                                                                                                               gCTTCGATGCCCAGTACCTCGGCGGCCTGGCGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTc                                        >  8:296498/1‑149 (MQ=255)                                                                                                                    gATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCgg                                    >  4:158681/1‑148 (MQ=255)                                                                                                                      tGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAg                                 <  8:129899/149‑1 (MQ=255)                                                                                                                          cAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAgg                             <  4:128046/149‑1 (MQ=255)                                                                                                                          cAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAgg                             >  3:55700/1‑149 (MQ=255)                                                                                                                               ccTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCcgcg                        <  5:57321/149‑1 (MQ=255)                                                                                                                                cTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAg                                                                                    <  8:270752/88‑1 (MQ=255)                                                                                                                                cTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAg                                                                                    >  7:270752/1‑88 (MQ=255)                                                                                                                                cTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCgg                       >  6:153488/1‑149 (MQ=255)                                                                                                                                 tCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGgt                      >  6:247103/1‑149 (MQ=255)                                                                                                                                 tCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGgt                      >  6:138687/1‑149 (MQ=255)                                                                                                                                   ggcggcCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGgt                      <  3:239191/147‑1 (MQ=255)                                                                                                                                     ctgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGACAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTcgc                                      <  4:257769/127‑1 (MQ=255)                                                                                                                                     ctgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGACAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTcgc                                      >  3:257769/3‑129 (MQ=255)                                                                                                                                      tgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGa                                   >  5:97340/2‑131 (MQ=255)                                                                                                                                      tgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGa                                   <  6:97340/130‑1 (MQ=255)                                                                                                                                      tgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGa                 >  6:259819/2‑149 (MQ=255)                                                                                                                                      tgcCTGACCCAGCTGCGTCAGGCAGGAACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCg                  >  3:402896/2‑148 (MQ=255)                                                                                                                                      ggcCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGa                 <  3:153921/149‑1 (MQ=255)                                                                                                                                      ggcCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGa                 <  3:158681/149‑1 (MQ=255)                                                                                                                                        ccTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGGAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATc               >  1:128977/1‑149 (MQ=255)                                                                                                                                        ccTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATc               >  1:274969/1‑149 (MQ=255)                                                                                                                                                aGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCcgcg       >  5:339780/1‑149 (MQ=255)                                                                                                                                                aGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCcgcg       >  5:334211/1‑149 (MQ=255)                                                                                                                                                 gCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCCGCGt      <  2:275689/149‑1 (MQ=255)                                                                                                                                                     ggTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGGAGGCCGCGGTGTCGATCTGGCCGCGTTg    <  3:145263/147‑1 (MQ=255)                                                                                                                                                     ggTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCCGCGTTgg   <  2:128977/148‑1 (MQ=255)                                                                                                                                                     ggTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCCGCGTTGGc  <  5:153488/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCAGGGCCAGGTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCCGCGTTGGC  >  NC_002947/5734631‑5734927

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