BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,116,278	C→T	7.9%	Q315Q (CAG→CAA) ‡	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,278	0	C	T	7.9%	308.6 / 10.5	101	Q315Q (CAG→CAA) ‡	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (36/57); new base T (4/4); total (40/61)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 102 total aligned reads displayed.

CGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATA  >  NC_000913/4116131‑4116427

cGCCCCGCTCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTTGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTg                                                                                                                                                     <  2:1404697/150‑1 (MQ=255)

cGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTg                                                                                                                                                     <  1:1310438/150‑1 (MQ=255)

 gCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGa                                                                                                                                                    >  1:1059592/1‑150 (MQ=255)

  ccccGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGaa                                                                                                                                                   <  1:672907/150‑1 (MQ=255)

    ccGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATg                                                                                                                                                 <  2:952906/150‑1 (MQ=255)

     cgcgcgCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTGATTTCATCGCGCAGCCACTGAATgg                                                                                                                                                <  2:1605428/150‑1 (MQ=255)

            aTACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGcc                                                                                                                                         <  2:1091449/150‑1 (MQ=255)

             tACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCt                                                                                                                                        >  2:734788/1‑150 (MQ=255)

             tACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCt                                                                                                                                        >  2:704548/1‑150 (MQ=255)

               cGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGc                                                                                                                                      >  2:534452/1‑150 (MQ=255)

               cGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCACATCCATTCAGTGGCTgcgcgatg                                                                                                                                      <  1:272919/150‑5 (MQ=255)

                 ggTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCa                                                                                                                                    <  1:197495/150‑1 (MQ=255)

                 ggTCACAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCa                                                                                                                                    >  2:1179528/1‑150 (MQ=255)

                  gTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCATCCATTCAGTGGCTgcgcgatgaaa                                                                                                                                   <  2:1225808/150‑8 (MQ=255)

                  gTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCATCCATTCAGTGGCTgcgcgatgaaa                                                                                                                                   <  1:1340086/150‑8 (MQ=255)

                  gTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCATCCATTCAGTGGCTgcgcgatgaaa                                                                                                                                   <  1:1491983/150‑8 (MQ=255)

                    cccAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATaa                                                                                                                                 >  1:112050/1‑150 (MQ=255)

                          aCGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCg                                                                                                                           <  1:802616/150‑1 (MQ=255)

                           cGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGc                                                                                                                          <  2:283444/150‑1 (MQ=255)

                           cGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGc                                                                                                                          <  2:669908/150‑1 (MQ=255)

                            ggCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGcc                                                                                                                         >  1:419449/1‑149 (MQ=255)

                             gcgcACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCAc                                                                                                                        <  1:1197865/150‑1 (MQ=255)

                                 aCCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTc                                                                                                                    >  2:228229/1‑150 (MQ=255)

                                  cccATCCCGGTAAATGCCGGAACCACATACACACCCTTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGACTCGTAGGCGTCGTTAATCAACTTCATTTCCTCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTcc                                                                                                                   <  2:1085171/150‑1 (MQ=255)

                                  cccAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTcc                                                                                                                   <  1:734788/150‑1 (MQ=255)

                                   ccAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTc                                                                                                                    <  1:1135218/148‑1 (MQ=255)

                                   ccAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTc                                                                                                                    >  2:1135218/1‑148 (MQ=255)

                                   ccAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTc                                                                                                                    >  2:1140840/1‑148 (MQ=255)

                                   ccAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTc                                                                                                                    <  1:1140840/148‑1 (MQ=255)

                                      gCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTTGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACg                                                                                                               <  1:1025111/150‑1 (MQ=255)

                                        ccGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCa                                                                                                             >  1:838651/1‑150 (MQ=255)

                                              aaTGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTc                                                                                                       >  1:1014287/1‑150 (MQ=255)

                                              aaTGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTc                                                                                                       <  2:486573/150‑1 (MQ=255)

                                               aTGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCa                                                                                                      >  1:245966/1‑150 (MQ=255)

                                               aTGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCa                                                                                                      <  1:1505622/150‑1 (MQ=255)

                                                 gCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACt                                                                                                    >  1:605284/1‑150 (MQ=255)

                                                            aTACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATGCGTAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCgg                                                                                         <  2:434709/150‑1 (MQ=255)

                                                              acacacCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCTACGCATAGTTCACTTCGCCAGTCGGGc                                                                                       <  1:624459/150‑1 (MQ=255)

                                                                 cacCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGc                                                                                    <  2:567433/150‑1 (MQ=255)

                                                                 cacCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGc                                                                                    >  2:298041/1‑150 (MQ=255)

                                                                  acCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCa                                                                                   <  1:1576945/150‑1 (MQ=255)

                                                                         gtgtTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtg                                                                            <  1:114424/150‑1 (MQ=255)

                                                                          ggtTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGGCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtgg                                                                           <  1:1594057/149‑1 (MQ=255)

                                                                           gtTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtggt                                                                          <  1:1426505/150‑1 (MQ=255)

                                                                            ttttGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtggtg                                                                         >  1:443724/1‑150 (MQ=255)

                                                                            ttttGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtggtg                                                                         <  1:1492880/150‑1 (MQ=255)

                                                                             tttGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtggtgg                                                                        <  1:941880/150‑1 (MQ=255)

                                                                             tttGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGAtggtgg                                                                        <  2:1014287/150‑1 (MQ=255)

                                                                              ttGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATggtggt                                                                       <  1:562057/150‑1 (MQ=255)

                                                                                gCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTca                                                                     <  2:1112289/150‑1 (MQ=255)

                                                                                   cTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCGTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCCCAGGCGATGGTGGTcagcc                                                                  >  2:846078/1‑149 (MQ=255)

                                                                                   cTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTcagca                                                                  >  1:642974/1‑150 (MQ=255)

                                                                                   cTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTcagca                                                                  >  2:1553419/1‑150 (MQ=255)

                                                                                   cTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTcagca                                                                  >  2:1239243/1‑150 (MQ=255)

                                                                                            cGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGtttt                                                         <  1:84054/150‑1 (MQ=255)

                                                                                            cGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGtttt                                                         >  1:882367/1‑150 (MQ=255)

                                                                                              aaaTATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCt                                                       <  1:1179528/150‑1 (MQ=255)

                                                                                              aaaTATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCt                                                       <  2:975565/150‑1 (MQ=255)

                                                                                              aaaTATCCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCt                                                       <  2:608005/150‑1 (MQ=255)

                                                                                               aaTATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTg                                                      <  2:1059592/150‑1 (MQ=255)

                                                                                                 tatTCGGATTCGTAGGCGTGGTTAATCATCTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCTACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGAt                                                    <  2:1454098/150‑1 (MQ=255)

                                                                                                 tatTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGAt                                                    >  2:1445729/1‑150 (MQ=255)

                                                                                                   ttCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGAttt                                                  >  2:1436391/1‑150 (MQ=255)

                                                                                                    tCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTc                                                 >  2:132258/1‑150 (MQ=255)

                                                                                                     cGGAATCGTAGGCGTCGTTAATCAACTTCATTTCCTCGCGCAGCCACTGAAAGGATGCGCCTGCCATAACCACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGGCCATTTCTGATTTCa                                                >  2:378933/1‑150 (MQ=255)

                                                                                                       ccaTCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCAcc                                              >  1:1123833/3‑150 (MQ=255)

                                                                                                        aaTCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGACATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCg                                             >  1:357547/1‑150 (MQ=255)

                                                                                                         aTCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGc                                            <  1:1514871/150‑1 (MQ=255)

                                                                                                           cGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGTCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCtt                                          <  2:357547/150‑1 (MQ=255)

                                                                                                             tAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCAAAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTtc                                        <  1:1137928/150‑1 (MQ=255)

                                                                                                              aGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTtct                                       <  2:77229/150‑1 (MQ=255)

                                                                                                               ggCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTTTGATTTCACCGCTTtctc                                      <  2:245966/150‑1 (MQ=255)

                                                                                                                  gtcgtTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGcc                                   <  2:1326655/150‑1 (MQ=255)

                                                                                                                  gtcgtTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGcc                                   <  1:529571/150‑1 (MQ=255)

                                                                                                                    cgtTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAg                                 >  1:1123905/1‑150 (MQ=255)

                                                                                                                    cgtTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAg                                 <  1:132258/150‑1 (MQ=255)

                                                                                                                    cgtTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAg                                 <  1:42571/150‑1 (MQ=255)

                                                                                                                      ttAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACTCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAgtg                               <  2:882367/150‑1 (MQ=255)

                                                                                                                       tAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAgtgt                              <  1:298041/150‑1 (MQ=255)

                                                                                                                       tAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAgtgt                              >  1:890393/1‑150 (MQ=255)

                                                                                                                       tAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAgtgt                              >  1:731806/1‑150 (MQ=255)

                                                                                                                        aaTCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGtt                             <  2:1123833/150‑1 (MQ=255)

                                                                                                                           cAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTtcat                          <  1:1239243/150‑1 (MQ=255)

                                                                                                                              cTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAg                       >  1:1075522/1‑150 (MQ=255)

                                                                                                                                tcatttcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCa                     <  1:704548/150‑1 (MQ=255)

                                                                                                                                tcatttcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCa                     <  1:1203211/150‑1 (MQ=255)

                                                                                                                                 catttcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCAt                    <  2:644474/150‑1 (MQ=255)

                                                                                                                                  atttcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATa                   <  1:527057/150‑1 (MQ=255)

                                                                                                                                    ttcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATaaa                 >  2:1313133/1‑150 (MQ=255)

                                                                                                                                     ggcaggcgcATCCATTCAGTGGCTGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCACAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAag                >  1:570266/6‑150 (MQ=255)

                                                                                                                                      catCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagc               >  2:560860/1‑150 (MQ=255)

                                                                                                                                       caggcgcATCCATTCAGTGGCTGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagca              >  1:1463378/4‑150 (MQ=255)

                                                                                                                                       caggcgcATCCATTCAGTGGCTGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagca              >  1:342441/4‑150 (MQ=255)

                                                                                                                                         cgcgCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagcagc            >  1:1175658/1‑150 (MQ=255)

                                                                                                                                          gcgcATCCATTCAGTGGCTGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGcc           >  2:426986/1‑150 (MQ=255)

                                                                                                                                               gCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGc      <  1:560860/150‑1 (MQ=255)

                                                                                                                                                 cACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCa    >  1:1248230/1‑150 (MQ=255)

                                                                                                                                                  aCTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATTGTGGTCAGCATGCCGTTTTCTGATTTCACCGCTTTCTCGCCATTGTTCATCAGCATAAAGCAGCCATTGCCAt   <  2:478393/150‑1 (MQ=255)

                                                                                                                                                  aCTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAG‑CAGTGCCATa  <  2:208448/150‑1 (MQ=255)

                                                                                                                                                   cTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATa  <  1:138396/150‑1 (MQ=255)

CGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATA  >  NC_000913/4116131‑4116427

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`