BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,643,236	A→G	6.6%	R649R (CGT→CGC)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,643,236	0	A	G	6.6%	295.4 / 9.5	107	R649R (CGT→CGC)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand): ref base A (42/56); new base G (2/5); total (46/61)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

Only 100 of 107 total aligned reads displayed.

GGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGC  >  NC_000913/3643088‑3643384

ggCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCCGAACCGCCAc                                                                                                                                                     >  2:907023/1‑150 (MQ=255)

      ggagacGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCa                                                                                                                                               >  1:686518/6‑150 (MQ=255)

       tctcCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAg                                                                                                                                              >  1:1687488/1‑150 (MQ=255)

        ctcCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGa                                                                                                                                             <  1:1137058/150‑1 (MQ=255)

        ctcCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCACCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCAGGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGCACCGCCACGGCTCAGa                                                                                                                                             >  1:230885/1‑150 (MQ=255)

          ccGGTGCCTGTTTCATCAATTAATCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAGCCGTGGCGGTTCAGaag                                                                                                                                           <  2:1653989/150‑2 (MQ=255)

            gtaggcTGTTTCATCAATTAAGCAGATTTTCACTGATTGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAAtgt                                                                                                                                         <  2:124774/145‑1 (MQ=255)

             gTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAGCCGTGGCGGTTCAGaagagc                                                                                                                                        <  1:449/150‑5 (MQ=255)

              tGCCTGTTTCATCTATTAAGCAGATTTTCACTGAATGAGCAGCCCTTAATGCCGTAATGCTCCAGCATCTCATCCAGCTGCGGTTCACGACCACGCAAGCGTTTGAACAGATCCATCGGCTCTTCTGAGCCGTGGCGGTTCAGaagagcc                                                                                                                                       <  2:1506532/150‑6 (MQ=255)

              tGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAGCCGTGGCGGTTCAGaagagcc                                                                                                                                       <  1:1628859/150‑6 (MQ=255)

              tGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAGCCGTGGCGGTTCAGaagagcc                                                                                                                                       <  2:910639/150‑6 (MQ=255)

              tGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGAGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCAACGGCTCAGAAttttt                                                                                                                                       >  2:1497405/1‑149 (MQ=255)

                   gTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAg                                                                                                                                  >  1:525754/1‑150 (MQ=255)

                       catcaATTAAGCAGATTTTCACTGAATGATCCGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCAAGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGacc                                                                                                                              >  2:1258737/1‑148 (MQ=255)

                       catcaATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAAc                                                                                                                              >  1:1171566/1‑150 (MQ=255)

                        atcaATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACg                                                                                                                             <  2:1473951/150‑1 (MQ=255)

                          caATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGGAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCCCGGCTCAGAATGTTGTCGAGGAACGAc                                                                                                                           >  1:234612/1‑150 (MQ=255)

                           aaTTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGCACAGATCCATCGGCTCTTCTGAAACGCCACGGATCAGAATGTTGTCGAGGAACGCCt                                                                                                                          >  2:152447/1‑150 (MQ=255)

                           aaTTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACt                                                                                                                          >  2:907054/1‑150 (MQ=255)

                           aaTTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACt                                                                                                                          <  1:18132/150‑1 (MQ=255)

                           aaTTAACCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACt                                                                                                                          <  1:24156/150‑1 (MQ=255)

                            aTTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGGCGAGGAACGACTg                                                                                                                         <  2:49520/150‑1 (MQ=255)

                              tAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACAACTGcc                                                                                                                       >  2:1381437/1‑150 (MQ=255)

                                aGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCg                                                                                                                     <  1:675718/150‑1 (MQ=255)

                                  cAGATTTTCACTTAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTATGAACCGCCACGTCTCAGAATGTTTTAGAGTCACGACTGCCCGGt                                                                                                                   >  1:1272845/1‑150 (MQ=255)

                                     aTTTTCACTGAATGCGCAGCCCTTAATGCCGTAATGCTCCAGCATCGCCTCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGACACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTc                                                                                                                <  1:471248/150‑1 (MQ=255)

                                     aTTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGGTCACGACCACGGAAGCGTTTGAACAGCTCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTc                                                                                                                <  1:471247/150‑1 (MQ=255)

                                       tttCACTGAATGATCCGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCAc                                                                                                              <  2:880658/150‑1 (MQ=255)

                                          cACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTTAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGt                                                                                                           >  2:1204370/1‑150 (MQ=255)

                                          cACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGt                                                                                                           >  2:642227/1‑150 (MQ=255)

                                           aCTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGtt                                                                                                          >  1:1348875/1‑150 (MQ=255)

                                             tGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGCACGACTGCCCGGTTTCACGGTTGa                                                                                                        >  1:1404803/1‑150 (MQ=255)

                                             tGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGa                                                                                                        <  2:1666372/150‑1 (MQ=255)

                                               aaTGATCAGCCCTAAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTACAAGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGATCAGAATGATGTCGAGGAACGACTGCCCGGTTTCACGGTTGaaa                                                                                                      <  1:876068/150‑1 (MQ=255)

                                                   aTCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGCATGTTGTCGAGGAGCGACTGCCCGGTTAAGCGGTGGAAAATGc                                                                                                  >  2:653150/1‑150 (MQ=255)

                                                    tCAGCCCTTACTGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGcc                                                                                                 <  1:347752/150‑1 (MQ=255)

                                                    tCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGcc                                                                                                 <  2:609885/150‑1 (MQ=255)

                                                    tCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGAGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGcc                                                                                                 <  2:1056168/150‑1 (MQ=255)

                                                     cAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACAACTGCCCGGTTTCACGGTTGAAAATGccc                                                                                                <  1:1381437/150‑1 (MQ=255)

                                                         ccTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGCCTGCCCGGTTTCACGGTTGAAAATGCCCTCtt                                                                                            >  1:1551821/1‑150 (MQ=255)

                                                            tAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCt                                                                                         <  2:529687/150‑1 (MQ=255)

                                                            tAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCt                                                                                         <  2:460882/150‑1 (MQ=255)

                                                                 ccGTAATGCTCCCGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCATAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAg                                                                                    >  2:824813/1‑150 (MQ=255)

                                                                 ccGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCCTCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAATGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAg                                                                                    >  1:1078017/1‑150 (MQ=255)

                                                                 ccGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCATAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAg                                                                                    >  2:484874/1‑150 (MQ=255)

                                                                 ccGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAg                                                                                    <  1:510819/150‑1 (MQ=255)

                                                                   gTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAgcg                                                                                  <  2:686518/150‑1 (MQ=255)

                                                                     aaTGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTGCACGGTTGAAAATGCCCTCTTCCTCAAAGcgcg                                                                                <  2:539391/150‑1 (MQ=255)

                                                                         cTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGaa                                                                            >  2:942043/1‑150 (MQ=255)

                                                                           ccAGCATCGCATCCAGCTGCGGTTCATGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAg                                                                          <  2:200302/150‑1 (MQ=255)

                                                                           ccAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGCACCGCCACGGCTCACAATGAAGTCGAGCAGCGACTGCCCGGATTCTCGGTTGAAAATGCCCACTACCTCAAAGCGCGAGAAAg                                                                          >  1:129903/1‑150 (MQ=255)

                                                                           ccAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAg                                                                          <  2:242815/150‑1 (MQ=255)

                                                                                 tcgcatcCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATTCACTCTTCCTCAAAGCGAGAGAAAGCAACTg                                                                    >  2:268260/1‑150 (MQ=255)

                                                                                 tcgcatcCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATATg                                                                    <  2:1153204/150‑1 (MQ=255)

                                                                                  cgcatcCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAACCGCGAGAAAGCATCTgc                                                                   >  2:83137/1‑150 (MQ=255)

                                                                                    catcCAGCTGCGGTTCACGACCACTGAAGAGTTTTAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGTTTTCACGGTTGAAATTGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTgccg                                                                 <  1:833733/150‑1 (MQ=255)

                                                                                    catcCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTgccg                                                                 <  1:531870/150‑1 (MQ=255)

                                                                                    catcCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTgccg                                                                 <  1:1669656/150‑1 (MQ=255)

                                                                                                 ttCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGc                                                    >  1:772259/1‑150 (MQ=255)

                                                                                                  tCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGcc                                                   <  1:580294/150‑1 (MQ=255)

                                                                                                   cACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCCCGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGccc                                                  >  1:1329597/1‑150 (MQ=255)

                                                                                                      gACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACTGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAACCcaca                                               <  2:1362314/150‑1 (MQ=255)

                                                                                                      gACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCcaca                                               >  1:1553807/1‑150 (MQ=255)

                                                                                                      gACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCcaca                                               <  1:1522943/150‑1 (MQ=255)

                                                                                                        ccACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAgg                                             <  1:1690323/150‑1 (MQ=255)

                                                                                                         cACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGt                                            <  1:441491/150‑1 (MQ=255)

                                                                                                             gAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCt                                        >  1:316226/1‑150 (MQ=255)

                                                                                                              aaGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTg                                       <  2:1496480/150‑1 (MQ=255)

                                                                                                               aGCGTTTGAACTGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTGGGCGAGGACCGACTGCCCGGTTTCACGGTTGAAAATGCCCTCATCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTgt                                      <  2:842131/150‑1 (MQ=255)

                                                                                                               aGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGGTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAACGCATCTGCCGCCACTACGTCAGCCCACAGGTAGCTgt                                      >  1:231806/1‑150 (MQ=255)

                                                                                                                  gTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGCGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCCCAGGTAGCTgtagt                                   >  1:1216710/1‑150 (MQ=255)

                                                                                                                     tGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAAc                                <  1:1297563/150‑1 (MQ=255)

                                                                                                                      gaccAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTTTCTAGTTACGACTGCCAGGTTACACGGTTTAATATGCCCTCTTCCTCAAAGCGGGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGAAGCTGTAGTAAcc                               <  1:1068374/147‑1 (MQ=255)

                                                                                                                      gAACAGATCCATCGGCTCTTCTGAACCGCCCCGGCTCAGAATGTTGTCGGGGCAGGACTGCCGGGTATCCGGGTTGAAAATGCCCACTTCCTCAAAGCGCGAGCAAGCATCTGCCGCCAGTACGACAGCCCCCAGGTAGCAGTATTAccc                               >  2:242801/1‑150 (MQ=255)

                                                                                                                       aaCAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCt                              <  2:1551821/150‑1 (MQ=255)

                                                                                                                       aaCAGATCCAACGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAAccg                              >  1:1443213/1‑149 (MQ=255)

                                                                                                                         cAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGc                            >  2:1527546/1‑150 (MQ=255)

                                                                                                                          aGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGCCTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCg                           >  1:514346/1‑150 (MQ=255)

                                                                                                                            aaccATCGGCTCTTCTGAACCGCCACGGCTCAGAATTTTGTCGACGAACGACTGCCCGGTTTCACGGGTGAAAATGCCCTCTTCCTCAAAGCGCGAGATCGCATCTGCCGCCATTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGc                         <  2:234612/148‑1 (MQ=255)

                                                                                                                             tCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCa                        <  2:906951/150‑1 (MQ=255)

                                                                                                                              ccATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCTCCAGTACGTCAGCCCACAGGTATCTGTTGTAACCTCCGGCAt                       >  2:1052067/1‑150 (MQ=255)

                                                                                                                              ccATCGGCTCTTCCGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCAt                       <  1:907023/150‑1 (MQ=255)

                                                                                                                                aTCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATaa                     <  2:1288589/150‑1 (MQ=255)

                                                                                                                                aTCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATaa                     <  1:679139/150‑1 (MQ=255)

                                                                                                                                 tCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAac                    <  2:1652523/150‑1 (MQ=255)

                                                                                                                                 tCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAac                    <  1:907054/150‑1 (MQ=255)

                                                                                                                                  cGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGATAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAacc                   >  2:352308/1‑150 (MQ=255)

                                                                                                                                  cGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAacc                   <  1:973317/150‑1 (MQ=255)

                                                                                                                                      tcttctGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCg               <  2:1348875/150‑1 (MQ=255)

                                                                                                                                       attctGAGCCGTGGCGGTTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCgg              >  1:1412425/2‑150 (MQ=255)

                                                                                                                                        ttctGAGCCGTGGCGGTTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCCCAGGTAGCTGTACTAACCTGCGGCATAACCCCCGGc             >  1:837760/1‑150 (MQ=255)

                                                                                                                                           tGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCCCAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGaa          >  1:1557860/1‑150 (MQ=255)

                                                                                                                                             aaCCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGaaaa        <  2:1687488/150‑1 (MQ=255)

                                                                                                                                                 gCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCAGAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATg    >  2:213194/1‑150 (MQ=255)

                                                                                                                                                   cACGGCTCAGATTGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGc  <  2:1463989/150‑1 (MQ=255)

                                                                                                                                                   cACGGCTCAGAATGTTGTCGAGGAATGACTGCCCGGTTTCACTGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGc  <  2:1078017/150‑1 (MQ=255)

                                                                                                                                                   cACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGc  <  2:512497/150‑1 (MQ=255)

                                                                                                                                                   cACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGc  <  1:1327736/150‑1 (MQ=255)

                                                                                                                                                   cACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGc  >  2:968297/1‑150 (MQ=255)

                                                                                                                                                   cACGGCTCAGAATGTTGGCGAGGAACGACTGCCCGGTTTCACGGGTGAAAATGCCCTCTTCCTCAAAGCGAGAGAAAGCATCTTCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGc  <  1:152447/150‑1 (MQ=255)

GGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGC  >  NC_000913/3643088‑3643384

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