BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,713,613	A→T	6.3%	L315Q (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,613	0	A	T	6.3%	341.2 / 10.0	111	L315Q (CTG→CAG)	paaK	phenylacetate‑CoA ligase
Reads supporting (aligned to +/- strand): ref base A (48/56); new base T (4/3); total (52/59)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 115 total aligned reads displayed.

GAATCTCATAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTCC  >  NC_002947/3713464‑3713759

gAATCTCATAAGGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGCTCTGAGTCGGGAACACATTGCCCCCGCGAATGATCAGCATGTTGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGgca                                                                                                                                                     <  2:1004453/150‑1 (MQ=255)

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

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

       aTAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCg                                                                                                                                              >  2:1332843/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:630756/150‑1 (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)

                                                       aGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCTGCTGCCgggca                                                                                                                                          >  1:1228500/1‑104 (MQ=255)

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

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

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

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

                                                           gggAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGc                                                                                          >  1:1353141/1‑150 (MQ=255)

                                                             gAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGa                                                                                        <  2:493169/150‑1 (MQ=255)

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

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

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

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

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

                                                                   aTTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCgg                                                                                                             >  1:324583/1‑123 (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                                                                                                                                             >  1:1291286/1‑84 (MQ=255)

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

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

                                                                                   aTCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCCGAGAGGTG‑AATACCAGCTCGCCCAGt                                                                 <  2:1553300/150‑1 (MQ=255)

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

                                                                                    tCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAgc                                                                                                                       <  2:662977/96‑1 (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                                                                              <  1:1151243/128‑1 (MQ=255)

                                                                                            tcgtcACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATAc                                                                              >  2:1151243/1‑128 (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                                                                              >  2:1199111/1‑123 (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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                    gCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGAt                <  1:487175/150‑1 (MQ=255)

                                                                                                                                      cGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAg              >  2:70437/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                            <  1:1862944/127‑1 (MQ=255)

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

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

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

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

GAATCTCATAAAGCTCAGAAAGCTGTTTTATTTTTAGCACCTGTTCCTCGATCTGAGTCGGGAACACATTGACCCCGCGAATGATCAGCATGTCGTCACTGCGGCCGGTGATCTTGCCGATCCGCCGCATCGGCCGAGCGGTGCCCGGCAGCAGCCGGGTGAGGTCGCGGGTGCGGTAGCGCACCATCGGCAATGCCTCCTTGCTCAGCGAGGTG‑AATACCAGCTCGCCCAGTTGGCCATCCGGCAGCACTTCGCCGGTGGCAGGGTCGATGATTTCGGGATAGAAGTGGTCTTCC  >  NC_002947/3713464‑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: `‑`