BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,643,230	A→G	6.7%	G651G (GGT→GGC)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,643,230	0	A	G	6.7%	299.6 / 9.6	105	G651G (GGT→GGC)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand): ref base A (47/51); new base G (2/5); total (49/56)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

Only 100 of 107 total aligned reads displayed.

CAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAA  >  NC_000913/3643081‑3643378

cAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATAGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGaa                                                                                                                                                      >  1:1153204/1‑150 (MQ=255)

 aGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAAc                                                                                                                                                     >  2:1522943/1‑150 (MQ=255)

  gATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAAcc                                                                                                                                                    >  2:1589155/1‑150 (MQ=255)

  gATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTACAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAAcc                                                                                                                                                    <  2:603/150‑1 (MQ=255)

    tAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGc                                                                                                                                                  >  2:1021238/1‑150 (MQ=255)

     aaGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGcc                                                                                                                                                 >  2:1690323/1‑150 (MQ=255)

       ggCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGcccc                                                                                                                                               >  2:1011028/1‑148 (MQ=255)

       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)

                 ccGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGTCTCTTCTGAACCGCCACGGCTCAGAAt                                                                                                                                     <  1:711306/150‑1 (MQ=255)

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

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

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

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

                     tGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAGCCGTGGCGGTTCAGaagagcc                                                                                                                                 <  1:1628859/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                                                                                                                    <  1:18132/150‑1 (MQ=255)

                                  aaTTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACt                                                                                                                    >  2:907054/1‑150 (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)

                                                      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:460882/150‑1 (MQ=255)

                                                                   tAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCt                                                                                   <  2:529687/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:1669656/150‑1 (MQ=255)

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

                                                                                                        ttCACGACCACGGAAGCGGTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGc                                              <  1:269086/150‑1 (MQ=255)

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

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

                                                                                                            cGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCcac                                          >  2:1667505/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)

                                                                                                                     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               <  1:679139/150‑1 (MQ=255)

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

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

                                                                                                                                        tCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAac              <  2:1652523/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)

CAGATAAGGCACCGTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAA  >  NC_000913/3643081‑3643378

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