BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,246,073	C→T	100%	G435D (GGC→GAC)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,246,073	0	C	T	99.1%	359.8 / ‑4.6	105	G435D (GGC→GAC)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (55/49); total (56/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 112 total aligned reads displayed.

GTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGGGCGGCC  >  NC_000913/1245976‑1246168

gTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGgtcg                                                                                               >  2:216286/1‑100 (MQ=255)

   gTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGtcgccg                                                                                            <  2:392168/100‑1 (MQ=255)

   gTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGtcgccg                                                                                            >  2:2191971/1‑100 (MQ=255)

   gTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGtcgccg                                                                                            >  2:1339392/1‑100 (MQ=255)

   gTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGtcgccg                                                                                            <  1:618032/100‑1 (MQ=255)

   gTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGtcgccg                                                                                            >  1:3263715/1‑100 (MQ=255)

    ttttGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCgg                                                                                           <  2:2042778/100‑1 (MQ=255)

    ttttGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCgg                                                                                           <  1:1980327/100‑1 (MQ=255)

     tttGTAATCCTTCTGTCGCGACCCACTGTATCGTTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCggg                                                                                          <  2:1995559/100‑1 (MQ=255)

     tttGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCggg                                                                                          <  2:1151836/100‑1 (MQ=255)

      ttGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGt                                                                                         <  1:1064086/100‑1 (MQ=255)

      ttGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGt                                                                                         >  2:2630695/1‑100 (MQ=255)

        gTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTg                                                                                       <  2:3187320/100‑1 (MQ=255)

          aaTCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgca                                                                                     >  1:111448/1‑100 (MQ=255)

          aaTCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgca                                                                                     <  1:3048597/100‑1 (MQ=255)

           aTCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgcag                                                                                    <  2:2130283/100‑1 (MQ=255)

           aTCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgcag                                                                                    >  2:2646664/1‑100 (MQ=255)

            tCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgcagc                                                                                   <  1:1860471/100‑1 (MQ=255)

            tCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgcagc                                                                                   <  2:1062760/100‑1 (MQ=255)

            tCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgcagc                                                                                   <  2:2802185/100‑1 (MQ=255)

             ccTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgcagca                                                                                  <  1:748429/100‑1 (MQ=255)

              cTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGcagcag                                                                                 >  2:1498234/1‑100 (MQ=255)

              cTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGcagcag                                                                                 >  2:550637/1‑100 (MQ=255)

               ttCTGTCGCGACCCACTGTAACGGTGCCCAGCCA‑TTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAt                                                                               <  2:1595027/100‑1 (MQ=255)

                tCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAt                                                                               <  2:347633/100‑1 (MQ=255)

                 cTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAtg                                                                              <  1:1491133/100‑1 (MQ=255)

                   gTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAtgtg                                                                            >  2:1776047/1‑100 (MQ=255)

                   gTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAtgtg                                                                            >  2:3033945/1‑100 (MQ=255)

                   gTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAtgtg                                                                            >  2:1540932/1‑100 (MQ=255)

                   gTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAtgtg                                                                            >  2:1531488/1‑100 (MQ=255)

                     cgcgACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGtt                                                                          >  2:457971/1‑100 (MQ=255)

                      gcgACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGttt                                                                         <  2:25489/100‑1 (MQ=255)

                      gcgACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGttt                                                                         <  1:646796/100‑1 (MQ=255)

                         aCCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCGATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcg                                                                      >  1:2304709/1‑100 (MQ=255)

                          cccACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgt                                                                     >  2:139263/1‑100 (MQ=255)

                           ccACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtc                                                                    >  1:3245625/1‑100 (MQ=255)

                           ccACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtc                                                                    >  2:1719772/1‑100 (MQ=255)

                            cACTGTACCGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtcg                                                                   <  2:2171286/100‑1 (MQ=255)

                            cACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtcg                                                                   <  1:501889/100‑1 (MQ=255)

                            cACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtcg                                                                   >  2:2961031/1‑100 (MQ=255)

                             aCTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgc                                                                  >  2:802707/1‑100 (MQ=255)

                              cTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgcc                                                                 >  1:562232/1‑100 (MQ=255)

                              cTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgcc                                                                 >  1:1917564/1‑100 (MQ=255)

                              cTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgcc                                                                 >  1:1922057/1‑100 (MQ=255)

                              cTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgc                                                                  >  2:815320/1‑99 (MQ=255)

                               tGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccg                                                                >  1:2347741/1‑100 (MQ=255)

                                 tAACGGTGCCCAGCCATTTGGCGCATCCCATTTTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCTGGGTTGCAGCAGATGTGTTTTcgtcgccgtc                                                              >  1:2454329/1‑100 (MQ=255)

                                 tAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccgtc                                                              <  1:809460/100‑1 (MQ=255)

                                  aaCGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccgtcg                                                             <  2:3263697/100‑1 (MQ=255)

                                    cGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCgtcgccgtcgcc                                                           <  1:460906/100‑1 (MQ=255)

                                      gTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGAAGATGTGTTTTCGTCGCCGTCGCCAt                                                         >  2:2771168/1‑100 (MQ=255)

                                         cccAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAtttt                                                      >  1:1160484/1‑100 (MQ=255)

                                          ccAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttg                                                     >  1:199168/1‑100 (MQ=255)

                                            aGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttgtt                                                   <  2:111447/100‑1 (MQ=255)

                                             gCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttgttg                                                  <  2:1353769/100‑1 (MQ=255)

                                              ccATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTgg                                                 <  2:1468950/100‑1 (MQ=255)

                                                aTTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgcg                                               >  1:1550829/1‑100 (MQ=255)

                                                 tttGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgcgc                                              <  1:2721780/100‑1 (MQ=255)

                                                 tttGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgcgc                                              <  2:2120389/100‑1 (MQ=255)

                                                  ttGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGcgcg                                             <  2:1719281/100‑1 (MQ=255)

                                                  ttGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGcgcg                                             <  2:619463/100‑1 (MQ=255)

                                                     gcgcATCCCATTGTTGCCCACTTTTCGGCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTGGCCGTCGCCATTTTGTTGGCGCGATc                                          <  2:1882142/100‑1 (MQ=255)

                                                     gcgcATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATc                                          >  1:862708/1‑100 (MQ=255)

                                                      cgcATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCt                                         >  2:674162/1‑100 (MQ=255)

                                                      cgcATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCt                                         <  1:708873/100‑1 (MQ=255)

                                                         aTCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTc                                      <  1:1436358/100‑1 (MQ=255)

                                                          tCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCg                                     >  1:2556976/1‑100 (MQ=255)

                                                            ccATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGc                                    <  2:1206853/99‑1 (MQ=255)

                                                            ccATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCt                                   <  1:2256825/100‑1 (MQ=255)

                                                               ttgttgCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGcc                                >  2:187078/1‑100 (MQ=255)

                                                               ttgttgCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGcc                                <  1:179954/100‑1 (MQ=255)

                                                                   tgCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCAt                            >  2:616624/1‑100 (MQ=255)

                                                                    gCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCAtt                           <  2:1844077/100‑1 (MQ=255)

                                                                     cccACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTg                          >  1:500412/1‑100 (MQ=255)

                                                                     cccACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTg                          <  1:2865371/100‑1 (MQ=255)

                                                                      ccACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGa                         >  1:2544460/1‑100 (MQ=255)

                                                                      ccACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGa                         <  1:2646676/100‑1 (MQ=255)

                                                                       cACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGAc                        >  2:837886/1‑100 (MQ=255)

                                                                         cTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGt                      >  1:1135001/1‑100 (MQ=255)

                                                                          ttttCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTa                     >  2:922977/1‑100 (MQ=255)

                                                                          ttttCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTa                     <  1:1164770/100‑1 (MQ=255)

                                                                           tttCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTAc                    >  2:3224496/1‑100 (MQ=255)

                                                                            ttCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGt                      >  2:2443910/1‑97 (MQ=255)

                                                                            ttCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCCTTGACGTAc                    >  2:2230093/1‑99 (MQ=255)

                                                                             tCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAg                  <  1:193249/100‑1 (MQ=255)

                                                                              cACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGc                 >  1:2741359/1‑100 (MQ=255)

                                                                              cACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGc                 >  2:1759571/1‑100 (MQ=255)

                                                                               aCCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCg                <  1:3266759/100‑1 (MQ=255)

                                                                                ccGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCgg               >  1:2530561/1‑100 (MQ=255)

                                                                                 cGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCggg              <  1:496232/100‑1 (MQ=255)

                                                                                  gACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGa             <  2:1850146/100‑1 (MQ=255)

                                                                                   aCGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGaa            >  2:2137463/1‑100 (MQ=255)

                                                                                       ggtgTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTAAAGCgg               >  2:2826156/1‑93 (MQ=255)

                                                                                          gtccAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcg     <  2:381952/97‑1 (MQ=255)

                                                                                          gtTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcg     >  2:1394340/1‑100 (MQ=255)

                                                                                          gtTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcg     >  1:3186136/1‑100 (MQ=255)

                                                                                          gtTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcg     >  1:1127307/1‑100 (MQ=255)

                                                                                            tCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcggc   >  1:1516773/1‑100 (MQ=255)

                                                                                            tCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcgg    <  1:2561439/99‑1 (MQ=255)

                                                                                             cAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGGGCGGcc  >  1:795956/1‑100 (MQ=255)

GTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGGGCGGCC  >  NC_000913/1245976‑1246168

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