BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,993,145	C→A	57.0%	R665S (CGT→AGT)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,145	0	C	A	57.0%	58.6 / 306.8	184	R665S (CGT→AGT)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base C (45/34); new base A (59/46); total (104/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

Only 100 of 187 total aligned reads displayed.

GATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAA  >  NC_000913/3992997‑3993290

 atgatCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCg                                                                                                                                                  >  1:461456/1‑149 (MQ=255)

  tgatCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGt                                                                                                                                                 <  2:532591/149‑1 (MQ=255)

   gatCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTc                                                                                                                                                >  4:534190/1‑149 (MQ=255)

    atCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCt                                                                                                                                               >  3:1109111/1‑149 (MQ=255)

        aaGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTcc                                                                                                                                           >  3:1017302/1‑149 (MQ=255)

         aGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCa                                                                                                                                          >  4:186227/1‑149 (MQ=255)

              tGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCAcc                                                                                                                                     >  4:680692/1‑149 (MQ=255)

                    cTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAg                                                                                                                               >  1:979566/1‑149 (MQ=255)

                     tATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGTGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAgg                                                                                                                              >  3:300441/1‑149 (MQ=255)

                     tATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAgg                                                                                                                              >  3:1432236/1‑149 (MQ=255)

                          tcGGCAAAATGCATCAGGACGCCGCACCGCTAGAGAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATACGTACTCGCGTGCAGCAACAGGTTTCTTAGTGTATTGAATTGCGTCTTTCCCGCACCCGCCAGGAAAcc                                                                                                                         <  4:901386/149‑1 (MQ=255)

                          tcGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGATATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAAcc                                                                                                                         >  1:282233/1‑149 (MQ=255)

                            ggCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCgg                                                                                                                       <  2:943106/149‑1 (MQ=255)

                            ggCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCgg                                                                                                                       >  2:425719/1‑149 (MQ=255)

                             gCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCggg                                                                                                                      >  3:1031237/1‑149 (MQ=255)

                               aaaaTGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCg                                                                                                                    >  1:303452/1‑149 (MQ=255)

                                     cATCAGGACGCCGCACCGCCAGATAGCGTGGTAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCaa                                                                                                              <  1:91017/149‑1 (MQ=255)

                                     cATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCaa                                                                                                              <  3:297383/149‑1 (MQ=255)

                                      aTCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAg                                                                                                             >  1:665517/1‑149 (MQ=255)

                                      aTCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAg                                                                                                             >  2:1022685/1‑149 (MQ=255)

                                       tCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAgg                                                                                                            >  1:263563/1‑149 (MQ=255)

                                       tCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAgg                                                                                                            >  1:550433/1‑149 (MQ=255)

                                        cAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGgc                                                                                                           >  4:452692/1‑149 (MQ=255)

                                           gACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCt                                                                                                        >  1:731909/1‑149 (MQ=255)

                                            aCGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTg                                                                                                       >  3:1397294/1‑149 (MQ=255)

                                               ccgcACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTgcgc                                                                                                    <  2:1285607/149‑1 (MQ=255)

                                               ccgcACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTgcgc                                                                                                    >  2:877827/1‑149 (MQ=255)

                                                cgcACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGcgcg                                                                                                   >  2:1284223/1‑149 (MQ=255)

                                                      gCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTg                                                                                             <  4:797998/149‑1 (MQ=255)

                                                       ccAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTgg                                                                                            <  2:665517/149‑1 (MQ=255)

                                                        cAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGt                                                                                           <  4:1031237/149‑1 (MQ=255)

                                                          gATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCa                                                                                         >  3:1014370/1‑149 (MQ=255)

                                                               cGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTCAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCt                                                                                    <  1:1022685/149‑1 (MQ=255)

                                                                  ggAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTggg                                                                                 >  1:1394626/1‑149 (MQ=255)

                                                                  ggAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTggg                                                                                 <  4:732535/149‑1 (MQ=255)

                                                                      gTCTTCTGTTATAGCGAGCATCTGCGCGGCTTAATTCGTACTAGCGTGCAGCAACTGGTTTCTGGGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGAGCTGCGCGTTTCTGGTCAAACCTGGGGgtt                                                                             <  2:676487/149‑1 (MQ=255)

                                                                      gTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgtt                                                                             <  3:308888/149‑1 (MQ=255)

                                                                       tcttctGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgttg                                                                            >  1:730990/1‑149 (MQ=255)

                                                                        cttctGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgttgt                                                                           <  4:603972/149‑1 (MQ=255)

                                                                         ttctGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTc                                                                                                                <  2:172412/111‑1 (MQ=255)

                                                                         ttctGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTc                                                                                                                >  1:172412/1‑111 (MQ=255)

                                                                            tGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGttctt                                                                       >  3:887803/1‑149 (MQ=255)

                                                                             gTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGttcttc                                                                      <  2:461456/149‑1 (MQ=255)

                                                                              ttATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGc                                                                                                                     >  3:645635/1‑101 (MQ=255)

                                                                                 taGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAAc                                                                  <  2:858558/149‑1 (MQ=255)

                                                                                   gCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGc                                                                >  1:383347/1‑149 (MQ=255)

                                                                                    ccAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTcc                                                                                                                                           >  1:1298509/1‑73 (MQ=255)

                                                                                    ccAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTcc                                                                                                                                           <  2:1298509/73‑1 (MQ=255)

                                                                                     cAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCt                                                              <  1:877827/149‑1 (MQ=255)

                                                                                     cAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgtt                                                                             <  4:245047/134‑1 (MQ=255)

                                                                                     cAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgtt                                                                             >  3:245047/1‑134 (MQ=255)

                                                                                        cATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGaa                                                           >  4:623251/1‑149 (MQ=255)

                                                                                        cATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAAGTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGaa                                                           <  2:282233/149‑1 (MQ=255)

                                                                                         aTCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCg                                                                                                                        >  1:1017732/1‑87 (MQ=255)

                                                                                         aTCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTg                                                                                             >  2:788651/1‑114 (MQ=255)

                                                                                              cgcgGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATc                                                     <  1:1264576/149‑1 (MQ=255)

                                                                                               gcgGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCg                                                    >  2:1248094/1‑149 (MQ=255)

                                                                                                   cTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTAc                                                <  2:1390634/149‑1 (MQ=255)

                                                                                                    ttAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACa                                               >  3:1308786/1‑149 (MQ=255)

                                                                                                    ttAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACa                                               >  3:1080770/1‑149 (MQ=255)

                                                                                                     tAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAg                                              <  3:212611/149‑1 (MQ=255)

                                                                                                        ttCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgttgtt                                                                          >  3:1283076/1‑118 (MQ=255)

                                                                                                        ttCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgttgtt                                                                          <  4:1283076/118‑1 (MQ=255)

                                                                                                        ttCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGaaa                                           <  3:1346758/149‑1 (MQ=255)

                                                                                                         tCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAAc                                          >  1:1468449/1‑149 (MQ=255)

                                                                                                          cGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACt                                         <  4:1014370/149‑1 (MQ=255)

                                                                                                            tACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTgg                                       >  3:427380/1‑149 (MQ=255)

                                                                                                               tCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGaaa                                           <  4:969657/142‑1 (MQ=255)

                                                                                                               tCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGaaa                                           >  3:969657/1‑142 (MQ=255)

                                                                                                               tCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGaaa                                    >  1:1034508/1‑149 (MQ=255)

                                                                                                               tCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGaaa                                    >  1:849260/1‑149 (MQ=255)

                                                                                                                 gcgTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAAc                                  <  2:730990/149‑1 (MQ=255)

                                                                                                                 gcgTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAAc                                  >  2:1478027/1‑149 (MQ=255)

                                                                                                                 gcgTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACt                                         >  3:175353/1‑142 (MQ=255)

                                                                                                                 gcgTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACt                                         <  4:175353/142‑1 (MQ=255)

                                                                                                                 gcgTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAAc                                  <  3:765903/149‑1 (MQ=255)

                                                                                                                     gcagcaACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTgg                                                                                            <  4:496692/87‑1 (MQ=255)

                                                                                                                       agcaACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGcc                               <  4:313844/146‑1 (MQ=255)

                                                                                                                       agcaACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATc                            <  3:472925/149‑1 (MQ=255)

                                                                                                                       agcaACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATc                            <  4:867929/149‑1 (MQ=255)

                                                                                                                        gcaACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGc                                >  1:1130595/1‑144 (MQ=255)

                                                                                                                            cTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAAAAACTGGAAAACGCCATCGAGtt                       >  1:812757/1‑149 (MQ=255)

                                                                                                                            cTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGtt                       <  4:1397294/149‑1 (MQ=255)

                                                                                                                             tGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGttt                      <  1:490547/149‑1 (MQ=255)

                                                                                                                             tGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGt                                                  <  2:247285/121‑1 (MQ=255)

                                                                                                                             tGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGt                                                  >  1:247285/1‑121 (MQ=255)

                                                                                                                              ggTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGtttt                     >  4:411854/1‑149 (MQ=255)

                                                                                                                               gTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTa                    >  2:226845/1‑149 (MQ=255)

                                                                                                                                tttCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTAt                   <  3:411854/149‑1 (MQ=255)

                                                                                                                                tttCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCa                              >  1:432454/1‑138 (MQ=255)

                                                                                                                                 ttCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATg                  <  2:145380/149‑1 (MQ=255)

                                                                                                                                  tCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATgg                 <  3:680692/149‑1 (MQ=255)

                                                                                                                                    tGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGgcg               >  2:250961/1‑149 (MQ=255)

                                                                                                                                          tATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGTGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGAttt         >  1:433886/1‑149 (MQ=255)

                                                                                                                                          tATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTACTGGTCAAACCTGGGGGTTGCTCTTCGAACGCCTGGATGTATAGGTACAGAAACTGGAAAACTCCATCGAGTTTTATGGCGCGAttt         >  1:1460277/1‑149 (MQ=255)

                                                                                                                                           attgaattgaGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTc        <  3:534190/149‑1 (MQ=255)

                                                                                                                                            ttGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCg       >  3:847788/1‑149 (MQ=255)

                                                                                                                                              gAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCa     <  2:264156/149‑1 (MQ=255)

                                                                                                                                                 ttGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCCGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATaa  >  2:483966/1‑149 (MQ=255)

                                                                                                                                                 ttGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTAt                                                      <  4:1103909/97‑1 (MQ=255)

GATGATCGAAGCCCTGAAAACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAA  >  NC_000913/3992997‑3993290

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