BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,494,883	(A)_7→8	coding (261/564 nt)	idnK →	D‑gluconate kinase, thermosensitive

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,494,876	1	.	A	94.7%	400.7 / 12.8	113	coding (254/564 nt)	idnK	D‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand): ref base . (4/2); new base A (58/49); total (62/51)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAG  >  NC_000913/4494777‑4494973

tATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAt‑t                                                                                                      <  1:2832918/100‑1 (MQ=255)

tATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAt‑t                                                                                                      >  1:46834/1‑100 (MQ=255)

tATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAt‑t                                                                                                      >  2:1181306/1‑100 (MQ=255)

 aTCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑a                                                                                                   <  2:288914/100‑1 (MQ=255)

 aTCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑a                                                                                                   <  2:2604719/100‑1 (MQ=255)

     gatgaAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑aaa‑aa                                                                                              >  1:2831814/1‑100 (MQ=255)

         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑aaa‑aaaaac                                                                                          <  1:1701943/100‑3 (MQ=255)

         aaGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑aaa‑aaaaac                                                                                          <  2:1391356/100‑3 (MQ=255)

           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑aaa‑aaaaacag                                                                                        <  1:976629/100‑5 (MQ=255)

           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑aaa‑aaaaacag                                                                                        <  1:3137576/100‑5 (MQ=255)

           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑aaa‑aaaaacag                                                                                        <  1:2941069/100‑5 (MQ=255)

           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑aaa‑aaaaacag                                                                                        <  1:2223697/100‑5 (MQ=255)

           gATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑aaa‑aaaaacag                                                                                        >  2:890395/1‑96 (MQ=255)

              cGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTAt                                                                                      >  2:2356094/1‑100 (MQ=255)

               gACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTAt                                                                                      >  1:1463967/1‑99 (MQ=255)

                 cTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGt                                                                                   >  2:2746547/1‑100 (MQ=255)

                   tCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGa                                                                                 >  1:3071842/1‑100 (MQ=255)

                     ccTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGata                                                                               <  1:1128093/100‑1 (MQ=255)

                      ctggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGatat                                                                              >  2:73943/1‑100 (MQ=255)

                        ggctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATAttt                                                                            <  2:1463964/100‑1 (MQ=255)

                         gctggAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATAtttt                                                                           <  1:858143/100‑1 (MQ=255)

                             gAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGc                                                                       >  1:162580/1‑100 (MQ=255)

                             gAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGc                                                                       <  1:2356101/100‑1 (MQ=255)

                              aaCGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCa                                                                      >  1:2322798/1‑100 (MQ=255)

                              aaCGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCa                                                                      >  1:1943780/1‑100 (MQ=255)

                               aCGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCaa                                                                     >  1:719379/1‑100 (MQ=255)

                               aCGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCaa                                                                     <  1:1187251/100‑1 (MQ=255)

                                cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCGT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAg                                                                    <  2:1834127/100‑1 (MQ=255)

                                cGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAg                                                                    <  2:46834/100‑1 (MQ=255)

                                  cTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAggg                                                                  >  2:3351163/1‑100 (MQ=255)

                                      aaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CATCA‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTa                                                                >  2:3247727/1‑98 (MQ=255)

                                      aaTGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGc                                                              >  2:738890/1‑100 (MQ=255)

                                      aaTGATGATTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGc                                                              >  1:1798572/1‑100 (MQ=255)

                                       atgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGcc                                                             <  1:533196/100‑1 (MQ=255)

                                        tgatgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                            >  1:2762508/1‑100 (MQ=255)

                                        cttccgatcttACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                            <  1:943816/90‑1 (MQ=255)

                                        cttccgatcttACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGccc                                                            <  1:3031193/90‑1 (MQ=255)

                                          atgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGccccc                                                          >  2:239306/1‑100 (MQ=255)

                                          atgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGccccc                                                          >  1:3340369/1‑100 (MQ=255)

                                           tgCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCa                                                         <  2:2930029/100‑1 (MQ=255)

                                            gCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCAt                                                        >  2:1439850/1‑100 (MQ=255)

                                             cTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATg                                                       <  2:2063496/100‑1 (MQ=255)

                                               tCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                     >  2:2735353/1‑100 (MQ=255)

                                                cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGtt                                                     <  2:2371168/99‑1 (MQ=255)

                                                cATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTc                                                    >  2:3211394/1‑100 (MQ=255)

                                                 aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                   <  1:2205449/100‑1 (MQ=255)

                                                 aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                   >  2:1372429/1‑100 (MQ=255)

                                                 aTACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCa                                                   >  1:2276853/1‑100 (MQ=255)

                                                  tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TAAAAA‑AAAACAGTATCTTGTTTTTTTATGCAAGGGTAGCCCCAATGt                                                      >  2:2980748/1‑97 (MQ=255)

                                                  tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGcccagatcggaag                                                  >  2:3031181/1‑90 (MQ=255)

                                                  tACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGcccagatcggaag                                                  >  2:943815/1‑90 (MQ=255)

                                                   aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TAAAAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAt                                                  >  2:2041647/1‑100 (MQ=255)

                                                   aCAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAtt                                                 <  1:1423515/100‑1 (MQ=255)

                                                    cAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCAttt                                                <  2:1805909/100‑1 (MQ=255)

                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                              <  2:3301111/100‑1 (MQ=255)

                                                      gTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTcc                                              <  1:996046/100‑1 (MQ=255)

                                                        cTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctc                                            >  2:2236551/1‑100 (MQ=255)

                                                         tttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCctct                                           <  2:3273692/100‑1 (MQ=255)

                                                          ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACTCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                          >  2:308752/1‑100 (MQ=255)

                                                          ttATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                          <  1:3097691/100‑1 (MQ=255)

                                                          cgatctaaaGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTg                                          <  1:1780167/94‑1 (MQ=255)

                                                           tataAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgg                                         <  2:1333020/100‑1 (MQ=255)

                                                            ataAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCGTGTTCATTTCCTCTGGt                                        >  1:631308/1‑100 (MQ=255)

                                                             taAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGtt                                       <  2:224292/100‑1 (MQ=255)

                                                              aaaaaGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTa                                      <  1:610839/100‑1 (MQ=255)

                                                               aaaaGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAg                                     <  1:3211407/100‑1 (MQ=255)

                                                                aaaGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTgagatcg                                    >  2:1780163/1‑94 (MQ=255)

                                                                aaaGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGa                                    >  2:745704/1‑100 (MQ=255)

                                                                 aaGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                                   <  1:2746557/100‑1 (MQ=255)

                                                                 aaGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                                   <  1:278691/100‑1 (MQ=255)

                                                                 aaGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                                   <  1:1373206/100‑1 (MQ=255)

                                                                    aaTGAAATAGGTTTTATTTTCTGTTCATAAAT‑A‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGAt                                   <  2:54979/97‑1 (MQ=255)

                                                                      tGAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGATAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGa                              >  2:607618/1‑100 (MQ=255)

                                                                       gAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                             >  2:99092/1‑100 (MQ=255)

                                                                       gAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                             >  2:3304448/1‑100 (MQ=255)

                                                                       gAAACAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGAc                             <  1:745704/100‑1 (MQ=255)

                                                                        aaaCAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACt                            >  1:1188899/1‑100 (MQ=255)

                                                                         aaCAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                           <  1:1185732/100‑1 (MQ=255)

                                                                         aaCAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTa                           >  2:2800271/1‑100 (MQ=255)

                                                                         aaCAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑AAAGAAGACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTATATGGCGACTa                           >  2:2234060/1‑100 (MQ=255)

                                                                          aCAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                          <  2:3071830/100‑1 (MQ=255)

                                                                          aCAGGATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTAt                          <  1:659199/100‑1 (MQ=255)

                                                                          aCAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATg                         >  1:23856/1‑100 (MQ=255)

                                                                              gATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaaa                      <  2:1335499/100‑1 (MQ=255)

                                                                              gATTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGaaa                      >  2:2802755/1‑100 (MQ=255)

                                                                              gATTTATTGTCTGTTCAT‑CAT‑T‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                     <  2:23856/100‑1 (MQ=255)

                                                                               ccgatctGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                     <  1:246759/94‑1 (MQ=255)

                                                                               aTTTATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAAc                     <  1:540671/100‑1 (MQ=255)

                                                                                 ttATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTa                   <  1:3184875/100‑1 (MQ=255)

                                                                                  tATTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAt                  >  1:887959/1‑100 (MQ=255)

                                                                                   aTTGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTAtt                 <  2:1188896/100‑1 (MQ=255)

                                                                                     tGTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtct               >  2:1167125/1‑100 (MQ=255)

                                                                                      gTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc              >  1:600912/1‑100 (MQ=255)

                                                                                      gTCTGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATtctc              <  2:2090875/100‑1 (MQ=255)

                                                                                         tGTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcg           >  1:3163264/1‑100 (MQ=255)

                                                                                          gTTCAT‑CAT‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTcgcgc          >  2:1141286/1‑100 (MQ=255)

                                                                                             cat‑cat‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGaa       >  2:1952120/1‑100 (MQ=255)

                                                                                              at‑cat‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAAt      >  1:1933659/1‑100 (MQ=255)

                                                                                               t‑cat‑TA‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATg     >  1:3325304/1‑100 (MQ=255)

                                                                                                     t‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAg  <  1:166998/98‑1 (MQ=255)

TATCTGATGAAGATCGACTTCCCTGGCTGGAACGCTTAAATGATGCTTCATACAGTCTTTATAAAAAGAATGAAACAGGATTTATTGTCTGTTCAT‑CAT‑T‑‑AAA‑AAAACAGTATCGTGATATTTTACGCAAGGGTAGCCCCCATGTTCATTTCCTCTGGTTAGATGGCGACTATGAAACTATTCTCGCGCGAATGCAG  >  NC_000913/4494777‑4494973

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