BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,732,026	A→T	5.8%	L196M (TTG→ATG) ‡	hutI ←	imidazolonepropionase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,732,026	0	A	T	5.8%	321.5 / 8.9	103	L196M (TTG→ATG) ‡	hutI	imidazolonepropionase
Reads supporting (aligned to +/- strand): ref base A (50/47); new base T (3/3); total (53/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

Only 100 of 103 total aligned reads displayed.

CCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCTT  >  NC_002947/5731878‑5732173

ccGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAg                                                                                                                                                    <  1:482299/150‑1 (MQ=255)

  gTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATAAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGc                                                                                                                                                  >  2:1225807/1‑150 (MQ=255)

   tGCAACGATGACAGTTGCTCGGAATCCATCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGcc                                                                                                                                                 <  2:1160371/150‑1 (MQ=255)

         gATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAAc                                                                                                                                           <  2:2396073/150‑1 (MQ=255)

          aTGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCTTGGCCGGcgaag                                                                                                                                          <  1:995673/150‑5 (MQ=255)

           tGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGGAACAt                                                                                                                                         <  2:94997/150‑1 (MQ=255)

           tGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACAt                                                                                                                                         <  1:282439/150‑1 (MQ=255)

            gACAGTTGTTCGGCATGCAGCTTGACTCAGACGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACAtc                                                                                                                                        <  1:1811325/150‑1 (MQ=255)

             aCAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCTTGGCCGGcgaa                                                                                                                                           <  2:103488/146‑4 (MQ=255)

              cAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACAtctc                                                                                                                                      <  2:660596/150‑1 (MQ=255)

               aGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCg                                                                                                                                     >  2:1951295/1‑150 (MQ=255)

               aGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAATTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCg                                                                                                                                     >  2:234588/1‑150 (MQ=255)

                    ttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTcaca                                                                                                                                <  1:1729829/150‑1 (MQ=255)

                    ttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTcaca                                                                                                                                >  2:2239624/1‑150 (MQ=255)

                    ttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTcaca                                                                                                                                >  1:682362/1‑150 (MQ=255)

                        gCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAtg                                                                                                                            >  1:1204998/1‑150 (MQ=255)

                           tGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTgc                                                                                                                         >  1:1776960/1‑150 (MQ=255)

                           tGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTgc                                                                                                                         <  2:682362/150‑1 (MQ=255)

                              aGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGcgcg                                                                                                                      >  2:1258572/1‑150 (MQ=255)

                                cTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCg                                                                                                                                                <  2:1183763/122‑1 (MQ=255)

                                cTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCg                                                                                                                                                >  1:1183763/1‑122 (MQ=255)

                                   gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACCGATGTGCGCGAtacc                                                                                                                 <  2:33849/150‑4 (MQ=255)

                                        gCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAgtcgt                                                                                                            <  1:2254212/150‑1 (MQ=255)

                                        gCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAgtcgt                                                                                                            <  1:471580/150‑1 (MQ=255)

                                         ccggcccaggtcgagcggctgttcatcaaggcgcgcgaactgggccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAgtcgtc                                                                                                           >  1:211538‑M2/47‑150 (MQ=255)

                                            gCCCAGTTCGCGCGCCTTGATGAACAGCAGAACGACGTGAGATAGCGAGCGGACGAGGTGTTCGCAGAAGGTAGCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACAACTTGTCACATATGTGCGCGATGTAGTCGTCGGc                                                                                                        >  1:2224318/1‑150 (MQ=255)

                                            GCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGc                                                                                                        <  2:2215383‑M1/150‑1 (MQ=255)

                                              CCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGccc                                                                                                      <  2:135841‑M1/150‑1 (MQ=255)

                                               CAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAgtcgt                                                                                                            <  1:2002535‑M1/143‑1 (MQ=255)

                                               CAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAgtcgt                                                                                                            >  2:2002535‑M1/1‑143 (MQ=255)

                                                aggtcgagcggctgttcatcaaggcgcgcgaactgggccTGCGAGCACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCcgt                                                                                                    >  2:88560‑M2/40‑149 (MQ=255)

                                                aggtcgagcggctgttcatcaaggcgcgcgaactgggccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGAt                                                                                                                    >  1:2235979‑M2/40‑134 (MQ=255)

                                                aggtcgagcggctgttcatcaaggcgcgcgaactgggccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGAt                                                                                                                    <  2:2235979‑M2/95‑1 (MQ=255)

                                                AGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGa                                                                                                    <  1:2358078‑M1/150‑1 (MQ=255)

                                                 GTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGAc                                                                                                   >  2:1976509‑M1/1‑150 (MQ=255)

                                                  ttCGCGCGCCTTGATGAAAAGCCGCTCGACCTGGGCCAGCGAGAACGCGAGGTGTTCGCACAAGGCATCCACCGCATCCACAAGGCCTTCGCCGGCCAAGACCGGCAACATCTCGTCCCAGATGTGCGCGATGTAGACGTCGGCCCGAcc                                                                                                  >  2:708248/1‑150 (MQ=255)

                                                     gcgcgcCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTTGCCGGCCTTGGCCGGcgaag                                                                                                                                          <  1:1786400/107‑5 (MQ=255)

                                                     gcgcgcCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTTGCCGGCCTTGGCCGGcgaag                                                                                                                                          >  2:1786400/1‑103 (MQ=255)

                                                         gcCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTAc                                                                                           >  2:1091531/1‑150 (MQ=255)

                                                          ccTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACt                                                                                          >  1:7934/1‑150 (MQ=255)

                                                           cTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTc                                                                                         >  2:1279242/1‑150 (MQ=255)

                                                              gatgaACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCgg                                                                                      >  1:777121/1‑150 (MQ=255)

                                                               atgaACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCggg                                                                                     <  1:2215383/150‑1 (MQ=255)

                                                               atgaACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCCCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCggg                                                                                     <  2:1854029/150‑1 (MQ=255)

                                                                   aCAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTgc                                                                                                                         >  2:1031470/1‑110 (MQ=255)

                                                                   aCAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTgc                                                                                                                         <  1:1031470/110‑1 (MQ=255)

                                                                    cAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAgtcgt                                                                                                            >  1:508758/1‑122 (MQ=255)

                                                                    cAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAgtcgt                                                                                                            <  2:508758/122‑1 (MQ=255)

                                                                     aGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCaaa                                                                               <  1:1225807/150‑1 (MQ=255)

                                                                     aGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCaaa                                                                               >  2:86141/1‑150 (MQ=255)

                                                                      gCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAg                                                                              >  1:498104/1‑150 (MQ=255)

                                                                       ccGCTCGCCCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTGGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGc                                                                             <  2:2314039/150‑1 (MQ=255)

                                                                         gCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTca                                                                                                                                  >  1:1118505/1‑95 (MQ=255)

                                                                         gCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTca                                                                                                                                  <  2:1118505/95‑1 (MQ=255)

                                                                                ctgggccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGcgccgc                                                                    >  1:1068255‑M2/8‑150 (MQ=255)

                                                                                 tgggccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCCAGATCTCCTCACAGCGGTGCTCGACGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCgccgcc                                                                   <  1:2101942‑M2/144‑1 (MQ=255)

                                                                                   ggccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTgcgccgccgg                                                                 <  2:2263256‑M2/146‑3 (MQ=255)

                                                                                   ggccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAg                                                                 <  2:1796636‑M2/146‑1 (MQ=255)

                                                                                   ggccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAg                                                                 >  1:1286486‑M2/5‑150 (MQ=255)

                                                                                   ggccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAg                                                                 <  2:2224318‑M2/146‑1 (MQ=255)

                                                                                     ccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGc                                                               >  1:1465082‑M2/3‑150 (MQ=255)

                                                                                     ccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGc                                                               >  1:142518‑M2/3‑150 (MQ=255)

                                                                                         ggagaaCGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGT                                                           <  2:1660541‑M1/149‑1 (MQ=255)

                                                                                           agaACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGG                                                         >  1:2013792‑M1/1‑150 (MQ=255)

                                                                                           agaACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGG                                                         <  1:2065040‑M1/150‑1 (MQ=255)

                                                                                            gaACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGC                                                        <  2:1286486‑M1/150‑1 (MQ=255)

                                                                                             aaCGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGgca                                                       <  1:155235/150‑1 (MQ=255)

                                                                                              aCGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGgcac                                                      <  1:234588/150‑1 (MQ=255)

                                                                                                 cgAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGgcacgca                                                   >  2:559426/1‑150 (MQ=255)

                                                                                                       gtTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTa                                                                                            >  1:778585/1‑103 (MQ=255)

                                                                                                       gtTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTa                                                                                            <  2:778585/103‑1 (MQ=255)

                                                                                                        ttCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAg                                            >  2:957598/1‑150 (MQ=255)

                                                                                                           gCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCgg                                         >  1:1444124/1‑150 (MQ=255)

                                                                                                                    ctccaccgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgtc                                <  1:1057519/149‑1 (MQ=255)

                                                                                                                    atccaccgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgtc                                <  1:1258572/150‑1 (MQ=255)

                                                                                                                    atccaccgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgtc                                <  2:142518/150‑1 (MQ=255)

                                                                                                                    agccaccgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCGAGGCGTCGGCCGCCGGGCACGTGGCACGCACCGTCAGGGGCAGCtcgtc                                >  1:2154107/3‑150 (MQ=255)

                                                                                                                     tccaccgcatgcaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgtcg                               >  1:1197356/1‑150 (MQ=255)

                                                                                                                     tccaccgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgtcg                               >  2:2012785/1‑150 (MQ=255)

                                                                                                                      ccaccgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCgg                              >  1:956823/1‑150 (MQ=255)

                                                                                                                        accgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGCCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGcc                            >  1:2253239/1‑150 (MQ=255)

                                                                                                                          cgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGACGTAGTCGTCGGCCCGACCAGCGTACTCCGGGGGCAAAGCGTGCGCCGCCTGGCACGTGGCACGCACCGTCAGCGGAAGCGCCTAGGCCAg                          >  1:885637/1‑150 (MQ=255)

                                                                                                                           gcatccaccAGGCCTTCGCCGGCCAATGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGGCTTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCCCCGTCAGCGGCAGCTCGTCGGCCAgg                         >  2:490502/1‑150 (MQ=255)

                                                                                                                            catccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAggc                        <  1:559426/150‑1 (MQ=255)

                                                                                                                                caccaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGgcggcg                    <  1:88560/150‑1 (MQ=255)

                                                                                                                                 accaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGcggcgg                   <  2:2065500/150‑1 (MQ=255)

                                                                                                                                 accaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGcggcgg                   >  1:2065500/1‑150 (MQ=255)

                                                                                                                                  ccaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCccgccgcctggccgacgagctgccgctgacggtgcgtg                  >  1:818530‑M1/1‑112 (MQ=255)

                                                                                                                                  ccaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCccgccgcctggccgacgagctgccgctgacggtgcgtg                  <  1:1456263‑M1/150‑39 (MQ=255)

                                                                                                                                  ccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGg                  >  2:1444124‑M2/1‑149 (MQ=255)

                                                                                                                                    aGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCccgccgcctggccgacgagctgccgctgacggtgcgtgcc                <  2:818530‑M1/150‑41 (MQ=255)

                                                                                                                                        cTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTc            >  1:888524/1‑150 (MQ=255)

                                                                                                                                         ttCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCaaa                                                                               <  1:1848745/82‑1 (MQ=255)

                                                                                                                                         ttCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCaaa                                                                               >  2:1848745/1‑82 (MQ=255)

                                                                                                                                          ttgCCGGCCTTGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCgc          >  1:413128/3‑150 (MQ=255)

                                                                                                                                          tCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGGGCGCCGCCAGGCACGTGGCccgccgcctggccgacgagctgccgctgacggtgcgtgccacgtgc          <  2:1407318‑M1/150‑47 (MQ=255)

                                                                                                                                            gCCGGCCTTGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCgcag        >  2:2000631/1‑150 (MQ=255)

                                                                                                                                               ggCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCATCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGCCGGCCAGGCGGCGGGCCACTCGCAGCAt     <  1:1530851/150‑1 (MQ=255)

                                                                                                                                                 ccAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCt   >  2:2052301/1‑150 (MQ=255)

                                                                                                                                                  cAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCtt  >  1:1448307/1‑150 (MQ=255)

CCGTGCAACGATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCTT  >  NC_002947/5731878‑5732173

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`