BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,550,259	A→C	6.4%	N206T (AAC→ACC)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,550,259	0	A	C	6.4%	414.0 / 13.6	126	N206T (AAC→ACC)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand): ref base A (59/58); new base C (5/3); total (64/62)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 134 total aligned reads displayed.

CTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGCCGAAGACGTTG  >  NC_002947/1550110‑1550405

cTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTcaa                                                                                                                                                    <  2:747816/150‑1 (MQ=255)

  ccATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTcaaca                                                                                                                                                  <  1:4681393/150‑1 (MQ=255)

        gTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCgg                                                                                                                                            >  1:463132/1‑150 (MQ=255)

         tGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGa                                                                                                                                           <  2:1483914/150‑1 (MQ=255)

         tGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGa                                                                                                                                           <  1:2649368/150‑1 (MQ=255)

             atcatcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACAcc                                                                                                                                       >  1:4257527/1‑150 (MQ=255)

             atcatcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACAcc                                                                                                                                       <  1:4638657/150‑1 (MQ=255)

             atcatcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACAcc                                                                                                                                       <  1:4638977/150‑1 (MQ=255)

                atcGCTGAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGacaagatg                                                                                                                                    >  2:3704644/1‑145 (MQ=255)

                atcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATg                                                                                                                                    >  2:3248376/1‑150 (MQ=255)

                atcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATg                                                                                                                                    <  1:4957426/150‑1 (MQ=255)

                 tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                                                                                   >  1:2598555/1‑150 (MQ=255)

                 tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                                                                                   >  2:1581516/1‑150 (MQ=255)

                 tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                                                                                   >  1:4838727/1‑150 (MQ=255)

                 tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                                                                                   >  2:4373908/1‑150 (MQ=255)

                 tcGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATgg                                                                                                                                   >  1:1405263/1‑150 (MQ=255)

                   gctgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGtt                                                                                                                                 <  1:708264/147‑1 (MQ=255)

                   gctgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGtt                                                                                                                                 >  2:1640099/4‑150 (MQ=255)

                   gctgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGtt                                                                                                                                 >  1:375046/4‑150 (MQ=255)

                    ccGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTg                                                                                                                                >  2:4870210/1‑150 (MQ=255)

                     tgAAGCGATGCAGAAGGTTGGCAACGAAGGCGTCATCACGGTTGAAGAAGCCAAGACCGCTGAAACCGAACTCGAAGTCGTTGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCATACTTCGTCACCAACCCAGACAAGATGGTTGc                                                                                                                               >  2:15247/2‑150 (MQ=255)

                     cGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGc                                                                                                                               >  1:754393/1‑150 (MQ=255)

                      gAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGcc                                                                                                                              <  2:357938/150‑1 (MQ=255)

                       aaGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCg                                                                                                                             >  1:1799413/1‑150 (MQ=255)

                          ccATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGc                                                                                                                          <  2:2207475/150‑1 (MQ=255)

                          ccATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGc                                                                                                                          >  1:4543638/1‑150 (MQ=255)

                           cATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCt                                                                                                                         >  1:2072622/1‑150 (MQ=255)

                               gAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGaa                                                                                                                     <  2:5061170/150‑1 (MQ=255)

                               gAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGaa                                                                                                                     <  2:3365230/150‑1 (MQ=255)

                               gAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGaa                                                                                                                     >  2:372604/1‑150 (MQ=255)

                                aaaaaGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTCTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGCTGCCGAGCTGGAAg                                                                                                                    <  1:198000/150‑1 (MQ=255)

                                    aGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGccc                                                                                                                <  2:2434435/150‑1 (MQ=255)

                                      tCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgc                                                                                                              <  2:2719932/150‑1 (MQ=255)

                                      tCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgc                                                                                                              <  2:1413316/150‑1 (MQ=255)

                                       cGGTTAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgct                                                                                                             <  1:4461255/150‑1 (MQ=255)

                                            aagaagGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCgctgctgc                                                                                                        <  1:1042920/150‑1 (MQ=255)

                                                 ggCGAGATCCCCGTTGAAGAAGGCTCGGGCCGGGAAAACGAACTGTCCGACGTAGAAGGCAAGCAGTTCGACCGCGGCAACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGCTTGCCGAGCTGGAAGGCCCGctgctgctgctg                                                                                                   >  2:412156/1‑150 (MQ=255)

                                                       aTCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGGACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGAc                                                                                             <  2:411710/150‑1 (MQ=255)

                                                         cACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACaa                                                                                           >  1:4585825/1‑150 (MQ=255)

                                                          aCCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAg                                                                                          >  1:2734522/1‑150 (MQ=255)

                                                          aCCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAg                                                                                          >  1:2734205/1‑150 (MQ=255)

                                                           ccGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGa                                                                                         <  2:2061347/150‑1 (MQ=255)

                                                           ccGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGa                                                                                         >  1:1020613/1‑150 (MQ=255)

                                                             gTTGAAGTAGGCTAGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGTCATGCAGTTCGACCGCGGCTACCTGTCGACGTACTTCGTCATCCAGCCGGACTCCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCAGCTGGTCGACAAGaaa                                                                                       <  2:841994/150‑1 (MQ=255)

                                                             gTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGaaa                                                                                       >  1:509875/1‑150 (MQ=255)

                                                                gaagaaGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAAtc                                                                                    <  1:877663/150‑1 (MQ=255)

                                                                 aagaagGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACACGAAAAtct                                                                                   >  2:1212196/1‑150 (MQ=255)

                                                                       gCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACa                                                                             >  2:135237/1‑150 (MQ=255)

                                                                       gCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACa                                                                             >  2:1517862/1‑150 (MQ=255)

                                                                        cTCGGGTCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACAt                                                                            <  2:670751/150‑1 (MQ=255)

                                                                         tCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATc                                                                           <  2:2169622/150‑1 (MQ=255)

                                                                         tCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTATCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATc                                                                           >  1:2278669/1‑150 (MQ=255)

                                                                          cGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATcc                                                                          <  2:2244103/150‑1 (MQ=255)

                                                                            ggCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGt                                                                        <  1:2593070/150‑1 (MQ=255)

                                                                                acgAAAACCAACTCTCCGTCCTAGAACGCAACCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCCGACACCATGCTAGCCGAGCTGGAAGGCCCGCTGCTGCTGCACGTCGACAACAAAATCTCCAACATCCGTGAgc                                                                    <  2:3490645/148‑1 (MQ=255)

                                                                                    aaaCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAgctgct                                                                <  1:837723/150‑1 (MQ=255)

                                                                                    aaaCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAgctgct                                                                >  2:346314/1‑150 (MQ=255)

                                                                                     aacgaacgGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGctgctg                                                               <  1:3787619/150‑1 (MQ=255)

                                                                                     aaCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGctgctg                                                               >  1:5181232/1‑150 (MQ=255)

                                                                                      aCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAg                                                                                                                           >  1:1376258/1‑89 (MQ=255)

                                                                                      aCGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCtgctgc                                                              <  2:1700604/150‑1 (MQ=255)

                                                                                      aCGAACTGGCAGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGGTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCACCATCCGTGAGCtgctgc                                                              <  1:3730760/150‑1 (MQ=255)

                                                                                          aCTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGctgctgctgctg                                                                                                   >  2:4974126/1‑109 (MQ=255)

                                                                                          aCTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGctgctgctgctg                                                                                                   <  1:4974126/109‑1 (MQ=255)

                                                                                            tGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTc                                                        <  1:4569508/150‑1 (MQ=255)

                                                                                            tGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTc                                                        >  1:817423/1‑150 (MQ=255)

                                                                                             gTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCt                                                       <  1:3978714/150‑1 (MQ=255)

                                                                                                cgtcgtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGAc                                                                                             >  1:2213476/1‑109 (MQ=255)

                                                                                                cgtcgtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGAc                                                                                             <  2:2213476/109‑1 (MQ=255)

                                                                                                cgtcgtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGa                                                    >  1:153331/1‑150 (MQ=255)

                                                                                                  tcgtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAg                                                  >  2:884232/1‑150 (MQ=255)

                                                                                                   cgtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGc                                                 <  2:2956017/150‑1 (MQ=255)

                                                                                                    gtAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGcc                                                <  1:3525968/150‑1 (MQ=255)

                                                                                                     tAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCg                                               >  2:5172331/1‑150 (MQ=255)

                                                                                                      aGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGt                                              >  2:4628122/1‑150 (MQ=255)

                                                                                                       gAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCAGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGtt                                             >  1:4239243/1‑150 (MQ=255)

                                                                                                        aaGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTg                                            <  2:2972965/150‑1 (MQ=255)

                                                                                                         aGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGc                                           >  1:2216835/1‑150 (MQ=255)

                                                                                                           gCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCa                                         >  2:466020/1‑150 (MQ=255)

                                                                                                             aTGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATg                                                                                                                                    >  1:2180961/1‑57 (MQ=255)

                                                                                                             aTGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAg                                       >  1:5202996/1‑150 (MQ=255)

                                                                                                             aTGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAg                                       <  1:4628122/150‑1 (MQ=255)

                                                                                                                  gTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAggccgg                                  <  2:1315697/150‑1 (MQ=255)

                                                                                                                   ttCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAggccggc                                 <  2:2512778/150‑1 (MQ=255)

                                                                                                                   ttCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAggccggc                                 >  2:5431162/1‑150 (MQ=255)

                                                                                                                    tCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAggccggcc                                <  1:466020/150‑1 (MQ=255)

                                                                                                                        ccGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGccc                            >  1:3696362/1‑150 (MQ=255)

                                                                                                                         cgcgGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGcccc                           >  2:824363/1‑149 (MQ=255)

                                                                                                                             gCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCActgc                       <  2:1405263/150‑1 (MQ=255)

                                                                                                                                      cGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTg              <  2:5091704/150‑1 (MQ=255)

                                                                                                                                       gCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGc             <  1:1212196/150‑1 (MQ=255)

                                                                                                                                       gCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGc             <  1:2010047/150‑1 (MQ=255)

                                                                                                                                        ccGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGcc            >  2:401866/1‑150 (MQ=255)

                                                                                                                                        ccGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGcc            >  2:3149089/1‑150 (MQ=255)

                                                                                                                                             cTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAACTCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCCCTGCTGATCGTTGCCGAAGa       >  2:296867/1‑150 (MQ=255)

                                                                                                                                             cTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGCCGAAGa       >  1:274088/1‑150 (MQ=255)

                                                                                                                                                cGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGCCGAAGACGt    >  2:363266/1‑150 (MQ=255)

                                                                                                                                                  tCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGCCGAAGACGTTg  >  2:2159026/1‑150 (MQ=255)

                                                                                                                                                     acaaGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGCCGAAGACGtt   <  1:1909917/146‑1 (MQ=255)

                                                                                                                                                     acaaGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGCCGAAGACGtt   >  2:1909917/1‑146 (MQ=255)

CTCCATCGGTGAAATCATCGCCGAAGCCATGGAAAAAGTCGGTAAAGAAGGCGTGATCACCGTTGAAGAAGGCTCGGGCCTGGAAAACGAACTGTCCGTCGTAGAAGGCATGCAGTTCGACCGCGGCTACCTGTCGCCGTACTTCGTCAACAAGCCGGACACCATGGTTGCCGAGCTGGAAGGCCCGCTGCTGCTGCTGGTCGACAAGAAAATCTCCAACATCCGTGAGCTGCTGCCAGTTCTGGAAGCCGTTGCCAAGGCCGGCCGCCCACTGCTGATCGTTGCCGAAGACGTTG  >  NC_002947/1550110‑1550405

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