BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,246,015	T→C	A454A (GCA→GCG)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,246,015	0	T	C	99.1%	391.1 / NA	106	A454A (GCA→GCG)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (68/37); total (68/38)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

TGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGG  >  NC_000913/1245870‑1246162

tggtgCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGccc                                                                                                                                                 >  1:1208393/1‑150 (MQ=255)

 ggtgCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCa                                                                                                                                                <  1:1470542/150‑1 (MQ=255)

      tGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCAtt                                                                                                                                          >  1:31082/1‑151 (MQ=255)

       gACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCAttt                                                                                                                                         <  2:37608/151‑1 (MQ=255)

       gACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCCGTCGCGACCCACTGTAACGGCGCCCAGCCAttt                                                                                                                                         <  2:1212625/151‑1 (MQ=255)

        aCATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTg                                                                                                                                        <  1:617636/151‑1 (MQ=255)

         catcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTgg                                                                                                                                       <  1:1617057/151‑1 (MQ=255)

          atcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGgc                                                                                                                                      <  2:1064837/151‑1 (MQ=255)

          atcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGgc                                                                                                                                      >  1:1252342/1‑151 (MQ=255)

           tcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGgcg                                                                                                                                     >  1:1331920/1‑151 (MQ=255)

           tcatATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGgc                                                                                                                                      >  2:901029/1‑150 (MQ=255)

            catATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGgcgc                                                                                                                                    >  2:672586/1‑151 (MQ=255)

                 ttttCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATcc                                                                                                                                <  1:1252467/150‑1 (MQ=255)

                  tttCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATccc                                                                                                                               >  1:1296726/1‑150 (MQ=255)

                     ccaccaGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCAttg                                                                                                                           >  1:1143833/1‑151 (MQ=255)

                      caccaGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCAttgt                                                                                                                          <  2:1299329/151‑1 (MQ=255)

                        ccaGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCAttgttg                                                                                                                        <  2:1072597/151‑1 (MQ=255)

                         caGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGc                                                                                                                       <  1:519297/151‑1 (MQ=255)

                            cTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGccc                                                                                                                     >  1:403081/1‑150 (MQ=255)

                             tttttCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCAc                                                                                                                   >  1:1443679/1‑151 (MQ=255)

                             tttttCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCAc                                                                                                                   >  2:1047280/1‑151 (MQ=255)

                              ttttCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACt                                                                                                                  >  1:1336135/1‑151 (MQ=255)

                                  ctcCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTc                                                                                                              >  1:227415/1‑151 (MQ=255)

                                    cccGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCAc                                                                                                            <  1:509062/151‑1 (MQ=255)

                                    cccGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCAc                                                                                                            <  1:1581677/151‑1 (MQ=255)

                                         tCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACg                                                                                                       >  1:1239029/1‑151 (MQ=255)

                                           aTAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtg                                                                                                     >  1:334351/1‑151 (MQ=255)

                                            tAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtgg                                                                                                    >  1:309578/1‑151 (MQ=255)

                                            tAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtgg                                                                                                    >  1:1281723/1‑151 (MQ=255)

                                             aGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtgg                                                                                                    <  1:16801/150‑1 (MQ=255)

                                              ggTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACgtggt                                                                                                   <  1:194265/150‑1 (MQ=255)

                                                tgtgCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGtt                                                                                                >  2:1670540/1‑151 (MQ=255)

                                                     tGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGgc                                                                                           >  1:1327013/1‑151 (MQ=255)

                                                      gAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGgcc                                                                                          <  2:1030836/151‑1 (MQ=255)

                                                          aTTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGccgccg                                                                                      <  2:573809/151‑1 (MQ=255)

                                                            tGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCggg                                                                                    >  1:1696722/1‑151 (MQ=255)

                                                            tGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCggg                                                                                    >  1:205765/1‑151 (MQ=255)

                                                             ggTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGt                                                                                   <  1:1615499/151‑1 (MQ=255)

                                                              gTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGtt                                                                                  <  2:1613452/151‑1 (MQ=255)

                                                              gTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGtt                                                                                  <  1:717068/151‑1 (MQ=255)

                                                                cAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgc                                                                                <  1:1321973/151‑1 (MQ=255)

                                                                 aGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgca                                                                               <  2:601418/151‑1 (MQ=255)

                                                                  ggAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgc                                                                                >  2:324090/1‑149 (MQ=255)

                                                                   gAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgcagc                                                                             >  2:112604/1‑151 (MQ=255)

                                                                    aaGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgcagca                                                                            <  2:197257/151‑1 (MQ=255)

                                                                    aaGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgcagc                                                                             >  1:1382640/1‑150 (MQ=255)

                                                                     aGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGcagcag                                                                           <  1:1090121/151‑1 (MQ=255)

                                                                     aGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGcagcag                                                                           <  2:1296820/151‑1 (MQ=255)

                                                                       tGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgcagc                                                                             >  1:438539/1‑147 (MQ=255)

                                                                        gCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGAtg                                                                        >  2:1285953/1‑151 (MQ=255)

                                                                        gCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGAtg                                                                        <  2:1533668/151‑1 (MQ=255)

                                                                         ccAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGAtgt                                                                       <  1:1259679/151‑1 (MQ=255)

                                                                          cAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGAtgtg                                                                      >  1:995124/1‑151 (MQ=255)

                                                                          cAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGAtgtg                                                                      >  1:1657524/1‑151 (MQ=255)

                                                                           aGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATgtgt                                                                     >  2:1346590/1‑151 (MQ=255)

                                                                           aGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGAGTTGCAGCAGATgtgt                                                                     >  1:1408588/1‑151 (MQ=255)

                                                                             cTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGttt                                                                   >  2:1102161/1‑151 (MQ=255)

                                                                              tAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGtttt                                                                  >  1:731705/1‑151 (MQ=255)

                                                                              tAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGtttt                                                                  >  2:750675/1‑151 (MQ=255)

                                                                               aaTGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCGTTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTtc                                                                 <  2:1692545/151‑1 (MQ=255)

                                                                                        cGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccgtc                                                        >  1:1607604/1‑151 (MQ=255)

                                                                                         gCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccgtcg                                                       >  2:1266222/1‑151 (MQ=255)

                                                                                           cACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCgtcgccgtcgcc                                                     >  1:892562/1‑151 (MQ=255)

                                                                                             ccTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAt                                                   >  1:958393/1‑151 (MQ=255)

                                                                                                ctTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAtttt                                                >  1:1455975/1‑151 (MQ=255)

                                                                                                  ttttGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAtttt                                                <  1:1540605/149‑1 (MQ=255)

                                                                                                  ttttGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAttt                                                 <  1:1047201/148‑1 (MQ=255)

                                                                                                     tGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTgg                                           >  2:1631017/1‑151 (MQ=255)

                                                                                                      gCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgc                                          >  2:531893/1‑151 (MQ=255)

                                                                                                       cccGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgcg                                         >  2:319689/1‑151 (MQ=255)

                                                                                                        ccGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgcgc                                        >  2:725144/1‑151 (MQ=255)

                                                                                                         cGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGcgcg                                       >  2:1664920/1‑151 (MQ=255)

                                                                                                          gTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGcgcg                                       >  2:335774/1‑150 (MQ=255)

                                                                                                           tAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGAt                                     >  1:1054148/1‑151 (MQ=255)

                                                                                                           tAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGAt                                     >  2:822452/1‑151 (MQ=255)

                                                                                                             gTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCt                                   >  2:838373/1‑151 (MQ=255)

                                                                                                             gTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCt                                   >  1:1231839/1‑151 (MQ=255)

                                                                                                               tttGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCttt                                 >  2:1601194/1‑151 (MQ=255)

                                                                                                                ttGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTc                                >  2:742019/1‑151 (MQ=255)

                                                                                                                 tGTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTc                                <  2:755488/150‑1 (MQ=255)

                                                                                                                  gTAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGc                              >  1:1451568/1‑151 (MQ=255)

                                                                                                                   tAATCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTTGCCATTTTGTTGGCGCGATCTTTca                               <  1:1084345/149‑2 (MQ=255)

                                                                                                                    aaTCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTg                            >  1:1309200/1‑151 (MQ=255)

                                                                                                                    aaTCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTg                            >  2:70243/1‑151 (MQ=255)

                                                                                                                     aTCCTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGc                           >  2:658170/1‑151 (MQ=255)

                                                                                                                       ccTTCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGcc                          >  2:748852/1‑150 (MQ=255)

                                                                                                                         ttCTGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGc                        >  2:609211/1‑150 (MQ=255)

                                                                                                                            tGTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTg                    >  2:1312971/1‑151 (MQ=255)

                                                                                                                             gTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTg                    <  1:112595/150‑1 (MQ=255)

                                                                                                                             gTCGCGACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTg                    <  2:1208486/150‑1 (MQ=255)

                                                                                                                                  gACCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTAc              >  2:436017/1‑151 (MQ=255)

                                                                                                                                   aCCCACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACa             <  2:31083/151‑1 (MQ=255)

                                                                                                                                    cccACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAg            <  2:266935/151‑1 (MQ=255)

                                                                                                                                     ccACTGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGc           >  2:476016/1‑151 (MQ=255)

                                                                                                                                         tGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCggg        >  2:360837/1‑150 (MQ=255)

                                                                                                                                         tGTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCggg        >  2:268800/1‑150 (MQ=255)

                                                                                                                                          gTAACGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGa       <  1:571457/150‑1 (MQ=255)

                                                                                                                                            aaCGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAAc     >  1:1063433/1‑150 (MQ=255)

                                                                                                                                             aCGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAAc     >  1:843697/1‑149 (MQ=255)

                                                                                                                                              cGGCGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAgg  >  1:853533/1‑151 (MQ=255)

TGGTGCTGACATCATATTTTTCCACCAGCTTTTTCTCCCGGTCATAGGTGTGCTGAACATTGGTCAGGAAGTGCCAGCTAATGTCCATCGCCACCTCTTTTTGCCCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGG  >  NC_000913/1245870‑1246162

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`