BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,639,323	T→A	5.1%	H418L (CAT→CTT)	hisS ←	histidyl tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,639,323	0	T	A	5.1%	1937.1 / 16.7	571	H418L (CAT→CTT)	hisS	histidyl tRNA synthetase
Reads supporting (aligned to +/- strand): ref base T (187/337); new base A (42/2); total (230/341)
Fisher's exact test for biased strand distribution p-value = 9.31e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 573 total aligned reads displayed.

ATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAA  >  NC_000913/2639174‑2639471

    gTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTTGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                                                                                                  <  2:9420042/150‑1 (MQ=255)

    gTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                                                                                                  <  2:5574464/150‑1 (MQ=255)

    gTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                                                                                                  <  2:1547373/150‑1 (MQ=255)

    gTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                                                                                                  <  1:4270094/150‑1 (MQ=255)

    gTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                                                                                                  <  2:7883357/150‑1 (MQ=255)

    gTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATgcgc                                                                                                                                                  <  2:6636251/150‑1 (MQ=255)

      gcgcCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCgg                                                                                                                                                <  1:1610009/150‑1 (MQ=255)

      gcgcCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCgg                                                                                                                                                <  2:1915269/150‑1 (MQ=255)

        gcCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTGTTGGTTCTCGGAAATTTGCACGCGGTCCGTCTCCGTCCTTTACCCAGTAAAGTGCGCAAAGGCGCGGCt                                                                                                                                              >  2:8491220/1‑150 (MQ=255)

         ccAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTa                                                                                                                                             <  1:1660979/150‑1 (MQ=255)

          cAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTGTTGGTTCTCGTAAAATTCCACCCGGTCCTTCTCCGTCCGTTACCCAGTAAAGTGCGCAAAATCGCGGCTAc                                                                                                                                            >  1:5079606/1‑150 (MQ=255)

          cAACGCCCAAAATAACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTGTTAGGTCTCGTAAATTTACAAGCAGGCCTTCTCCGTCCTTAACCCATTAAAGAGCGCAAATTcgcgccaac                                                                                                                                            >  2:2237564/1‑145 (MQ=255)

                ccAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTGTTCGGTGTCGTAAAATTACACCCGGTCCGTATCCTTACATAACCCCCTAACGTGCGCAAAATCGCGGCGACCCGATc                                                                                                                                      >  1:7633847/1‑150 (MQ=255)

                 cAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTGGTTGTAGTAAAATTACACCCGGTCCTTATACATAAGTAAACCACTAAAGTGCGCAAAATCGCGGCTACGCAAtac                                                                                                                                     >  1:2088459/1‑148 (MQ=255)

                 cAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATcc                                                                                                                                     <  1:773192/150‑1 (MQ=255)

                   aaaTCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTATCGTAAAATTGCACGCGGTCCGTCTCCGTCCGTAACACAGAAACCTGCGCACATGCGCGGCGACGCTATCCTg                                                                                                                                   >  2:316149/1‑150 (MQ=255)

                    aaTCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTgc                                                                                                                                  <  1:3058322/150‑1 (MQ=255)

                     aTCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTgcg                                                                                                                                 <  2:2041186/150‑1 (MQ=255)

                     aTCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTgcg                                                                                                                                 <  2:9075578/150‑1 (MQ=255)

                        aCCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCaa                                                                                                                              <  1:1838268/150‑1 (MQ=255)

                          cccAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTACATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAAAGTGCGCAAATGCGCGGCTACGCTATCCTGCGCaaat                                                                                                                            >  2:7850711/1‑148 (MQ=255)

                          cccAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACt                                                                                                                            <  1:5164260/150‑1 (MQ=255)

                          cccAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTGTGGGGTGTCGTAAATTTACACGCGGGGCGTATACATACATAAACCAGTCAAGTGCGCAAAATCCCGGCGACGCTATCCAGCCCAAca                                                                                                                            >  1:7297470/1‑149 (MQ=255)

                            cAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCATAACCCAGTAAAGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGc                                                                                                                          >  2:623662/1‑150 (MQ=255)

                            cAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTGTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCGTCCATAACCCAGTAACGGGCGCCAATGCGCGGCTACGCGAGCCTGCGCCACGGc                                                                                                                          >  2:6502518/1‑150 (MQ=255)

                             aaCAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGcc                                                                                                                         <  1:3712232/150‑1 (MQ=255)

                             aaCAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTGTTAGTTCTCGGAAATTTACACCCTGTCCTTCTACGTACATAACCCAGTAACCTGCGCAAAATCGCGGCTAAGCTATACTTCGCAACAGcc                                                                                                                         >  2:141827/1‑150 (MQ=255)

                                  ccAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTg                                                                                                                    <  1:8317397/150‑1 (MQ=255)

                                   cAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTGTCGTAAATTTCCACGCTGTCCTTCTCCGTCCATAAACCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAAATGCCGTTTGc                                                                                                                   >  2:8712049/1‑150 (MQ=255)

                                         tttGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCa                                                                                                             <  2:3676111/150‑1 (MQ=255)

                                         tttGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCa                                                                                                             <  2:8916965/150‑1 (MQ=255)

                                         tttGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTACACCCTGTCCTTCACCTTCCTTAACCCAGTAAAGTGCGCAAATGCGCGGCTACCCAATCCAGCGCAACATCCCTTTGCTCACCa                                                                                                             >  2:5231454/1‑150 (MQ=255)

                                         tttGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTAGTTCTCGTAAATTTCCAACCTGTCCTTCTCCGTCCTTAAACCAGTAACGTGCGCAAATGCGCGGCTACGCTATACAGCGCAACAGCCGGTTGCAAACCa                                                                                                             >  2:4088232/1‑150 (MQ=255)

                                           tGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCaga                                                                                                           <  1:8789158/150‑1 (MQ=255)

                                            gCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCagag                                                                                                          <  1:6587918/150‑1 (MQ=255)

                                            gCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCagag                                                                                                          <  2:852544/150‑1 (MQ=255)

                                             ccATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGGTCTCGTAAAATTCCACCCAGTCCTTCTCCTTCCGTAACCCACTAAACTGCGCAAAAGCGCCGCTACCCGATCCTGCGCAACTGCCGGTTGCTCAccccagc                                                                                                         >  2:420775/1‑145 (MQ=255)

                                                     aGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATcc                                                                                                 <  2:1686271/150‑1 (MQ=255)

                                                       cAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCtt                                                                                               <  2:3782781/150‑1 (MQ=255)

                                                       cAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCtt                                                                                               <  1:862703/150‑1 (MQ=255)

                                                       cAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCtt                                                                                               <  1:8121026/150‑1 (MQ=255)

                                                       cAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCtt                                                                                               <  2:8831289/150‑1 (MQ=255)

                                                       cAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCtt                                                                                               <  2:8281026/150‑1 (MQ=255)

                                                       cAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCtt                                                                                               <  2:3402010/150‑1 (MQ=255)

                                                       cAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTAGTAAATTTACACGCTGTCCTTAACCTTACTTAAACCAGTAAAGTGCGCAAAATCGCGGCTACGCTATCCAGCGCAACAGCCGTTTTCTCACCAGAGCGCAAAAccat                                                                                               >  2:5967413/1‑148 (MQ=255)

                                                        aaaaaaGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTc                                                                                              <  1:1141266/150‑1 (MQ=255)

                                                             aGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCactac                                                                                         <  1:9478034/150‑1 (MQ=255)

                                                              gCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCactact                                                                                        <  1:4039412/150‑1 (MQ=255)

                                                              gCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAAATTCCACGCGGTCCTTCTCCATCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCAACAGAGCGCAAATACATCactaca                                                                                        >  2:699429/1‑149 (MQ=255)

                                                               cGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTActg                                                                                       <  1:6330277/150‑1 (MQ=255)

                                                                 tttAACCGCTTCTACCTGGTCGTGTTGGTTGGCGTAAATTTACAACCGGGCCATCACCGTCCGGAACCCAGGAACGGGCGCCCAAGCGCGGCGACCCGATCCTGCCCCACCGCCGGTTGGTCCCCCGCGCCCCAATCCGTCACAAcagcg                                                                                     >  2:4383547/1‑149 (MQ=255)

                                                                  ttAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTActgctg                                                                                    <  1:8336380/150‑1 (MQ=255)

                                                                   tAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCtgt                                                                                   <  2:9498996/150‑1 (MQ=255)

                                                                   tAACCGCTTCTACCTGGTCGTTTTCGTTATCGTAAAATTCCACCCGGTCCGTATCCATCCTTAAACCACTAAACTGCGCAAAAGCGCGGCGACGCAAACAAGCCCAACAGCCCTTTTCTCACCACACCCCAAAACCTTAACAAAAGCtgt                                                                                   >  2:7743015/1‑150 (MQ=255)

                                                                     aCCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGc                                                                                 <  1:6888131/150‑1 (MQ=255)

                                                                      ccGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGcc                                                                                <  1:3423576/150‑1 (MQ=255)

                                                                      ccGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGcc                                                                                <  2:6760376/150‑1 (MQ=255)

                                                                           tCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGGAACGTGCGCAACGGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAg                                                                           <  1:6259499/150‑1 (MQ=255)

                                                                             tACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGcc                                                                         <  2:5049179/150‑1 (MQ=255)

                                                                              aCCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCa                                                                        <  1:1372764/150‑1 (MQ=255)

                                                                                     cGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTTTCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGa                                                                 <  1:1018634/150‑1 (MQ=255)

                                                                                        tttCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTc                                                              <  1:8491282/150‑1 (MQ=255)

                                                                                        tttCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTc                                                              >  2:3532069/1‑150 (MQ=255)

                                                                                          tCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCAc                                                            >  1:7810738/1‑150 (MQ=255)

                                                                                             ttCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCa                                                         <  1:4637801/150‑1 (MQ=255)

                                                                                              tctcGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAg                                                        <  1:7487853/150‑1 (MQ=255)

                                                                                               ctcGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGc                                                       <  1:3960971/150‑1 (MQ=255)

                                                                                               ctcGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGc                                                       >  1:4781557/1‑150 (MQ=255)

                                                                                               ctcGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGc                                                       <  1:5515105/150‑1 (MQ=255)

                                                                                                tcGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGca                                                      <  1:4000845/150‑1 (MQ=255)

                                                                                                tcGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGca                                                      >  1:1962867/1‑150 (MQ=255)

                                                                                                 cGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGcac                                                     >  2:1006151/1‑150 (MQ=255)

                                                                                                 cGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGGGCGCAAATGCGCGGCTACGCGATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGcac                                                     <  2:5515525/150‑1 (MQ=255)

                                                                                                     aaTTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACcaca                                                 >  1:5739738/1‑150 (MQ=255)

                                                                                                       tttCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGc                                               >  2:2690936/1‑150 (MQ=255)

                                                                                                       tttCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGAGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGc                                               >  1:2766304/1‑150 (MQ=255)

                                                                                                         tCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCaa                                             >  1:9209286/1‑150 (MQ=255)

                                                                                                            aCGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAAcgc                                          <  2:4815530/150‑1 (MQ=255)

                                                                                                                         ccttccttAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCAttt                             <  1:3326202/150‑1 (MQ=255)

                                                                                                                         ccttccttAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCAttt                             >  2:7487946/1‑150 (MQ=255)

                                                                                                                          cttccttAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTa                            >  2:2949562/1‑150 (MQ=255)

                                                                                                                           ttccttAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCACCGCGGGCCCCCCATTTAt                           >  1:6177025/1‑150 (MQ=255)

                                                                                                                             ccttAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCa                         <  1:2237564/150‑1 (MQ=255)

                                                                                                                              cttAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAg                        >  2:7588343/1‑150 (MQ=255)

                                                                                                                              cttAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAg                        >  1:3552111/1‑150 (MQ=255)

                                                                                                                                tAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCa                      <  2:2469200/150‑1 (MQ=255)

                                                                                                                                tAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCa                      <  2:3956515/150‑1 (MQ=255)

                                                                                                                                     cAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGc                 >  2:8032403/1‑150 (MQ=255)

                                                                                                                                     cAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGc                 >  1:2406347/1‑150 (MQ=255)

                                                                                                                                      aGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCa                >  2:4936917/1‑150 (MQ=255)

                                                                                                                                          aCGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACt            <  2:5661321/150‑1 (MQ=255)

                                                                                                                                          aCGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACt            <  1:741431/150‑1 (MQ=255)

                                                                                                                                            gtgtGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGt          <  1:2931459/150‑1 (MQ=255)

                                                                                                                                             tGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGtt         <  2:4601133/150‑1 (MQ=255)

                                                                                                                                             tGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATAAGCACGGGCAAACTGtt         <  2:1140333/150‑1 (MQ=255)

                                                                                                                                              gcgcAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGttt        <  2:3688310/150‑1 (MQ=255)

                                                                                                                                                 cAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCtt     >  2:581302/1‑150 (MQ=255)

                                                                                                                                                    aTGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTaaa  <  1:2518380/150‑1 (MQ=255)

                                                                                                                                                    aTGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTaaa  <  2:7289247/150‑1 (MQ=255)

                                                                                                                                                    aTGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTaaa  <  2:3245264/150‑1 (MQ=255)

ATCAGTGCGCCAACGCCCAAAATCACCCCAACAGCCAGTGCTTTGCCATTTTCAGCAAAAAAGCGTTTAACCGCTTCTACCTGGTCGTTTTCGTTCTCGTAAATTTCCACGCTGTCCTTCTCCTTCCTTAACCCAGTAACGTGCGCAAATGCGCGGCTACGCTATCCTGCGCAACTGCCGTTTGCTCACCAGAGCGCAAATCCTTCACTACTGCTGTGCCGTTAGCCACTTCAGACTCACCCAGCACCACAGCAACGCGGGCACCCCATTTATCAGCACGGGCAAACTGTTTCTTAAA  >  NC_000913/2639174‑2639471

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