BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,109	T→G	94.7%	Q38P (CAA→CCA)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,109	0	T	G	94.7%	330.2 / 4.1	113	Q38P (CAA→CCA)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base T (3/3); new base G (53/53); total (57/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

Only 100 of 113 total aligned reads displayed.

GTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGA  >  NC_000913/4116965‑4117252

gTTTCACGCTGGTTCGTAATAACGATGGCTGCAATTTGAGCGGAGGTGATATCGGCTTGCGCCAGCACTTCTCCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTc                                                                                                                                               <  2:1087656/149‑1 (MQ=255)

gTTTCACGCTGGTTCGTAAGACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTc                                                                                                                                               <  2:95530/149‑1 (MQ=255)

 tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACGTGGTTTTGGGTAGGTTGGCGCa                                                                                                                                              >  2:449976/1‑149 (MQ=255)

 tttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATAGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCCCa                                                                                                                                              >  2:239658/1‑149 (MQ=255)

  ttCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATAGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaa                                                                                                                                             <  2:1475065/149‑1 (MQ=255)

   tCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCaaa                                                                                                                                            <  1:1563336/149‑1 (MQ=255)

       gCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCCAATTcc                                                                                                                                        >  1:105596/1‑148 (MQ=255)

       gCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCCGGCTCTGGGTGGATTGGCCCAAATTcg                                                                                                                                        >  2:512208/1‑149 (MQ=255)

        cTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCTAACCTGGTTTTGGGTAGATTGGCTCAAATTcg‑c                                                                                                                                      <  1:1161587/149‑1 (MQ=255)

        cTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCCCAAATTcg‑c                                                                                                                                      >  2:943898/1‑149 (MQ=255)

         tGGTTCGTAATACCGAGAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTcg‑cg                                                                                                                                     <  2:1768348/149‑1 (MQ=255)

             tCGTAATACCGATAGCTGCAATTTGCTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGGGTTCCATTGGGTCGGGTTCTACCCAACCGCGTTTTGGGTAGTGTCGCTCAAGTGCG‑CGCTGc                                                                                                                                 >  2:209033/1‑149 (MQ=255)

               gTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGATCTACCCAACCTGATTTTGGGTAGATTGGCTCAAATTCG‑CGCGGCGa                                                                                                                               <  1:1166113/149‑1 (MQ=255)

                tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGac                                                                                                                              >  1:1156030/1‑149 (MQ=255)

                 aagaCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGaca                                                                                                                             <  1:846944/146‑1 (MQ=255)

                 aaTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAACTCG‑CGCTGCGaca                                                                                                                             <  2:1101438/149‑1 (MQ=255)

                  aTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGacac                                                                                                                            <  1:962572/149‑1 (MQ=255)

                   tACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACg                                                                                                                           >  1:277276/1‑149 (MQ=255)

                   tACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACg                                                                                                                           >  2:749249/1‑149 (MQ=255)

                    aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGGTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑c                                                                                                                         <  1:708279/149‑1 (MQ=255)

                       gATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTaa                                                                                                                      <  1:1486938/149‑1 (MQ=255)

                        aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGCCACG‑CTAAt                                                                                                                     >  2:1804910/1‑149 (MQ=255)

                        aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAAt                                                                                                                     <  1:1453640/149‑1 (MQ=255)

                        aTAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAAt                                                                                                                     <  2:1218666/149‑1 (MQ=255)

                           gCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGCAGCTTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGat                                                                                                                  >  2:869455/1‑149 (MQ=255)

                                 aTTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGc                                                                                                            <  2:179243/149‑1 (MQ=255)

                                      aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGt                                                                                                       >  2:1331578/1‑149 (MQ=255)

                                      aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCGCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGt                                                                                                       >  2:238101/1‑149 (MQ=255)

                                        cGGAACTGATATCGGCTTTCGCCAGCGCTTCCACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGa                                                                                                     <  2:526185/149‑1 (MQ=255)

                                        cGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGCGTCGTGTTCTACCCAACCTGGGTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTACTGACATTGGCATCCTGa                                                                                                     >  2:532086/1‑149 (MQ=255)

                                          gAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATc                                                                                                   >  1:246914/1‑149 (MQ=255)

                                          gAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATc                                                                                                   >  1:760842/1‑149 (MQ=255)

                                           aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATcc                                                                                                  <  2:227588/149‑1 (MQ=255)

                                               gATACCGGCTTTCGCCAGCACTTCTCCCAGCGTGGGGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTa                                                                                              <  1:1625950/149‑1 (MQ=255)

                                                atatCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTAc                                                                                             >  2:886924/1‑149 (MQ=255)

                                                  atCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGa                                                                                           <  1:1728166/149‑1 (MQ=255)

                                                  atCGGCTTTCGCCAGCACTTCGACCAGCGTGGAGCTTTGGGTGGCCCAGTTTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGa                                                                                           <  1:618726/149‑1 (MQ=255)

                                                   tCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCT‑CGCTGCGACACG‑CTAATGATATTGGCATGGGGATCAATTACGAc                                                                                          >  1:1660089/1‑149 (MQ=255)

                                                   tCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGAc                                                                                          >  1:369461/1‑149 (MQ=255)

                                                    cGGCTTTCGCCAGCACTTCTACCATCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGAcc                                                                                         <  1:1107153/149‑1 (MQ=255)

                                                     ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑TGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACcg                                                                                        <  2:1530395/149‑1 (MQ=255)

                                                     ggCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACcg                                                                                        >  1:1383044/1‑149 (MQ=255)

                                                      gCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGAAGATTGGCTCAATTTCG‑CGCTCCGACACG‑CTAATGATATTGGCATCGTCGTCCATTCCGAcccc                                                                                       >  2:1547651/1‑147 (MQ=255)

                                                       cTTTCGCCAGCACTTCTACCAGCGTGGCGCTTTGGGTTGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACcgca                                                                                      >  1:1795878/1‑148 (MQ=255)

                                                        tttCGCCAGCACTTCTACCAGAGAGGGGCGTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCGTTACGACcgcgc                                                                                     >  2:839961/1‑149 (MQ=255)

                                                        tttCGCCAGCACTTCCACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTGTGGGTAGATTTGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACcgcgc                                                                                     <  1:446097/149‑1 (MQ=255)

                                                          tCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCgg                                                                                   >  1:1576370/1‑149 (MQ=255)

                                                            gCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGa                                                                                 <  2:1083462/149‑1 (MQ=255)

                                                             ccAGCACTTCTACCTGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACGACTGTTG‑TTGTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                                <  2:1661952/149‑1 (MQ=255)

                                                             ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAAGGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                                <  2:1083837/149‑1 (MQ=255)

                                                             ccAGCACCTCTACCAGCGTGGCGCTTTGGGTGGCCCAGTTTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                                <  2:1717447/149‑1 (MQ=255)

                                                              cAGCACTTCTACCAGCGTGGAGCTATGGGTGGCCCAGATTTCCATTGGGTCGTATTGTACCCACCCTGGTTGTGGGTAGATCGGCACAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc                                                                               >  2:816949/1‑149 (MQ=255)

                                                                gCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATGTTGGCATCGTGATCCATTACGACCGCGGGGGGGCtg                                                                             >  1:6251/1‑149 (MQ=255)

                                                                 cACTTCTCCCAGCGGGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATACTGGCATCGTGATCCATTACGACCGCGCGGGAGCtgg                                                                            <  2:879494/149‑1 (MQ=255)

                                                                   cTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCtggtg                                                                          >  2:1528182/1‑149 (MQ=255)

                                                                    ttCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCGATTACGACCGCGCGGGAGCtggggg                                                                         >  1:1101665/1‑149 (MQ=255)

                                                                    tgcTACCAGCGGGGAGCTGTGGGTGGCCCGGCTTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCCAATTCG‑CGCTGCGGCCCG‑CTACTGATATTGGCATCGTGATCGATTACGACCGAGGGGGAGCtggagg                                                                         >  2:166560/3‑149 (MQ=255)

                                                                       tACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACCCG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGc                                                                      >  2:1005143/1‑149 (MQ=255)

                                                                        aCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGcc                                                                     <  1:1445219/149‑1 (MQ=255)

                                                                        aCCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGcc                                                                     <  1:544519/149‑1 (MQ=255)

                                                                         ccAGCGTGGGGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGccc                                                                    <  2:784268/149‑1 (MQ=255)

                                                                          cAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCt                                                                   <  1:1455522/149‑1 (MQ=255)

                                                                            gCGTGGCGCTTTGGGTGGCCCACATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTTCCCTgg                                                                 <  2:630029/149‑1 (MQ=255)

                                                                               tGGGGCTTTGGGGGGCCCAGGTTTCCAGTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCg                                                              <  2:910390/149‑1 (MQ=255)

                                                                               tGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCg                                                              >  1:311850/1‑149 (MQ=255)

                                                                               tGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGCGGTGCCCTGGTCg                                                              >  1:515964/1‑149 (MQ=255)

                                                                               tGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTAGAAATCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCg                                                              >  1:1796672/1‑149 (MQ=255)

                                                                               tGGAGCTTATGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCGGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCg                                                              <  2:567759/149‑1 (MQ=255)

                                                                                  aGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAgc                                                           <  2:281373/149‑1 (MQ=255)

                                                                                         ggTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGa                                                    >  1:832193/1‑149 (MQ=255)

                                                                                             ggccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGatatt                                                >  2:1595974/3‑148 (MQ=255)

                                                                                             gCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGatata                                                >  1:629600/1‑149 (MQ=255)

                                                                                               ccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGCCACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCCACGatctatc                                              >  2:1126806/1‑148 (MQ=255)

                                                                                                   aTTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAtttttt                                          >  1:629448/1‑149 (MQ=255)

                                                                                                    tttCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATAttttttt                                         <  2:1784480/149‑1 (MQ=255)

                                                                                                      tCCATTGGGTGGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTc                                       <  2:348191/149‑1 (MQ=255)

                                                                                                      tCCATTGGGTCGTGTTCTCCCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTc                                       <  1:1349850/149‑1 (MQ=255)

                                                                                                       ccATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCa                                      <  1:1133755/149‑1 (MQ=255)

                                                                                                       ccATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTAGAAACTGC‑TGACT‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCa                                      <  1:1320065/149‑1 (MQ=255)

                                                                                                       ccATTGGGTCGTGTTCTACCCAACCTGGTTTGGGGTAGACTTGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGGGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCa                                      <  1:20679/149‑1 (MQ=255)

                                                                                                         aTTGGGTTGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGt                                    <  2:277276/149‑1 (MQ=255)

                                                                                                          ttGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTc                                   >  1:508598/1‑149 (MQ=255)

                                                                                                             ggTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTccagcacg                                >  1:253455/1‑142 (MQ=255)

                                                                                                             ggTCGTGTTCTACCCAACCTGGTTGTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGAGCTATTTTTTTTCAGCCATg                                >  2:88405/1‑149 (MQ=255)

                                                                                                                  tgtTCTACCCTACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATAGTGTTCCATTACGCCCGCGGGTGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTaa                           <  1:809612/149‑1 (MQ=255)

                                                                                                                  tgtTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAGATTCA‑TGCTGCGACACG‑CTAATGATATTGGCATCGTTATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTaa                           >  2:757372/1‑149 (MQ=255)

                                                                                                                     tCTACCCAACCTGGTTTTGGGTAGCTTGGCTCAAGTTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCTGTCATGGTTAATTg                        >  1:914589/1‑149 (MQ=255)

                                                                                                                      cTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTACGTCATGTTTAAttct                       >  2:697434/1‑147 (MQ=255)

                                                                                                                       tACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTc                      >  2:699494/1‑149 (MQ=255)

                                                                                                                        aCCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTcc                     <  2:1796672/149‑1 (MQ=255)

                                                                                                                         cccAACCTGGTTTTGGGTAGATTTGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCACCGATATATTTTTTTTCAGTCATGTTTAATTGTccc                    >  2:854710/1‑149 (MQ=255)

                                                                                                                          ccAACCTGGTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGGAACGATATATTTTTTTTCAGTCGTGTTTAATTGTCCCg                   >  2:1388976/1‑149 (MQ=255)

                                                                                                                           cAACCTGGGTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGGCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGt                  <  1:1201493/149‑1 (MQ=255)

                                                                                                                               cttgTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGTCATCGTGTTCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTc              <  1:243992/146‑1 (MQ=255)

                                                                                                                               ccggTTTTGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTc              <  1:247414/147‑1 (MQ=255)

                                                                                                                                tGGTTTTGGGTAGATTGGCTCAATTTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCa             >  1:1004480/1‑149 (MQ=255)

                                                                                                                                   ttttGGGTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTTCGTCGTGTTTAATTGTCCCGTTGTCatgt          >  2:885363/1‑147 (MQ=255)

                                                                                                                                         gTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCTAATTatcga    >  2:893991/1‑145 (MQ=255)

                                                                                                                                         gTAGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGGCATATTACAGc    >  1:960789/1‑149 (MQ=255)

                                                                                                                                           aGATTGGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCCTGTTTAATTTGCGCGTAGTCATATTACAgaga  >  2:1053683/1‑149 (MQ=255)

GTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCG‑CGCTGCGACACG‑CTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGA  >  NC_000913/4116965‑4117252

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