BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,116,406	C→T	7.1%	M272I (ATG→ATA)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,116,406	0	C	T	7.1%	324.6 / 8.1	110	M272I (ATG→ATA)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base C (41/61); new base T (2/6); total (43/67)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

TCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTT  >  NC_000913/4116258‑4116554

tcctttcatCGCGCAGCCACTCAATGGAGGCGCCTGCCATAACCACCGTCCCTTCCACCGCATAGTTCACTTCGCCAGTCGGGTCGCAGGCGACGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGc                                                                                                                                                      <  1:302971/146‑1 (MQ=255)

 catttcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCa                                                                                                                                                     >  1:13463/1‑149 (MQ=255)

  atttcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCAt                                                                                                                                                    <  1:969391/149‑1 (MQ=255)

  atctcatCGCGCAGCCACTGAATGGAGGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCAt                                                                                                                                                    <  1:969525/146‑1 (MQ=255)

     tcatCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATaaa                                                                                                                                                 <  1:25976/149‑1 (MQ=255)

     ccgtCGCGCCGCCACTGACTGGTTGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGAAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATaac                                                                                                                                                 <  2:757664/146‑2 (MQ=255)

       cccgcgCAGCCACTGAAGGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGTCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagc                                                                                                                                               <  2:242588/147‑1 (MQ=255)

       atCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAagc                                                                                                                                               >  2:584940/1‑149 (MQ=255)

            gcAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCACGGGTTCATCAGCATAAAGCAGCCa                                                                                                                                          >  2:240941/1‑149 (MQ=255)

                 cACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTCCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGcc                                                                                                                                     <  2:141123/149‑1 (MQ=255)

                  aCTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCa                                                                                                                                    <  1:322378/149‑1 (MQ=255)

                   cTGAATGTATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCAt                                                                                                                                   >  2:462113/1‑149 (MQ=255)

                     gAATGGATGCGCCTGCCTTAAACCCCGCACCTTCCACCGCATAGTTCCCTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAg                                                                                                                                 <  2:1133514/149‑1 (MQ=255)

                     gAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAg                                                                                                                                 <  2:1179070/149‑1 (MQ=255)

                        tGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGgtg                                                                                                                              <  1:1152325/149‑1 (MQ=255)

                           aTGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGttc                                                                                                                           >  1:6387/1‑149 (MQ=255)

                             gcgcCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCATTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGttctt                                                                                                                         <  1:1368464/149‑1 (MQ=255)

                               gcCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCTCTTGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCg                                                                                                                       >  1:773153/1‑149 (MQ=255)

                                 cTGCCATAAACACCGCACCTTCCACCGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGcc                                                                                                                     <  1:175431/149‑1 (MQ=255)

                                     cATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTTAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCGGTGTTCATCAGCATAAAGCAGCCAATGCCATAGGTGTTCTTCGCCATcc                                                                                                                 >  1:276871/1‑149 (MQ=255)

                                     cATAAACACCGAACCTTCCAACGCATAGTTCAATTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGACGTTTTCTGATTTCACCTCTTGCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATcc                                                                                                                 <  1:1146670/149‑1 (MQ=255)

                                       tAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCct                                                                                                               >  1:206742/1‑149 (MQ=255)

                                          aacccctaCCTGCCAACGCAGAGTTCACTTCGCCCGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCCCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCcttct                                                                                                            <  2:363594/142‑1 (MQ=255)

                                           cacCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAACATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCcttctt                                                                                                           <  1:930649/149‑1 (MQ=255)

                                             ccGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTc                                                                                                         <  2:1446815/149‑1 (MQ=255)

                                               gCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCAc                                                                                                       <  1:44319/149‑1 (MQ=255)

                                                cACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCCTAGCGGTCCTTCCCCATCCcttcctccccg                                                                                                      >  2:572378/1‑142 (MQ=255)

                                                cACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGCTGCTCCTCGCCATCCCTTCTTTCACg                                                                                                      <  2:1103302/149‑1 (MQ=255)

                                                 cccgtCCACCGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGc                                                                                                     <  1:4869/145‑1 (MQ=255)

                                                   cTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCCGGGCCATCGGTGTTCTTCGCCATCCCTTCTTCCACGcac                                                                                                   >  2:8179/1‑149 (MQ=255)

                                                      ccAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGCGTTCTTCGCCATCCCTACTTTCACgaccacc                                                                                                >  2:178822/1‑143 (MQ=255)

                                                       cAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTGCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACt                                                                                               <  1:442564/149‑1 (MQ=255)

                                                        aaCGCATAGTTCACGTCGCCAATCGGACCGCAGGCGATGGTGGGCAGCAGGCCGGTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCCCGCACACCTg                                                                                              >  2:855882/1‑149 (MQ=255)

                                                          cGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGAc                                                                                            <  1:1035645/149‑1 (MQ=255)

                                                           gCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCCCAACTGccc                                                                                           >  1:4022/1‑149 (MQ=255)

                                                            cctcgTTCCCTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCa                                                                                          <  2:911294/145‑1 (MQ=255)

                                                            cATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCa                                                                                          <  1:299400/149‑1 (MQ=255)

                                                                gTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCAATATTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAAGTGTTCTTCGCCACCCCTTCTTTCACGCCAACCTGACCAAACa                                                                                      >  2:1096889/1‑149 (MQ=255)

                                                                 ttCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAg                                                                                     >  2:908722/1‑149 (MQ=255)

                                                                    aCTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAgcgc                                                                                  <  2:804141/149‑1 (MQ=255)

                                                                    aCTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAgcgc                                                                                  >  1:1355191/1‑149 (MQ=255)

                                                                    aCTGCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAgcgc                                                                                  <  2:941279/149‑1 (MQ=255)

                                                                      ttCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCgg                                                                                >  1:1417613/1‑149 (MQ=255)

                                                                       tCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGc                                                                               <  2:1424529/149‑1 (MQ=255)

                                                                        cGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACACCTGACCAACCCGCCCgggc                                                                              >  1:624823/1‑147 (MQ=255)

                                                                           cAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCctgc                                                                           <  2:1050439/149‑1 (MQ=255)

                                                                            cgTCGGGCCGCAGGCGATGGGGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCctgct                                                                          <  1:108329/148‑1 (MQ=255)

                                                                             gTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCctgctg                                                                         <  2:206742/149‑1 (MQ=255)

                                                                              tCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTgg                                                                        <  2:437518/149‑1 (MQ=255)

                                                                                gggCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTGTGATGTAACCGCTANCCCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTATTNAACGCACAACTCACCAACCAGCGCGGCCGCCTGGTc                                                                      >  2:813117/1‑149 (MQ=255)

                                                                                 ggCCGCAGGCGATGGTGGTCAGCAGGCCATTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCa                                                                     >  2:1328004/1‑149 (MQ=255)

                                                                                  gCCGCGGGCGGTGGGGGTCAGCAGGGCGTTTTCTGGTTTCACCGCCTTCTCGCCAGTGTTCATCAGCATAAGGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCctgctgaccgt                                                                    <  2:543529/149‑6 (MQ=255)

                                                                                   ccGCAGTCGATGGTGGTCAGCAGCTCGTTTTCTGATTTCCCCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCCCAACTGACCAAACAGCGCGGCCTGCTGGTCAcc                                                                   <  2:1275373/149‑1 (MQ=255)

                                                                                     gCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCgg                                                                 <  2:982980/149‑1 (MQ=255)

                                                                                     gCAGGCGATGGTGGTCAGCAGGCAGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCgg                                                                 >  1:768654/1‑149 (MQ=255)

                                                                                            aTGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATccc                                                          >  1:604548/1‑149 (MQ=255)

                                                                                              ggtggtCACCAGGCCGTTTTCTGTTTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCgg                                                        <  2:635178/149‑1 (MQ=255)

                                                                                               gtggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGga                                                       <  1:1273954/149‑1 (MQ=255)

                                                                                               gtggtCAGCACGCCGTTTTCTGCTTTCACCGCTTTCTCTCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTCCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGga                                                       <  2:732548/149‑1 (MQ=255)

                                                                                                tggtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGgag                                                      >  2:900754/1‑149 (MQ=255)

                                                                                                  gtCAGCGGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGAt                                                    <  2:13463/149‑1 (MQ=255)

                                                                                                  gtCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGAt                                                    >  1:512207/1‑149 (MQ=255)

                                                                                                   tCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGAtt                                                   >  1:463223/1‑149 (MQ=255)

                                                                                                   tCAGCAGGCCGTTTTCTGATTTCACCGCCTTCTCGCCAGTGTTCATCAGTATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGAtt                                                   <  1:282610/149‑1 (MQ=255)

                                                                                                   tCAGCAGGCCGTTTTCTGATTTCACCGCCTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGAtt                                                   <  2:278441/149‑1 (MQ=255)

                                                                                                   ccagaagGCTGTATACTGAGTTCAACGCTTTCTCGCCAGTTCTCACCAGCACAAAGCAGCCATTGCCATATGTGTTCTTCCCCATCCCCTCGTCCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGAtt                                                   <  2:767387/148‑1 (MQ=255)

                                                                                                    cagcagGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGTGATCCCGGAGATTg                                                  <  1:462113/149‑1 (MQ=255)

                                                                                                    cagcagGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTg                                                  <  1:101093/149‑1 (MQ=255)

                                                                                                       cagGCCGTTTTCTCATTTCACCGCTTTCTCCCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCCTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGaa                                               <  1:2703/149‑1 (MQ=255)

                                                                                                        ctcccGTTTTTTTATTGCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCCGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAAt                                              <  1:1322211/146‑1 (MQ=255)

                                                                                                         ggCCGATTTCTGCTTTGACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCTCGGGCTGCTGGTCACCGGCGATCCCGCGGATTGGaaga                                             >  2:174566/1‑147 (MQ=255)

                                                                                                            cGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGCTTGGAATAcgc                                          >  1:1131011/1‑149 (MQ=255)

                                                                                                             gTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCCGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATAcgcg                                         >  2:533614/1‑149 (MQ=255)

                                                                                                               tttCTGAGTTCACCGCTTTCCCGCCAGTGTTCATCCGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGCATAcgcggg                                       <  2:176283/149‑3 (MQ=255)

                                                                                                                ttCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTgc                                      >  1:492357/1‑149 (MQ=255)

                                                                                                                  cTGATTTCACCGCTTTCTCGCCAGTGTTGATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGCTTGGAATACGCGTgccg                                    >  1:874296/1‑149 (MQ=255)

                                                                                                                   tGATTTCACCGCTTTCTCGCCAGTGTTCATCCGCCTCACGCCGCCAGCGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTgccgc                                   <  1:1124097/149‑1 (MQ=255)

                                                                                                                   tGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCATGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGCGATTGGAATACGCGTgccgc                                   >  2:1413761/1‑149 (MQ=255)

                                                                                                                   tGATTTCACCGCTTTCCCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTgccgc                                   <  1:1448584/149‑1 (MQ=255)

                                                                                                                      tttCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCtt                                <  1:908722/149‑1 (MQ=255)

                                                                                                                       ttCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCttt                               <  1:442775/149‑1 (MQ=255)

                                                                                                                       ttCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCttt                               >  1:412189/1‑149 (MQ=255)

                                                                                                                         cACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGTGATTCAAAGACGCGTgacgcgtaggc                             <  2:1041939/149‑7 (MQ=255)

                                                                                                                         cACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGCATACGCGTGCCGCCTCTgc                             >  2:1215838/1‑149 (MQ=255)

                                                                                                                           ccGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccg                           <  1:1291079/149‑1 (MQ=255)

                                                                                                                            cGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGCATACGCGTGCCGCCTTTgccgc                          >  2:718499/1‑149 (MQ=255)

                                                                                                                            cGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccgc                          >  1:349236/1‑149 (MQ=255)

                                                                                                                            cGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTAGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccga                          >  1:824659/1‑148 (MQ=255)

                                                                                                                               tttCTCGCTAGTGGTCTACAGCATAAAGCAGCCCGGGCCCTCAGTGGTCTTCGCCATCCCTTCTGTCACGCACAGCTGACAATACAGTGCGGCGTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCaa                       <  2:1168241/149‑1 (MQ=255)

                                                                                                                               tttCCCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCaa                       <  1:819214/149‑1 (MQ=255)

                                                                                                                                ttCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTgccacccgt                      <  2:903190/149‑4 (MQ=255)

                                                                                                                                    cGCCAGTGTTCATCACCATCAAGCAGCCAGCGCCATAGGTGTTCTTCGCCATCCCTCCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTa                  <  1:564820/149‑1 (MQ=255)

                                                                                                                                    cGCCAGTGGTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACGGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTa                  <  2:1161520/149‑1 (MQ=255)

                                                                                                                                    agctaGTTTTCATCAGTATAAAGCAGCCAGTTCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTa                  <  2:773153/145‑1 (MQ=255)

                                                                                                                                          tgtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCATGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGa            <  2:961768/149‑1 (MQ=255)

                                                                                                                                          tgtTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCAGCCTTTGCCGCCAAGGTTAGTCTGa            >  1:475311/1‑149 (MQ=255)

                                                                                                                                            ttCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAcc          <  2:37996/149‑1 (MQ=255)

                                                                                                                                            ttCATCAGCATAAAGCAGCCAGTGCCATAGGCGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGAcc          <  1:833339/149‑1 (MQ=255)

                                                                                                                                                  aGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATAc    >  2:1032719/1‑149 (MQ=255)

                                                                                                                                                    cATACAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGCAGATTCGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACtt  <  2:1399114/149‑1 (MQ=255)

TCATTTCATCGCGCAGCCACTGAATGGATGCGCCTGCCATAAACACCGCACCTTCCAACGCATAGTTCACTTCGCCAGTCGGGCCGCAGGCGATGGTGGTCAGCAGGCCGTTTTCTGATTTCACCGCTTTCTCGCCAGTGTTCATCAGCATAAAGCAGCCAGTGCCATAGGTGTTCTTCGCCATCCCTTCTTTCACGCACAACTGACCAAACAGCGCGGCCTGCTGGTCACCGGCGATCCCGGAGATTGGAATACGCGTGCCGCCTTTGCCGCCAATGTTAGTCTGACCGTATACTT  >  NC_000913/4116258‑4116554

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