BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,643,240	C→T	6.6%	S648N (AGC→AAC)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,643,240	0	C	T	6.6%	322.1 / 7.3	106	S648N (AGC→AAC)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand): ref base C (42/55); new base T (2/5); total (45/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 = 6.97e-01

Only 100 of 108 total aligned reads displayed.

GTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAA  >  NC_000913/3643094‑3643389

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

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

  ctcCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGa                                                                                                                                                  <  1:1137058/150‑1 (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)

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

                 catcaATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAAc                                                                                                                                   >  1:1171566/1‑150 (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)

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

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

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

                                                aGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCct                                                                                                    <  2:1315746/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)

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

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

                                                                          atcgcatcCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCt                                                                          <  1:1497405/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)

                                                                                           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)

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

                                                                                                                          aTCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATaa                          <  2:1288589/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)

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

                                                                                                                                                  cTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGaaa  >  2:1483049/1‑150 (MQ=255)

                                                                                                                                                  cTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGaaa  >  1:510003/1‑150 (MQ=255)

GTCTCCGGTGCCTGTTTCATCAATTAAGCAGATTTTCACTGAATGATCAGCCCTTAATGCCGTAATGCTCCAGCATCGCATCCAGCTGCGGTTCACGACCACGGAAGCGTTTGAACAGATCCATCGGCTCTTCTGAACCGCCACGGCTCAGAATGTTGTCGAGGAACGACTGCCCGGTTTCACGGTTGAAAATGCCCTCTTCCTCAAAGCGCGAGAAAGCATCTGCCGCCAGTACGTCAGCCCACAGGTAGCTGTAGTAACCTGCGGCATAACCACCGGCGAAAATATGGCTGAAA  >  NC_000913/3643094‑3643389

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