BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,991,844	A→T	81.8%	D231V (GAT→GTT)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,991,844	0	A	T	81.8%	361.5 / 66.4	137	D231V (GAT→GTT)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base A (8/17); new base T (55/57); total (63/74)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

Only 100 of 137 total aligned reads displayed.

CTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCC  >  NC_000913/3991699‑3991985

ctTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTtct                                                                                                                                            <  2:36738/149‑1 (MQ=255)

 ttGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGAtctc                                                                                                                                           >  1:399542/1‑149 (MQ=255)

    aCGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCgg                                                                                                                                         <  4:785439/148‑1 (MQ=255)

       aaTTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGc                                                                                                                                     <  1:1012675/149‑1 (MQ=255)

       aaTTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGc                                                                                                                                     >  2:1047918/1‑149 (MQ=255)

       aaTTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTgg                                                                                                                                      <  4:845476/148‑1 (MQ=255)

        aTTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGc                                                                                                                                     <  3:405464/148‑1 (MQ=255)

         ttttATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCtt                                                                                                                                   >  2:782136/1‑149 (MQ=255)

             aTCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGc                                                                                                                               >  2:977202/1‑149 (MQ=255)

              tCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCt                                                                                                                              >  2:821391/1‑149 (MQ=255)

                 tACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGATGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGc                                                                                                                           >  2:80047/1‑149 (MQ=255)

                 tACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGc                                                                                                                           >  1:450312/1‑149 (MQ=255)

                     gccgTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCttt                                                                                                                        >  3:88116/1‑148 (MQ=255)

                      ccgTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTc                                                                                                                       <  4:1052214/148‑1 (MQ=255)

                        gTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTctgc                                                                                                                    <  2:618749/149‑1 (MQ=255)

                         tGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTctgc                                                                                                                    <  3:876440/148‑1 (MQ=255)

                            gTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGaa                                                                                                                <  1:821391/149‑1 (MQ=255)

                              ctcGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAag                                                                                                               <  3:418054/148‑1 (MQ=255)

                              ctcGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAag                                                                                                               <  3:1007658/148‑1 (MQ=255)

                              ctcGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAag                                                                                                               <  3:219376/148‑1 (MQ=255)

                                cGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGt                                                                                                            <  1:367520/149‑1 (MQ=255)

                                      tAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACttt                                                                                                       <  3:432535/148‑1 (MQ=255)

                                      tAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTg                                                                                                      <  1:1015017/149‑1 (MQ=255)

                                        aGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTgg                                                                                                     <  3:524739/148‑1 (MQ=255)

                                         gCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTgg                                                                                                     <  2:734629/147‑1 (MQ=255)

                                         gCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTg                                                                                                   <  2:878702/149‑1 (MQ=255)

                                          cGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTg                                                                                                      >  3:34396/1‑145 (MQ=255)

                                          cGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTg                                                                                                   >  3:708330/1‑148 (MQ=255)

                                          cGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGc                                                                                                  >  1:763466/1‑149 (MQ=255)

                                           gTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACttt                                                                                                       >  3:692452/1‑143 (MQ=255)

                                           gTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGc                                                                                                  >  3:652795/1‑148 (MQ=255)

                                            tATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGcc                                                                                                 >  4:358816/1‑148 (MQ=255)

                                            tATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCa                                                                                                <  1:222049/149‑1 (MQ=255)

                                               tCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGcc                                                                                             <  1:246804/149‑1 (MQ=255)

                                                cTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCt                                                                                            >  2:122097/1‑149 (MQ=255)

                                                 tgtgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCtt                                                                                           >  2:29751/1‑149 (MQ=255)

                                                  gtgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGTTTTACGCGCAGGGCGTGCTGAGGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCtt                                                                                           <  3:180787/148‑1 (MQ=255)

                                                  gtgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCttt                                                                                          <  2:1043841/149‑1 (MQ=255)

                                                  gtgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCtt                                                                                           >  3:816387/1‑148 (MQ=255)

                                                  gtgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCtt                                                                                           >  3:309470/1‑148 (MQ=255)

                                                  gtgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCtt                                                                                           >  3:911007/1‑148 (MQ=255)

                                                   tgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTg                                                                                                                                       >  1:154278/1‑103 (MQ=255)

                                                   tgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTg                                                                                                                                       <  2:154278/103‑1 (MQ=255)

                                                   tgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTg                                                                                         <  2:421660/149‑1 (MQ=255)

                                                    ggAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTgg                                                                                        <  2:348879/149‑1 (MQ=255)

                                                     gAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTgg                                                                                        <  4:722643/148‑1 (MQ=255)

                                                         atGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTct                                                                                                                      <  4:139353/114‑1 (MQ=255)

                                                         atGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTct                                                                                                                      >  3:139353/1‑114 (MQ=255)

                                                          tGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCttt                                                                                          >  1:165199/1‑141 (MQ=255)

                                                                 gTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAag                                                                            <  3:227143/148‑1 (MQ=255)

                                                                 gTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAaga                                                                           <  1:953167/149‑1 (MQ=255)

                                                                    cgacgaAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGt                                                                         <  4:276334/148‑1 (MQ=255)

                                                                     gacgaAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTa                                                                        <  3:428426/148‑1 (MQ=255)

                                                                     gacgaAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTAt                                                                       >  2:697937/1‑149 (MQ=255)

                                                                      acgaAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTAt                                                                       >  4:666465/1‑148 (MQ=255)

                                                                        gaagaaGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCg                                                                     <  3:375510/148‑1 (MQ=255)

                                                                        gaagaaGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGa                                                                    <  1:404407/149‑1 (MQ=255)

                                                                          agaagaGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGAt                                                                   <  4:88116/148‑1 (MQ=255)

                                                                            aagaGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGa                                                                                                                 <  4:932462/100‑1 (MQ=255)

                                                                            aagaGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGa                                                                                                                 >  3:932462/1‑100 (MQ=255)

                                                                             agagCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTcc                                                                <  4:652795/148‑1 (MQ=255)

                                                                             agagCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTccc                                                               <  1:130147/149‑1 (MQ=255)

                                                                                    tACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACa                                                         >  4:151941/1‑148 (MQ=255)

                                                                                    tACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACa                                                         >  4:564961/1‑148 (MQ=255)

                                                                                       gacgacTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGTAGCTCTACAAGAGTATCGATTCCCCATACAAAGc                                                     >  1:694162/1‑149 (MQ=255)

                                                                                         cgacTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCg                                                    <  4:362782/148‑1 (MQ=255)

                                                                                           acTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAag                                                                            >  4:166994/1‑122 (MQ=255)

                                                                                                gtgATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGaa                                            >  2:937830/1‑149 (MQ=255)

                                                                                                  gatgaCGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGc                                                                                                  >  3:776173/1‑93 (MQ=255)

                                                                                                  gatgaCGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTCGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGc                                                                                                  <  4:776173/93‑1 (MQ=255)

                                                                                                      aCGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAacac                                       >  4:457474/1‑148 (MQ=255)

                                                                                                       cGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActg                                     <  2:848590/149‑1 (MQ=255)

                                                                                                        gCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgc                                    >  2:610218/1‑149 (MQ=255)

                                                                                                         cTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgc                                    >  3:713524/1‑148 (MQ=255)

                                                                                                                cgcAGGGCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGa                             >  4:31319/1‑148 (MQ=255)

                                                                                                                 gcAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGaa                            <  3:662929/148‑1 (MQ=255)

                                                                                                                     ggCGTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTa                       <  1:748833/149‑1 (MQ=255)

                                                                                                                        gTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTGTCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAtt                     <  3:944367/148‑1 (MQ=255)

                                                                                                                        gTGCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTc                    >  2:468212/1‑149 (MQ=255)

                                                                                                                          gCTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCt                  <  1:29751/149‑1 (MQ=255)

                                                                                                                           cTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCt                  <  3:109888/148‑1 (MQ=255)

                                                                                                                           cTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCt                  <  4:88617/148‑1 (MQ=255)

                                                                                                                           cTGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTg                 <  2:476668/149‑1 (MQ=255)

                                                                                                                           cTGACGCCAAATGAATGGCTGGTTCTCGGTAGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTg                 <  1:697937/149‑1 (MQ=255)

                                                                                                                           cTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTg                 <  2:280928/149‑1 (MQ=255)

                                                                                                                            tGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTTCCAGCCTTTTGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAGCACAGCTGCTGGAAGCCTATTCCTgg                <  2:164787/149‑1 (MQ=255)

                                                                                                                            tGACGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTg                 <  3:151941/148‑1 (MQ=255)

                                                                                                                              aCGCCAAATGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTggg               >  4:574854/1‑148 (MQ=255)

                                                                                                                                gCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATa                                                           >  2:1015379/1‑102 (MQ=255)

                                                                                                                                gCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATa                                                           <  1:1015379/102‑1 (MQ=255)

                                                                                                                                   aaaTGAATGGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAc          >  3:468947/1‑148 (MQ=255)

                                                                                                                                   aaaTGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAcc         <  1:658375/149‑1 (MQ=255)

                                                                                                                                    aatgaatgGCTGGTTGTCGGTGGCTCAAGCTCGCTTTGTGCTGAAGAGTACTTTGGTGCCAGCCTGTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTGCTGAAAACACTGCTGCTGGAAGCCTATTCCTCGGAATAcc         >  3:658159/1‑148 (MQ=255)

                                                                                                                                    aatgaatgGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc        >  1:895269/1‑149 (MQ=255)

                                                                                                                                    aatgaatgGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAcc         >  4:915561/1‑148 (MQ=255)

                                                                                                                                    aatgaatgGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAcc         >  3:53989/1‑148 (MQ=255)

                                                                                                                                     atgaatgGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc        >  4:376004/1‑148 (MQ=255)

                                                                                                                                     atgaatgGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCg       >  2:757943/1‑149 (MQ=255)

                                                                                                                                     atgaatgGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCg       >  1:353632/1‑149 (MQ=255)

                                                                                                                                          tgGCTGGTTCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAAccc  <  2:353632/149‑1 (MQ=255)

CTTGACGAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCC  >  NC_000913/3991699‑3991985

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