BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,245,862	(C)_6→7	coding (1515/1698 nt)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,245,856	1	.	C	96.3%	373.9 / 6.5	108	coding (1521/1698 nt)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand): ref base . (1/3); new base C (45/59); total (46/62)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 122 total aligned reads displayed.

GTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAA  >  NC_000913/1245757‑1245950

gTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCAccgccg                                                                                                   >  2:927489/1‑100 (MQ=255)

gTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCAccgccg                                                                                                   >  1:2266084/1‑100 (MQ=255)

 tCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACcgccg‑‑c                                                                                                <  1:1399813/100‑1 (MQ=255)

  cacaCGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCgccg‑‑c‑c                                                                                              >  2:2222871/1‑100 (MQ=255)

  cacaCGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCgccg‑‑c‑c                                                                                              >  1:2118987/1‑100 (MQ=255)

  cacaCGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCgccg‑‑c‑c                                                                                              <  1:2057515/100‑1 (MQ=255)

    cacGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑‑c‑ccc                                                                                            >  1:721307/1‑100 (MQ=255)

      cGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑‑c‑ccccc                                                                                          <  1:101441/100‑1 (MQ=255)

          tGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑‑c‑ccccccgtt                                                                                      <  2:1990833/100‑5 (MQ=255)

           gCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑‑c‑ccccccgttc                                                                                     >  1:1700614/1‑95 (MQ=255)

            ctctTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTcc                                                                                     >  1:33934/1‑100 (MQ=255)

            ctctTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTcc                                                                                     <  2:2172588/100‑1 (MQ=255)

             tctTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑‑C‑CCCCCGTTCCgg                                                                                   <  1:1332024/100‑1 (MQ=255)

                 tCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCggtgg                                                                                <  1:820621/100‑1 (MQ=255)

                   gggCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGgtggtg                                                                              >  2:2519126/1‑100 (MQ=255)

                      cAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTg                                                                           >  1:19544/1‑100 (MQ=255)

                       agtgcAAAGCCATCGTGTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑‑C‑CCCCCGTTCCGGTGGTGCTg                                                                           <  1:1976989/94‑1 (MQ=255)

                        gATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGAc                                                                         >  1:2784154/1‑100 (MQ=255)

                        gATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGAc                                                                         >  2:1414166/1‑100 (MQ=255)

                        gATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGAc                                                                         <  1:2132567/100‑1 (MQ=255)

                         aTCAAATCCAGCCTTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGAca                                                                        <  2:1989890/100‑1 (MQ=255)

                         aTCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGAca                                                                        <  1:1664119/100‑1 (MQ=255)

                         aTCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGAca                                                                        <  1:1839426/100‑1 (MQ=255)

                         aTCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGAca                                                                        >  2:501904/1‑100 (MQ=255)

                         aTCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGAca                                                                        <  2:1990344/100‑1 (MQ=255)

                          tCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGAcat                                                                       >  2:1453425/1‑100 (MQ=255)

                           cAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGAcatc                                                                      >  2:2559867/1‑100 (MQ=255)

                             aaTCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGAcatcat                                                                    >  1:2686061/1‑100 (MQ=255)

                               tCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑‑C‑CCCCCGTTCCGGTGGTGCTGACATCATAtt                                                                 >  1:1072825/1‑100 (MQ=255)

                                 cAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATAttt                                                                >  1:2208550/1‑100 (MQ=255)

                                  aGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATAtttt                                                               <  2:1082742/100‑1 (MQ=255)

                                   gCATTTTCATCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATAttttt                                                              <  2:509106/100‑1 (MQ=255)

                                   gCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATAttttt                                                              <  2:597726/100‑1 (MQ=255)

                                   gCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATAttttt                                                              >  2:300036/1‑100 (MQ=255)

                                   gCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATAttttt                                                              >  1:303502/1‑100 (MQ=255)

                                    cATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTc                                                             <  2:857868/100‑1 (MQ=255)

                                    cATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTc                                                             <  2:857550/100‑1 (MQ=255)

                                     atttaaaGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTcc                                                            <  1:1147084/94‑1 (MQ=255)

                                     aTTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTcn                                                            <  2:1607900/100‑2 (MQ=255)

                                     aTTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTcc                                                            <  2:934418/100‑1 (MQ=255)

                                     aTTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTcc                                                            >  1:660832/1‑100 (MQ=255)

                                      ttttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTcca                                                           <  1:1408034/100‑1 (MQ=255)

                                      ttttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTcca                                                           <  1:2091915/100‑1 (MQ=255)

                                       tttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTccac                                                          <  1:440089/100‑1 (MQ=255)

                                       tttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTccac                                                          <  2:2827441/100‑1 (MQ=255)

                                       tttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTccac                                                          <  1:880385/100‑1 (MQ=255)

                                        ttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTccacc                                                         <  2:657130/100‑1 (MQ=255)

                                        ttCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTccacc                                                         >  2:1435574/1‑100 (MQ=255)

                                         tCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTccacca                                                        <  1:976250/100‑1 (MQ=255)

                                             cGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCtt                                                    >  2:2182466/1‑100 (MQ=255)

                                              gTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCttt                                                   >  2:1755197/1‑100 (MQ=255)

                                               tCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCtttt                                                  <  2:1872944/100‑1 (MQ=255)

                                                 aCGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTtc                                                <  1:815106/100‑1 (MQ=255)

                                                   gCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTtctc                                              <  1:2884394/100‑1 (MQ=255)

                                                     cATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTccc                                            <  1:1530277/100‑1 (MQ=255)

                                                      aTTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCg                                           <  2:1600954/100‑1 (MQ=255)

                                                       ttGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCgg                                          >  1:941358/1‑100 (MQ=255)

                                                        tGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGt                                         >  2:486191/1‑100 (MQ=255)

                                                        tGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGt                                         <  2:1819432/100‑1 (MQ=255)

                                                         ggTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGt                                         <  2:2588488/99‑1 (MQ=255)

                                                         ggTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTc                                        >  1:885515/1‑100 (MQ=255)

                                                         ggTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTc                                        <  2:2784133/100‑1 (MQ=255)

                                                          gTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCa                                       <  1:1518604/100‑1 (MQ=255)

                                                          gTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCa                                       >  1:2187000/1‑100 (MQ=255)

                                                            ccAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATa                                     >  1:1564854/1‑100 (MQ=255)

                                                            ccAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATa                                     >  2:663556/1‑100 (MQ=255)

                                                                  aaaGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGtgtg                               >  1:1970825/1‑100 (MQ=255)

                                                                  aaaGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGtgtg                               >  1:946300/1‑100 (MQ=255)

                                                                   aaGCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGc                              <  1:300037/100‑1 (MQ=255)

                                                                    ccgatctcTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCt                             <  2:136289/93‑1 (MQ=255)

                                                                     gCCATCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTg                            <  2:1831076/100‑1 (MQ=255)

                                                                        aTCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAAc                         <  2:2192742/100‑1 (MQ=255)

                                                                         tCCTGTAATGGATATTCGCCACCGCCGCCC‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAAc                         <  2:721305/100‑1 (MQ=255)

                                                                         tCCTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACa                        <  1:245959/100‑1 (MQ=255)

                                                                          ncTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACAt                       >  2:161/2‑100 (MQ=255)

                                                                          ccTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACAt                       <  2:2686041/100‑1 (MQ=255)

                                                                           cTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCtagatcgg                      >  1:136289/1‑93 (MQ=255)

                                                                           cTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACAtt                      <  2:2787019/100‑1 (MQ=255)

                                                                           cTGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACAtt                      >  1:912794/1‑100 (MQ=255)

                                                                            tGTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTg                     <  2:788542/100‑1 (MQ=255)

                                                                             tccgatcttatTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTgg                    <  1:519078/92‑1 (MQ=255)

                                                                             gTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTgg                    >  2:1874380/1‑100 (MQ=255)

                                                                             gTAATGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTgg                    >  1:144335/1‑100 (MQ=255)

                                                                               aaTGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTc                  <  1:2699928/100‑1 (MQ=255)

                                                                               aaTGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTc                  <  1:2222884/100‑1 (MQ=255)

                                                                               aaTGGATATTCGCCACCGCCG‑‑CGCCCCCGTTCCGGTGTTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTc                  <  2:1947627/100‑1 (MQ=255)

                                                                                aTGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCa                 >  1:2710104/1‑100 (MQ=255)

                                                                                 tGGTTATTCGCCCCCGCCGC‑C‑CCCCCGTTGCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTc                  <  2:33934/98‑1 (MQ=255)

                                                                                 tGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAg                <  1:2633101/100‑1 (MQ=255)

                                                                                 tGGATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAg                <  1:165835/100‑1 (MQ=255)

                                                                                   gATATTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGa              >  1:1619479/1‑100 (MQ=255)

                                                                                    atatTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGa              <  1:1059809/99‑1 (MQ=255)

                                                                                     tatTCGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTggagatcgga            >  2:519076/1‑92 (MQ=255)

                                                                                         cGCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGcc        >  2:1041990/1‑100 (MQ=255)

                                                                                          gCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGcc        <  2:16244/99‑1 (MQ=255)

                                                                                          gCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCa       >  1:2736861/1‑100 (MQ=255)

                                                                                          gCCACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCa       >  1:1647621/1‑100 (MQ=255)

                                                                                            cACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGc     <  2:2208537/100‑1 (MQ=255)

                                                                                            cACCGCCGC‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCAAGc     >  1:1214987/1‑100 (MQ=255)

                                                                                               cgccgc‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTaa  <  1:1414171/100‑1 (MQ=255)

GTCACACGGTTGCTCTTTCGGGCAGATCAAATCCAGCATTTTCAGCGTCACGCCATTGGTCCAGCCAAAGCCATCCTGTAATGGATATTCGCCACCGCCG‑‑C‑CCCCCGTTCCGGTGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAA  >  NC_000913/1245757‑1245950

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