BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,566,212	C→G	9.4%	G124R (GGT→CGT) ‡	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,566,212	0	C	G	9.4%	320.2 / 19.5	105	G124R (GGT→CGT) ‡	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base C (46/49); new base G (5/5); total (51/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 106 total aligned reads displayed.

AGCCAGTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGG  >  NC_000913/3566064‑3566360

agccagTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCg                                                                                                                                                     <  1:449597/150‑1 (MQ=255)

   cagTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                                                                                                                                  >  2:1658682/1‑150 (MQ=255)

   cagTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                                                                                                                                  >  2:408981/1‑150 (MQ=255)

     gtagtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTAcc                                                                                                                                                <  2:1669520/150‑1 (MQ=255)

       agtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCGAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGa                                                                                                                                              >  2:1198360/1‑150 (MQ=255)

       agtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGATTCAAGGAAGCAAGCCGCCAGACGTCCCCGGCCACCGTTACCgc                                                                                                                                              >  1:102511/1‑149 (MQ=255)

        gtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAg                                                                                                                                             >  2:1325676/1‑150 (MQ=255)

        gtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAg                                                                                                                                             <  1:614001/150‑1 (MQ=255)

         tCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAgg                                                                                                                                            <  1:13902/150‑1 (MQ=255)

            ggCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCt                                                                                                                                         >  1:1636943/1‑150 (MQ=255)

              cGACTCTATCTGGCTACCGTTAACGATGATCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTgg                                                                                                                                       <  2:856076/150‑1 (MQ=255)

                aCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATAGCGTATTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCTAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTggga                                                                                                                                     <  1:1279048/150‑2 (MQ=255)

                aCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTggga                                                                                                                                     <  1:1223626/150‑2 (MQ=255)

                aCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTggga                                                                                                                                     <  1:1709580/150‑2 (MQ=255)

                aCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTggga                                                                                                                                     <  1:1703084/150‑2 (MQ=255)

                aCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGt                                                                                                                                     <  2:528781/150‑1 (MQ=255)

                 ctctTTCTGGCTACCGTTTACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGCGGCCTGGGTc                                                                                                                                    >  2:560515/1‑150 (MQ=255)

                   ctTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCAt                                                                                                                                  >  1:653089/1‑150 (MQ=255)

                        tGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttct                                                                                                                             <  1:471182/150‑1 (MQ=255)

                        tGGCTACCGTTAACGATGGTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttct                                                                                                                             <  1:698891/150‑1 (MQ=255)

                          gCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCCCCGTTACCGAGGCCTGGGTCATTttcttc                                                                                                                           >  2:1107968/1‑150 (MQ=255)

                          gCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttcttc                                                                                                                           <  2:1594404/150‑1 (MQ=255)

                          gCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttcttc                                                                                                                           <  1:548793/150‑1 (MQ=255)

                           cTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTg                                                                                                                                        <  1:1324496/136‑1 (MQ=255)

                           cTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTg                                                                                                                                        >  2:1324496/1‑136 (MQ=255)

                           cTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtca                                                                                                                          >  2:546367/1‑150 (MQ=255)

                             aCCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtcatc                                                                                                                        <  2:1273036/150‑1 (MQ=255)

                              ccGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCCCCGTTACCGAGGGCTGGGTCATTTTAGtcatca                                                                                                                       >  2:1523504/1‑150 (MQ=255)

                              ccGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtcatca                                                                                                                       >  1:126056/1‑150 (MQ=255)

                                gTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTAACGAGGCCTGGGTCATTttcttaaaaaat                                                                                                                     >  2:1499144/1‑143 (MQ=255)

                                 ttAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCatcaatc                                                                                                                    <  1:888766/150‑1 (MQ=255)

                                  tAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCatcaatca                                                                                                                   <  1:546367/150‑1 (MQ=255)

                                    aCGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGc                                                                                                                 <  2:213204/150‑1 (MQ=255)

                                        tGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCtt                                                                                                             <  1:1325676/150‑1 (MQ=255)

                                          ttCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCg                                                                                                           <  2:935920/150‑1 (MQ=255)

                                            cTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGTCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCGGGGTCATTTTCTTCGTCGATCAGCTCTTCgag                                                                                                         >  1:664106/1‑150 (MQ=255)

                                             tGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCgaga                                                                                                        <  1:105542/150‑1 (MQ=255)

                                                 tGAACATACCGTAGTCATAGCGGAGGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTa                                                                                                    <  1:702440/150‑1 (MQ=255)

                                                   aaCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGAAATCCGCCAGACGTCCGAGGCCTCGGGTACCGAGGGCTGGGTGATTTTATTCATAAATCAGGTCTTCGAGATTTTAc                                                                                                  >  2:1372042/1‑150 (MQ=255)

                                                    aCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAcc                                                                                                 <  2:1566810/150‑1 (MQ=255)

                                                    aCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAcc                                                                                                 <  1:1039068/150‑1 (MQ=255)

                                                     cATACCGTTGTCATTGCGGATGCCGTATCCGCGCCGCGGCATCCCTAACGTCGCAAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAccc                                                                                                <  1:843290/150‑1 (MQ=255)

                                                     cATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAccc                                                                                                >  2:567831/1‑150 (MQ=255)

                                                       tACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGCGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCa                                                                                              >  2:1059306/1‑150 (MQ=255)

                                                            tAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCt                                                                                         <  1:650175/150‑1 (MQ=255)

                                                               tCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTcc                                                                                      <  1:952350/150‑1 (MQ=255)

                                                                cATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCa                                                                                     <  1:1658682/150‑1 (MQ=255)

                                                                cATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTATTCATCAATCAGCTCTTCGAGTTTTAACCCCATCGCTTCCa                                                                                     >  2:569220/1‑150 (MQ=255)

                                                                  tAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCCGGCCCCCGGTGCCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGAGTTAACCCCCTCGCTTCCAgg                                                                                   >  1:599783/1‑149 (MQ=255)

                                                                   aGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTg                                                                                  <  1:138459/150‑1 (MQ=255)

                                                                   aGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGGCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTg                                                                                  <  2:126056/150‑1 (MQ=255)

                                                                   aGCGGATGACGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTg                                                                                  <  1:1499144/150‑1 (MQ=255)

                                                                     cggcggCCGTCACCGCGCCCCGGCCACCCTAACGTCTCCAGAGAATCAACGAAGCACGCCGCCAGACGTCCCACGCCACCGTTACCGTGGCCTTGGTCATTTTCTTCCTCACTCAGCTCATCGAGATTTAACCCCATCGCTTCCAGTgcg                                                                                <  2:121188/145‑1 (MQ=255)

                                                                     cGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTgcg                                                                                <  2:592026/150‑1 (MQ=255)

                                                                                 ccGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCCCCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCt                                                                    >  1:659275/1‑150 (MQ=255)

                                                                                  cgcgCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCGCTTCGAGAGTTAACCCCATCGCTTCCAGGGCGCCCTGGACATcgt                                                                   >  1:667805/1‑148 (MQ=255)

                                                                                     gcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGt                                                                >  2:428745/1‑150 (MQ=255)

                                                                                        ccGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTaaa                                                             >  2:202342/1‑150 (MQ=255)

                                                                                         cGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCACGCCACCCTTACCGCGGCCTGGGTCAATTAATACATCAATCAGCACTTCGAAAATTAACCCCATAGCGTCCACAGCCCCCTGTACATCTTCGTaaaa                                                            >  2:1081158/1‑149 (MQ=255)

                                                                                            cAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACTTCCCAGGCCACCGTTACCGAGGCCTGGGGCAGTTTCGTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCATCCAGTGCGCCCTGTACAACTTCGTAAATTcg                                                         >  2:1047661/1‑149 (MQ=255)

                                                                                               cccTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAg                                                      <  1:849016/150‑1 (MQ=255)

                                                                                               cccTAACGTCGCCAGAGAATCAAGGAAGAAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTACATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAg                                                      >  2:1111186/1‑150 (MQ=255)

                                                                                                   aaCGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGac                                                  <  1:1623750/150‑1 (MQ=255)

                                                                                                   aaCGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGac                                                  >  1:1218032/1‑150 (MQ=255)

                                                                                                     cGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTTCCGAGGCCTGGGTCATTTTCTTCATAAATCAGCTCGTCGAGATTTAACCCCATCGCTTACAGTGCGCCCTGTACCTCTTCGTAAATTCCTAGCGacca                                                >  2:591480/1‑148 (MQ=255)

                                                                                                      gTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGacaac                                               <  2:450302/150‑1 (MQ=255)

                                                                                                      gTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGacaac                                               <  1:312417/150‑1 (MQ=255)

                                                                                                         gCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATg                                            <  2:653089/150‑1 (MQ=255)

                                                                                                          ccAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATgg                                           <  1:1143474/150‑1 (MQ=255)

                                                                                                          ccAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCCTTACCGCGGCCTGGGGCATTTTCGTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTTCGCCCGGTACAACGTCGTAAATTCCTAGCGACAACAAgg                                           >  2:1590780/1‑150 (MQ=255)

                                                                                                           cAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGATGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGAGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGc                                          <  2:1335343/150‑1 (MQ=255)

                                                                                                             gagaATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGt                                        <  1:1059306/150‑1 (MQ=255)

                                                                                                                aaTCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTgg                                     >  1:188932/1‑150 (MQ=255)

                                                                                                                   cAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCCTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGgaga                                  >  2:1658828/1‑150 (MQ=255)

                                                                                                                     aGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGc                                <  1:408981/150‑1 (MQ=255)

                                                                                                                     aGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGc                                <  2:188932/150‑1 (MQ=255)

                                                                                                                      ggAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGTGAGCg                               >  1:62462/1‑150 (MQ=255)

                                                                                                                       gAAGCAAGCAGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCATCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGt                              <  1:258281/150‑1 (MQ=255)

                                                                                                                          gcaagcCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTATTCAAAAATCACCTCTTCGAGATTTAACCCCAACGCTTCCAGTGCGCCCAGTACATCTTCGTAAATTACTAGCGACAACATCGCGTTGGAGAGCGTACg                           >  2:1686778/1‑150 (MQ=255)

                                                                                                                            aagcCGCCAGACGTCCCAGGCCACCGTTAACGACGCCTGGGTCATTTTGTTCAACAATAAGCTATTAGAGATTTTAACCCATAGCTTACACTGAGCCCGGTTAAAATTAGTAAAATTATAGCGACAACAAGGAGGTTGGGAGCGTaaagc                         >  2:1422052/1‑146 (MQ=255)

                                                                                                                              gccgccAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGAGTTCCCCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAAGTCCTAGCGACAACATGGCGGTGGCGAGCGTACGGCCa                       >  1:17985/1‑150 (MQ=255)

                                                                                                                                cgccAGTCGTCCCAGGCCACCGTTACCGAGTCCGGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATAGCTTCCAGTGCGCCCTGTACATGTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAAt                     <  1:1081158/150‑1 (MQ=255)

                                                                                                                                cgccAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCATCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAAt                     >  2:728484/1‑150 (MQ=255)

                                                                                                                                   aaaatgacccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaa                  >  2:783651/8‑150 (MQ=255)

                                                                                                                                     gACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaaaa                >  2:1293841/1‑150 (MQ=255)

                                                                                                                                     gACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaaaa                <  1:1210660/150‑1 (MQ=255)

                                                                                                                                       cGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTTGAGAGCGTACGGCCAATCAAAAAc              <  1:1686778/150‑1 (MQ=255)

                                                                                                                                        gTCCCAGGCCACCGTTACCGAGGCCTGGGACATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACt             >  2:525677/1‑150 (MQ=255)

                                                                                                                                         acccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTc            >  2:375299/2‑150 (MQ=255)

                                                                                                                                          cccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTcc           >  2:1417567/1‑150 (MQ=255)

                                                                                                                                           ccAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCa          <  1:567831/150‑1 (MQ=255)

                                                                                                                                           ccAGGCCACCGTTACCGAGGCCTGGGTCATTTTCATCAACAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAAATCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCa          <  1:1372042/150‑1 (MQ=255)

                                                                                                                                            cAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCAt         >  1:1088976/1‑150 (MQ=255)

                                                                                                                                            cAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCCTCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGAATTAGGAGAGCGTACGGCCAATCAAAACCTCCag         >  2:282727/1‑149 (MQ=255)

                                                                                                                                            cAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCAt         >  2:1209644/1‑150 (MQ=255)

                                                                                                                                            cAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCAt         >  1:904067/1‑150 (MQ=255)

                                                                                                                                             cggCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCGCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATc        <  2:1303303/149‑1 (MQ=255)

                                                                                                                                               gCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGa      >  2:493680/1‑150 (MQ=255)

                                                                                                                                                 cACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACa    >  1:45877/1‑150 (MQ=255)

                                                                                                                                                   ccGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAgg  >  1:608953/1‑150 (MQ=255)

AGCCAGTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGG  >  NC_000913/3566064‑3566360

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