BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	549,796	C→A	6.1%	D156E (GAC→GAA)	tufB →	elongation factor Tu‑B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	549,796	0	C	A	6.1%	557.9 / 17.4	173	D156E (GAC→GAA)	tufB	elongation factor Tu‑B
Reads supporting (aligned to +/- strand): ref base C (86/76); new base A (4/7); total (90/83)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

TCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCTGA  >  NC_002947/549647‑549942

tCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGAc                                                                                                                                                    >  1:4986312/1‑150 (MQ=39)

 cGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGAcc                                                                                                                                                   >  2:2143582/1‑150 (MQ=39)

 cGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGAcc                                                                                                                                                   >  1:1626241/1‑150 (MQ=39)

  ggCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACct                                                                                                                                                  >  1:1666801/1‑150 (MQ=39)

        cgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTc                                                                                                                                            >  1:2884663/1‑150 (MQ=39)

        cgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTc                                                                                                                                            >  1:194483/1‑150 (MQ=39)

         gATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTcc                                                                                                                                           <  1:2452379/150‑1 (MQ=39)

          aTGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCAGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCa                                                                                                                                          <  1:591917/150‑1 (MQ=255)

          aTGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCAGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCa                                                                                                                                          <  1:591708/150‑1 (MQ=255)

           tGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCAc                                                                                                                                         <  2:350656/150‑1 (MQ=39)

            ggcccGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTcgaag                                                                                                                                        >  2:5028980/4‑146 (MQ=255)

              tCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTa                                                                                                                                      <  1:2300129/150‑1 (MQ=39)

                cGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACg                                                                                                                                    >  1:914205/1‑150 (MQ=39)

                  aTGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGAc                                                                                                                                  <  1:4097632/150‑1 (MQ=39)

                       gCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTccc                                                                                                                             >  2:1557656/1‑150 (MQ=38)

                       gCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTccc                                                                                                                             <  2:2006377/150‑1 (MQ=38)

                                tGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGcgacga                                                                                                                    <  1:3175683/150‑1 (MQ=2)

                                tGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGcgacga                                                                                                                    <  1:4744152/150‑1 (MQ=2)

                                   gCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACctgct                                                                                                                                               >  2:78147/1‑120 (MQ=39)

                                   gCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGacac                                                                                                                 >  2:4007078/1‑150 (MQ=37)

                                      caTCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAGCTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGcgacgaagttcc                                                                                                              >  1:1283129/1‑144 (MQ=2)

                                       aTCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCg                                                                                                             >  1:4951601/1‑150 (MQ=37)

                                       aTCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGcgacgaagttcca                                                                                                             >  1:738686/1‑143 (MQ=2)

                                        tCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa                                                                                                            >  2:3089818/1‑150 (MQ=37)

                                        tCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa                                                                                                            <  2:4314122/150‑1 (MQ=37)

                                        tCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa                                                                                                            >  2:1035901/1‑150 (MQ=37)

                                               cgctcGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGCGACGAAGTTCCAatcatcaa                                                                                                     <  2:3112262/146‑2 (MQ=2)

                                                  ccGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCgg                                                                                                  >  2:4554927/1‑150 (MQ=34)

                                                   cGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGt                                                                                                 >  1:4881103/1‑150 (MQ=34)

                                                     tCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTc                                                                                               <  2:1487984/150‑1 (MQ=34)

                                                                 tCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGc                                                                                   >  2:2653622/1‑150 (MQ=34)

                                                                    gTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCt                                                                                <  1:3960075/150‑1 (MQ=37)

                                                                       cATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGa                                                                             <  1:5299957/150‑1 (MQ=37)

                                                                       cATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGa                                                                             <  1:1097917/150‑1 (MQ=37)

                                                                         tCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                                                                           <  1:210207/150‑1 (MQ=37)

                                                                         tCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                                                                           <  2:2142348/150‑1 (MQ=37)

                                                                             ggTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCaa                                                                       <  2:1188997/150‑1 (MQ=37)

                                                                             ggTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCaa                                                                       <  2:1188918/150‑1 (MQ=37)

                                                                                cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                                    <  2:430204/150‑1 (MQ=37)

                                                                                cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                                    <  1:3607060/150‑1 (MQ=37)

                                                                                cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                                    <  1:2562115/150‑1 (MQ=37)

                                                                                cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                                    <  2:3158865/150‑1 (MQ=37)

                                                                                cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                                    >  2:1825209/1‑150 (MQ=37)

                                                                                cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                                                                    <  2:2532607/150‑1 (MQ=37)

                                                                                 ttCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgac                                                                   >  2:4398249/1‑150 (MQ=37)

                                                                                   ccTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacga                                                                 <  2:2359316/150‑1 (MQ=37)

                                                                                   ccTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacga                                                                 <  2:2678890/150‑1 (MQ=37)

                                                                                   ccTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacga                                                                 >  2:5429590/1‑150 (MQ=37)

                                                                                    cTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacgac                                                                >  2:3781273/1‑150 (MQ=37)

                                                                                    cTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacgac                                                                <  2:1337406/150‑1 (MQ=37)

                                                                                     tGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACa                                                               >  1:1626727/1‑150 (MQ=37)

                                                                                     tGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACa                                                               >  2:1467009/1‑150 (MQ=37)

                                                                                       aaCAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAAc                                                             >  1:4101877/1‑150 (MQ=37)

                                                                                         cAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGa                                                           >  1:4601648/1‑150 (MQ=37)

                                                                                          aaGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTAc                                                                                                                                     <  2:2509401/75‑1 (MQ=39)

                                                                                              ccGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                                      <  1:1206205/150‑1 (MQ=35)

                                                                                               cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                                     >  1:1632586/1‑150 (MQ=35)

                                                                                               cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                                     >  2:3855933/1‑150 (MQ=35)

                                                                                               cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                                     <  2:3577833/150‑1 (MQ=35)

                                                                                                 aCCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                                                                           >  2:3421160/1‑126 (MQ=35)

                                                                                                 aCCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                                                                           <  1:3421160/126‑1 (MQ=35)

                                                                                                    tGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGccagg                                                >  1:1276765/1‑146 (MQ=21)

                                                                                                     ggTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTAc                                               >  1:218997/1‑150 (MQ=34)

                                                                                                           cgacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGt                                         <  1:2084506/150‑1 (MQ=34)

                                                                                                           cgacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGt                                         <  2:4962985/150‑1 (MQ=34)

                                                                                                           cgacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGt                                         <  2:3380190/150‑1 (MQ=34)

                                                                                                            gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTTTCCACCTACGACTTCCCAGGCGACGACACTc                                                                                                               <  1:4633362/79‑1 (MQ=21)

                                                                                                            gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatcatc                                                                                                       <  2:3683676/87‑1 (MQ=33)

                                                                                                            gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                                      <  1:2991824/136‑1 (MQ=34)

                                                                                                            gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                                      >  2:2991824/1‑136 (MQ=34)

                                                                                                            gacgCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTc                                        >  1:1389951/1‑150 (MQ=35)

                                                                                                                 tGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCTGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                                   >  2:1323727/1‑150 (MQ=37)

                                                                                                                 tGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                                   >  1:1122967/1‑150 (MQ=37)

                                                                                                                     ctgctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                                      <  2:882663/127‑1 (MQ=35)

                                                                                                                     ctgctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg                                                      >  1:882663/1‑127 (MQ=35)

                                                                                                                      tgctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTgg                              <  2:585500/150‑1 (MQ=37)

                                                                                                                       gctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGt                             >  2:2070401/1‑150 (MQ=37)

                                                                                                                        ctgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg                                                     >  2:2280471/1‑125 (MQ=35)

                                                                                                                         tgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCaagaa                                   >  1:1013607/1‑142 (MQ=37)

                                                                                                                         tgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAg                           >  2:794988/1‑150 (MQ=37)

                                                                                                                         tgGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAg                           >  1:3710944/1‑150 (MQ=37)

                                                                                                                          ggAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGa                          >  1:712814/1‑150 (MQ=37)

                                                                                                                           gAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGaa                         <  1:3686878/150‑1 (MQ=37)

                                                                                                                            aaCTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGaaa                        <  2:2755705/150‑1 (MQ=35)

                                                                                                                                  tCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTgg                  >  1:794577/1‑150 (MQ=34)

                                                                                                                                  tCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTgg                  >  1:4278227/1‑150 (MQ=34)

                                                                                                                                  tCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTgg                  >  1:1472407/1‑150 (MQ=34)

                                                                                                                                   cGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGa                 <  1:2070401/150‑1 (MQ=34)

                                                                                                                                   cGAGATGGAAGTTCGCGACCTGCTGTACACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGa                 >  1:2520787/1‑150 (MQ=21)

                                                                                                                                    gagaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGAt                >  1:2621994/1‑150 (MQ=35)

                                                                                                                                      gaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGc              <  1:4394230/150‑1 (MQ=37)

                                                                                                                                      gaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGc              >  1:4103344/1‑150 (MQ=37)

                                                                                                                                      gaTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGc              <  1:2143582/150‑1 (MQ=37)

                                                                                                                                       aTGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGt                                         >  1:851410/1‑122 (MQ=32)

                                                                                                                                            aGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACAt        <  2:4907852/150‑1 (MQ=255)

                                                                                                                                               tCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATccc     <  2:4986312/150‑1 (MQ=255)

                                                                                                                                                cgcgACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCt    <  2:1934870/150‑1 (MQ=255)

                                                                                                                                                cgcgACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCt    <  2:1666801/150‑1 (MQ=255)

                                                                                                                                                  cgACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCTGa  <  1:2737200/150‑1 (MQ=255)

                                                                                                                                                  cgACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCTGa  <  1:4548422/150‑1 (MQ=255)

TCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGGCACTACCGCTGTCAAGAAGCTGGTAGAAACTCTGGATGCCTACATCCCTGA  >  NC_002947/549647‑549942

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