BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,713,616	A→T	6.4%	L314Q (CTG→CAG)	paaK ←	phenylacetate‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,713,616	0	A	T	6.4%	338.2 / 10.1	111	L314Q (CTG→CAG)	paaK	phenylacetate‑CoA ligase
Reads supporting (aligned to +/- strand): ref base A (48/55); new base T (4/3); total (52/59)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 114 total aligned reads displayed.

TCTCATAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTCC  >  NC_002947/3713467‑3713759

tctcATAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGTCCGAGCGGTGCCCGgcagca                                                                                                                                                  >  2:1392065/1‑150 (MQ=255)

tctcATAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGgcagca                                                                                                                                                  >  2:711434/1‑150 (MQ=255)

  tcATAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCagcagc                                                                                                                                                <  2:1590682/150‑1 (MQ=255)

    aTAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCg                                                                                                                                              >  2:430470/1‑150 (MQ=255)

     tAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCgg                                                                                                                                             <  2:1042435/150‑1 (MQ=255)

       aaGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGt                                                                                                                                           <  1:282283/150‑1 (MQ=255)

          cTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAg                                                                                                                                        <  1:430470/150‑1 (MQ=255)

               aaaGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCAGGTCGTCACTGCGGCCGGTGATCTTGCCGAGCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCGGCCGGGTGAGGTcgc                                                                                                                                   <  2:165021/150‑1 (MQ=255)

                 aGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCgg                                                                                                                                 <  2:1301868/150‑1 (MQ=255)

                 aGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCgg                                                                                                                                 <  1:1728067/150‑1 (MQ=255)

                          aTTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAg                                                                                                                        <  1:726208/150‑1 (MQ=255)

                          aTTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGAATCGGCAGAGCGGTGCCCTGCAGCAGCCGGGTGAGTTCGCGGGTGCGGTAg                                                                                                                        >  2:208112/1‑150 (MQ=255)

                          aTTTGTTGCACCTGTTCCTGGATATGAGTCGTGAACACATTGACCCCGCGAATGAGATGCATGTCGTTACTGCTGCCGGTGCTCTTGCAGAACAGATGCAGCGGCGGAGGGGGGCCCGGCAGAAGAAGGGGGAGGCCACGGGTGCGGTAg                                                                                                                        >  2:890042/1‑150 (MQ=255)

                           tttttAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAgc                                                                                                                       <  2:888451/150‑1 (MQ=255)

                                 gCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCa                                                                                                                 >  2:373697/1‑150 (MQ=255)

                                 gCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTCCAAACCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGAGCGGCAGAGGGTGCGGTAGCGCACCa                                                                                                                 >  1:1553300/1‑150 (MQ=255)

                                  cACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCAt                                                                                                                >  2:1396113/1‑150 (MQ=255)

                                  cACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCAt                                                                                                                >  2:516951/1‑150 (MQ=255)

                                   aCCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATc                                                                                                               >  1:1326284/1‑150 (MQ=255)

                                     cTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCgg                                                                                                             <  2:763193/150‑1 (MQ=255)

                                      tGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGc                                                                                                            >  2:1471083/1‑150 (MQ=255)

                                         tCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAAt                                                                                                         <  1:789523/150‑1 (MQ=255)

                                         tCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAAt                                                                                                         <  1:630756/150‑1 (MQ=255)

                                           cTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGc                                                                                                       >  1:380420/1‑150 (MQ=255)

                                           cTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGc                                                                                                       >  1:11788/1‑150 (MQ=255)

                                            tCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCaa                                                                                                          >  1:760193/1‑146 (MQ=255)

                                            tCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCaa                                                                                                          <  2:760193/146‑1 (MQ=255)

                                              gATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATCAAGTCACAGGGGCCGGTGATCTTGCCTATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGcctc                                                                                                    >  1:515852/1‑150 (MQ=255)

                                               aTCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGGATGTCGTCACTGCGGCCGGTGTTCTTGCCGATCGGCCGCATCGGCCGCGAGGTGCACGGCAGCAGCCGGGTGAGGTCGCGGGTGTTGTAGCGCACCATCGGCAATGcctcc                                                                                                   <  1:731185/150‑1 (MQ=255)

                                                tCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGcctcct                                                                                                  >  1:1195528/1‑150 (MQ=255)

                                                   gAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGc                                                                                               >  1:362815/1‑150 (MQ=255)

                                                   gAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGc                                                                                               >  1:1644763/1‑150 (MQ=255)

                                                    aGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCTGCTGCCgggca                                                                                                                                          <  2:1228500/106‑3 (MQ=255)

                                                      tCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATc                                                                                                               <  2:1612973/131‑1 (MQ=255)

                                                      tCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATc                                                                                                               >  1:1612973/1‑131 (MQ=255)

                                                      tCGGGAACACATCGCCCCCGCGAATGATCAGCATGTCGTCACTGCCGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGGCCAGAAGCAGCCGAGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCa                                                                                            <  2:130234/150‑1 (MQ=255)

                                                        gggAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCtt                                                                                                 >  2:1422280/1‑143 (MQ=255)

                                                        gggAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCtt                                                                                                 <  1:1422280/143‑1 (MQ=255)

                                                          acacacATTGAGCACGCGAATGTTCAGCATGTCGTCACTGCGTCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGctg                                                                                        >  2:503086/3‑148 (MQ=255)

                                                           aaCACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAg                                                                                       <  2:218373/150‑1 (MQ=255)

                                                           aaCACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAg                                                                                       <  1:1466313/150‑1 (MQ=255)

                                                            acacATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGAATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAgg                                                                                      <  2:533240/150‑1 (MQ=255)

                                                             cacaTTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGt                                                                                     <  1:208112/150‑1 (MQ=255)

                                                              acaTTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTg                                                                                    >  2:1320672/1‑150 (MQ=255)

                                                               caTTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑a                                                                                  >  1:874823/1‑150 (MQ=255)

                                                                aTTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCgg                                                                                                             <  2:324583/123‑1 (MQ=255)

                                                                  tGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATa                                                                               <  1:1332843/150‑1 (MQ=255)

                                                                    aCCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGGCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATAcc                                                                             >  1:1263328/1‑150 (MQ=255)

                                                                     ccccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCa                                                                            <  1:711434/150‑1 (MQ=255)

                                                                     ccccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCa                                                                            <  2:599439/150‑1 (MQ=255)

                                                                      cccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAg                                                                           >  1:1163670/1‑150 (MQ=255)

                                                                      cccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAg                                                                           >  1:994470/1‑150 (MQ=255)

                                                                      cccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAg                                                                           >  1:553174/1‑150 (MQ=255)

                                                                       ccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCgg                                                                                                                                             <  2:1291286/84‑1 (MQ=255)

                                                                       ccGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCgg                                                                                                                                             >  1:1291286/1‑84 (MQ=255)

                                                                            aaTGATCAGCATGTCGTCACTGCGGCCGGTGATCTGGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGcc                                                                     <  1:1585670/150‑1 (MQ=255)

                                                                              tGATCAGCATGGCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCa                                                                   <  1:559936/150‑1 (MQ=255)

                                                                                 tCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGa                                                                                                                                         >  1:1339055/1‑78 (MQ=255)

                                                                                 tCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGa                                                                                                                                         <  2:1339055/78‑1 (MQ=255)

                                                                                 tCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAgc                                                                                                                       <  2:662977/96‑1 (MQ=255)

                                                                                 tCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAgc                                                                                                                       >  1:662977/1‑96 (MQ=255)

                                                                                 tCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATACGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGtt                                                                >  1:365548/1‑150 (MQ=255)

                                                                                   accATTTCTTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTTCCCTGCAGCAGCCGGGTTCTCCTGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTgg                                                              <  1:890042/148‑1 (MQ=255)

                                                                                    gCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGc                                                             <  2:553174/150‑1 (MQ=255)

                                                                                         tcgtcACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATAc                                                                              >  2:1151243/1‑128 (MQ=255)

                                                                                         tcgtcACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATAc                                                                              <  1:1151243/128‑1 (MQ=255)

                                                                                         tcgtcACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATcc                                                        >  2:126366/1‑150 (MQ=255)

                                                                                         tcgtcACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATcc                                                        <  1:1320672/150‑1 (MQ=255)

                                                                                           gtcACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCgg                                                      <  2:380420/150‑1 (MQ=255)

                                                                                            tcACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGGGCGGTAGCGCACCATCGGCAACGAAAACTTGCGCGGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGgc                                                     >  1:1069391/1‑150 (MQ=255)

                                                                                              aCTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATAc                                                                              <  1:1199111/123‑1 (MQ=255)

                                                                                              aCTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGgcag                                                   >  2:1084909/1‑150 (MQ=255)

                                                                                                    gCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTc                                             >  2:1522474/1‑150 (MQ=255)

                                                                                                         tGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCgg                                        <  1:126366/150‑1 (MQ=255)

                                                                                                           aTCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTGGGCCATCCGGCAGCACTTCGCCGGTg                                      <  2:1163670/150‑1 (MQ=255)

                                                                                                              ttGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAAt                                                                                                         <  1:340294/81‑1 (MQ=255)

                                                                                                              ttGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAAt                                                                                                         >  2:340294/1‑81 (MQ=255)

                                                                                                                gCCGATGCGGCGGATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTg                                                               <  2:1639255/120‑1 (MQ=255)

                                                                                                                    aTCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCg                             >  1:1446678/1‑150 (MQ=255)

                                                                                                                     tCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGNCga                            <  2:1195528/150‑1 (MQ=255)

                                                                                                                      ccgccgCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCgat                           <  2:11788/150‑1 (MQ=255)

                                                                                                                      ccgccgCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCgat                           <  1:738227/150‑1 (MQ=255)

                                                                                                                        gccgcATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCa                                                                            <  2:720954/99‑1 (MQ=255)

                                                                                                                        gccgcATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCgatga                         >  2:1146600/1‑150 (MQ=255)

                                                                                                                        gccgcATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCATTGCTCAGCGAGGTG‑AATACCAGAACGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCgatga                         <  1:1471083/150‑1 (MQ=255)

                                                                                                                         ccgcATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCgat                           <  2:911307/147‑1 (MQ=255)

                                                                                                                             aTCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCg                    <  1:1146600/150‑1 (MQ=255)

                                                                                                                             aTCGGCCGAGCGGTGCCCGGCAGCAGCCGAAGCAGACCGCGGGTGCGGTAGCGTACCATGGGCAATGCCTCCTTAC‑CAAAGCGTTGTATTTGCATGTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCg                    >  1:1876455/1‑150 (MQ=255)

                                                                                                                                 gCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGAt                >  1:1666764/1‑150 (MQ=255)

                                                                                                                                    gAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGa             >  2:567073/1‑150 (MQ=255)

                                                                                                                                      gCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAg           <  1:70437/150‑1 (MQ=255)

                                                                                                                                        ggTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATTGCTTCTTGCTCCGCGATGGG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTg         <  2:1401697/150‑1 (MQ=255)

                                                                                                                                        ctcacccGGCTGCTGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTg                                                               >  2:358506/5‑96 (MQ=255)

                                                                                                                                        ctcacccGGCTGCTGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTg                                                               <  1:358506/92‑1 (MQ=255)

                                                                                                                                         tcacccGGCTGCTGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCga                            >  2:1862944/4‑130 (MQ=255)

                                                                                                                                         tcacccGGCTGCTGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCga                            <  1:1862944/127‑1 (MQ=255)

                                                                                                                                          tGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGt       <  2:362815/150‑1 (MQ=255)

                                                                                                                                           acccGGCTGCTGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTc      >  1:1107724/2‑150 (MQ=255)

                                                                                                                                             ccGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCtt    <  2:1263328/150‑1 (MQ=255)

                                                                                                                                               ggCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTcc  <  1:373697/150‑1 (MQ=255)

TCTCATAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTCC  >  NC_002947/3713467‑3713759

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