BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,993,338	Δ1 bp	54.1%	coding (2186/2547 nt)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,338	0	T	.	54.1%	198.9 / 278.3	111	coding (2186/2547 nt)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base T (24/27); new base . (34/26); total (58/53)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

Only 100 of 111 total aligned reads displayed.

TTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGA  >  NC_000913/3993195‑3993482

ttCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACcg                                                                                                                                              >  1:24724/1‑148 (MQ=255)

      tCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAGATTACCGGCggtgg                                                                                                                                       >  3:392198/1‑149 (MQ=255)

      tCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCggtg                                                                                                                                        >  1:216376/1‑148 (MQ=255)

      tCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCggtggt                                                                                                                                      <  2:770693/148‑1 (MQ=255)

       cAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCggtgg                                                                                                                                       >  1:909586/1‑148 (MQ=255)

        aaaCCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGgtggtg                                                                                                                                     >  4:35101/1‑149 (MQ=255)

        aaaCCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGa                                                                                                                                   >  3:888421/1‑149 (MQ=255)

         aaCCTGGGGGTTGTTCTTCGATCGCCTGAATATATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCATGTCAAA‑‑ACCGGCGGTGGTGGAc                                                                                                                                  >  3:899582/1‑149 (MQ=255)

         aaCCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGtggtgg                                                                                                                                    >  4:635710/1‑149 (MQ=255)

         aaCCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGa                                                                                                                                   <  1:304737/148‑1 (MQ=255)

              gggggTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCt                                                                                                                              >  2:387800/1‑148 (MQ=255)

                 ggTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTg                                                                                                                           <  1:620026/148‑1 (MQ=255)

                  gttgttCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTg                                                                                                                           >  3:660977/1‑149 (MQ=255)

                    tgttCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCCGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGcc                                                                                                                         <  4:660977/149‑1 (MQ=255)

                     gttCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGcc                                                                                                                         <  1:961748/148‑1 (MQ=255)

                     gttCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCa                                                                                                                        >  3:686083/1‑149 (MQ=255)

                      ttcttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCa                                                                                                                        >  2:965835/1‑148 (MQ=255)

                       tcttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGc                                                                                                                      >  3:325605/1‑149 (MQ=255)

                        cttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGaa                                                                                                                   >  4:431914/1‑149 (MQ=255)

                         ttcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGa                                                                                                                    >  4:476440/1‑149 (MQ=255)

                          tcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGa                                                                                                                    >  1:368524/1‑148 (MQ=255)

                          tcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAg                                                                                                                  >  2:945987/1‑148 (MQ=255)

                          tcGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAg                                                                                                                  <  2:83755/148‑1 (MQ=255)

                              aCGCCTTAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAggg                                                                                                                >  2:518743/1‑148 (MQ=255)

                              aCGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGa                                                                                                               <  3:241770/149‑1 (MQ=255)

                              aCGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGat                                                                                                              >  2:863539/1‑148 (MQ=255)

                               cGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGat                                                                                                              <  4:877477/149‑1 (MQ=255)

                                 ccTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGatc                                                                                                             <  1:518743/148‑1 (MQ=255)

                                 ccTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGatcat                                                                                                           >  1:685734/1‑148 (MQ=255)

                                  cTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATcc                                                                                                         >  3:112562/1‑149 (MQ=255)

                                   tGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACTGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGatcat                                                                                                           >  2:822233/1‑148 (MQ=255)

                                   tGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGatcat                                                                                                           <  2:664097/148‑1 (MQ=255)

                                   tGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGCAGGGATCATcc                                                                                                         >  2:23010/1‑148 (MQ=255)

                                    gAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAg                                                                                                       >  4:639297/1‑149 (MQ=255)

                                      aTGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGt                                                                                                      >  2:695433/1‑148 (MQ=255)

                                          aTCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttt                                                                                                 >  3:36383/1‑149 (MQ=255)

                                           tCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATcc                                                                                                         >  1:643000/1‑140 (MQ=255)

                                           tCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATcc                                                                                                         <  2:643000/140‑1 (MQ=255)

                                           tCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGCTCttt                                                                                                 <  2:685734/148‑1 (MQ=255)

                                            cGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCAGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttttt                                                                                               <  4:33034/149‑1 (MQ=255)

                                             ggTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttt                                                                                                 <  2:606211/148‑1 (MQ=255)

                                             ggTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttt                                                                                                 <  1:835198/148‑1 (MQ=255)

                                                 cAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCga                                                                                            >  3:745923/1‑149 (MQ=255)

                                                    aaaCTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCgaagaa                                                                                        >  2:382332/1‑148 (MQ=255)

                                                      aCTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCgaagaa                                                                                        <  1:960668/148‑1 (MQ=255)

                                                        tGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACg                                                                                     >  4:863580/1‑149 (MQ=255)

                                                                     cGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAAt                                                                         >  2:230153/1‑148 (MQ=255)

                                                                        gTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCttt                                                                   <  4:602874/149‑1 (MQ=255)

                                                                          tttATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCtt                                                                    <  1:487569/148‑1 (MQ=255)

                                                                            tATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTaa                                                                 <  3:476440/149‑1 (MQ=255)

                                                                             aTGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTaa                                                                 <  2:842036/148‑1 (MQ=255)

                                                                                gcgcGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAatat                                                              <  1:169824/148‑1 (MQ=255)

                                                                                   cgATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTa                                                           >  2:767879/1‑148 (MQ=255)

                                                                                      tttCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATtct                                                     <  4:788493/149‑1 (MQ=255)

                                                                                         cGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTcgcg                                                   <  2:604510/148‑1 (MQ=255)

                                                                                           cATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGGCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcga                                                  <  4:779981/149‑1 (MQ=255)

                                                                                             tAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcgac                                                 >  2:958141/1‑148 (MQ=255)

                                                                                             tAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTcgacga                                               <  1:290970/148‑1 (MQ=255)

                                                                                             tAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaa                                              >  4:391739/1‑149 (MQ=255)

                                                                                              aaCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaaa                                             >  3:201879/1‑149 (MQ=255)

                                                                                              aaCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaa                                              <  1:292639/148‑1 (MQ=255)

                                                                                               aCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttt                                                                                                 <  2:214963/96‑1 (MQ=255)

                                                                                               aCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCttt                                                                                                 >  1:214963/1‑96 (MQ=255)

                                                                                                cAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaa                                              <  2:909586/148‑1 (MQ=255)

                                                                                                cAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGaa                                              <  1:965835/148‑1 (MQ=255)

                                                                                                  aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAg                                            <  2:623586/148‑1 (MQ=255)

                                                                                                  aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAg                                            >  1:71374/1‑148 (MQ=255)

                                                                                                  aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGaa                                                                                                                   <  2:945310/75‑1 (MQ=255)

                                                                                                  aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGaa                                                                                                                   >  1:945310/1‑75 (MQ=255)

                                                                                                  aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATc                                                             <  4:476753/129‑1 (MQ=255)

                                                                                                  aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATc                                                             >  3:476753/1‑129 (MQ=255)

                                                                                                  aaCTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCaa                                         >  3:948764/1‑149 (MQ=255)

                                                                                                     tGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCaa                                         <  1:767879/148‑1 (MQ=255)

                                                                                                      gCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCgg                                     <  3:289436/149‑1 (MQ=255)

                                                                                                      gCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCg                                      <  1:23010/148‑1 (MQ=255)

                                                                                                             cTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTg                                 <  1:940647/148‑1 (MQ=255)

                                                                                                             cTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTg                                 <  2:24724/148‑1 (MQ=255)

                                                                                                              tGTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAAACGGGTTGAGGt                             >  4:769357/1‑149 (MQ=255)

                                                                                                               gTCAGTGCAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGt                             >  2:647231/1‑148 (MQ=255)

                                                                                                                  aGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGta                            <  1:132379/148‑1 (MQ=255)

                                                                                                                  aGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGta                            <  2:216376/148‑1 (MQ=255)

                                                                                                                      cAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcacc                     <  4:360474/149‑1 (MQ=255)

                                                                                                                      cAGGTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcac                      >  1:39385/1‑148 (MQ=255)

                                                                                                                         gTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATcacca                    <  4:366712/149‑1 (MQ=255)

                                                                                                                         gTTGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACt                  >  4:668331/1‑149 (MQ=255)

                                                                                                                         gTTGAAACCAATCACGTCAAA‑‑ACCGGC‑‑‑GGTGGAGGGCTTTGTCACCGAAGAGAGCATCCAGCACTTTTTCGAAGAAACGCAAGACGACAAGAGCGTCAATATCTACATTCTCGTCGAAAGAAACAAGCGTGAGGTATATTACCACTGct               >  4:900308/1‑148 (MQ=255)

                                                                                                                          ttGAAACCGATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAg                                                                                                                  <  1:612070/52‑1 (MQ=255)

                                                                                                                          ttGAAACCGATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAg                                                                                                                  >  2:612070/1‑52 (MQ=255)

                                                                                                                          ttGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACt                  >  1:917494/1‑148 (MQ=255)

                                                                                                                          ttGAAACCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACt                  >  2:321469/1‑148 (MQ=255)

                                                                                                                             aaaCCAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGa              >  4:615731/1‑149 (MQ=255)

                                                                                                                                ccAATCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAgg           <  4:147478/149‑1 (MQ=255)

                                                                                                                                 cAATCACGTCAAA‑‑ACCGGCGGTGGTGTACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGgc          >  4:30780/1‑149 (MQ=255)

                                                                                                                                  aaTCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGgca         <  4:875850/149‑1 (MQ=255)

                                                                                                                                   aTCACGTCAAA‑‑ACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCg                                      <  2:798930/119‑1 (MQ=255)

                                                                                                                                        gTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAgg                                                                                                                 <  4:467159/41‑1 (MQ=255)

                                                                                                                                        gTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAgg                                                                                                                 >  3:467159/1‑41 (MQ=255)

                                                                                                                                         tCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCaaa    <  3:664919/149‑1 (MQ=255)

                                                                                                                                          cAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCaaa    >  1:87065/1‑148 (MQ=255)

                                                                                                                                            aaTTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAaga  <  2:763742/148‑1 (MQ=255)

TTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGCATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGA  >  NC_000913/3993195‑3993482

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