BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,392,394	C→G	6.4%	noncoding (35/76 nt)	glyV →	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,394	0	C	G	6.4%	411.3 / 16.4	125	noncoding (35/76 nt)	glyV	tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base C (50/67); new base G (4/4); total (54/71)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 126 total aligned reads displayed.

AAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC  >  NC_000913/4392246‑4392540

aaaaaCGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                                                                                                                                    >  4:510811/1‑149 (MQ=255)

    aCGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAg                                                                                                                                                >  4:742357/1‑149 (MQ=255)

    aCGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAg                                                                                                                                                >  2:511350/1‑149 (MQ=255)

     cccAAAAAAGTTTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCTAAgg                                                                                                                                               >  1:962486/3‑149 (MQ=255)

      gCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGt                                                                                                                                              <  1:376863/149‑1 (MQ=255)

        aaaaaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCg                                                                                                                                            >  3:145160/1‑149 (MQ=255)

         aaaaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                           <  2:1014951/149‑1 (MQ=255)

          aaaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                          >  2:623143/1‑149 (MQ=255)

          aaaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                          <  4:394080/149‑1 (MQ=255)

          aaaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                          <  4:376419/149‑1 (MQ=255)

           aaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                                                                                                                                         >  2:4028/1‑149 (MQ=255)

           aaaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                                                                                                                                         >  3:946231/1‑149 (MQ=255)

            aaaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGACCCCGACCTTGGCAAGGTCgtgct                                                                                                                                        <  2:177705/149‑5 (MQ=255)

             aaaCTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTc                                                                                                                                       >  1:13981/1‑149 (MQ=255)

                 ttttttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgtgc                                                                                                                                         >  3:246963/1‑140 (MQ=255)

                 ttttttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                          >  3:298626/1‑142 (MQ=255)

                 ttttttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                          <  4:298626/142‑1 (MQ=255)

                 ttttttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                                                   <  2:116331/149‑1 (MQ=255)

                  tttttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAg                                                                                                                                  <  4:542992/149‑1 (MQ=255)

                    tttGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGtt                                                                                                                                <  2:183941/149‑1 (MQ=255)

                       ggggggTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                             <  3:421306/149‑1 (MQ=255)

                       ggggggTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                             <  3:997952/149‑1 (MQ=255)

                       ggggggTTGCAGAGGGAAAGATTCCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGTCCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                             <  4:731751/149‑1 (MQ=255)

                         ggggTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                                                                                                                           >  2:370035/1‑149 (MQ=255)

                         ggggTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGACCCCGACCTTGGCAAGGTCgtgctc                                                                                                                                       <  1:586078/137‑6 (MQ=255)

                                 agagGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                                                                   <  2:253647/149‑1 (MQ=255)

                                 agagGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTGGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                                                                   >  2:883792/1‑149 (MQ=255)

                                  gagGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                                                                                                  <  2:121311/149‑1 (MQ=255)

                                  gagGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                                                                                                  <  2:66243/149‑1 (MQ=255)

                                  gagGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                                                                                                  <  1:402495/149‑1 (MQ=255)

                                     ggAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTCCGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGc                                                                                                               >  1:686144/1‑149 (MQ=255)

                                      gAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCt                                                                                                              >  3:1022173/1‑149 (MQ=255)

                                       aaaGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                                                             >  4:828734/1‑149 (MQ=255)

                                       aaaGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                                                             >  2:327728/1‑149 (MQ=255)

                                        aaGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgtgc                                                                                                                                         <  3:251267/120‑4 (MQ=255)

                                        aaGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgtgc                                                                                                                                         >  4:251267/1‑117 (MQ=255)

                                        aaGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                             <  3:44760/132‑1 (MQ=255)

                                        aaGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                             >  4:44760/1‑132 (MQ=255)

                                            tttCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaaaa                                                                                                        <  1:253104/149‑1 (MQ=255)

                                               ctcGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAttt                                                                                                     <  4:376454/149‑1 (MQ=255)

                                                tcGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTg                                                                                                    >  2:948419/1‑149 (MQ=255)

                                                    ataATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaaaa                                                                                                <  3:165200/149‑1 (MQ=255)

                                                      aaTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGt                                                                                              <  3:150487/149‑1 (MQ=255)

                                                       aTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTTTGAAAAGtg                                                                                             >  2:301559/1‑149 (MQ=255)

                                                       aTGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAgag                                                                                             >  4:162887/1‑147 (MQ=255)

                                                        tGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAACAGtgc                                                                                            >  2:584000/1‑149 (MQ=255)

                                                          cgcCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtgctg                                                                                          <  3:742357/149‑1 (MQ=255)

                                                          cgcCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGtgctg                                                                                          <  2:920834/149‑1 (MQ=255)

                                                              tCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCaaa                                                                                      <  2:205686/149‑1 (MQ=255)

                                                                ccGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGc                                                                                    >  4:918030/1‑149 (MQ=255)

                                                                  gTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcac                                                                                  <  2:85359/149‑1 (MQ=255)

                                                                  gTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcac                                                                                  >  2:388770/1‑149 (MQ=255)

                                                                   tAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcaca                                                                                 <  1:950375/149‑1 (MQ=255)

                                                                     acgacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGACCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTGTCGTTTCCCGATCCAAAATTTGAAAAGTGCTGCAAAGCACAGa                                                                               <  4:13472/149‑1 (MQ=255)

                                                                      cgacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGac                                                                              <  2:482869/149‑1 (MQ=255)

                                                                       gacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGacc                                                                             >  3:23347/1‑149 (MQ=255)

                                                                       gacgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGacc                                                                             >  1:812687/1‑149 (MQ=255)

                                                                        acgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGacca                                                                            >  2:991183/1‑149 (MQ=255)

                                                                        acgCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGacca                                                                            >  3:475996/1‑149 (MQ=255)

                                                                                 tGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGAACACCCAAGCgg                                                                   <  3:839156/149‑1 (MQ=255)

                                                                                    gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGaa                                                                <  1:59630/149‑1 (MQ=255)

                                                                                    gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGaa                                                                >  2:351890/1‑149 (MQ=255)

                                                                                    gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGaa                                                                <  3:789688/149‑1 (MQ=255)

                                                                                    gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTCCCAAGGTGGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGaa                                                                <  2:758928/149‑1 (MQ=255)

                                                                                       aaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAg                                                             <  1:51518/149‑1 (MQ=255)

                                                                                       aaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGGGCTGCAAAGCACAGACCACCCAAGCGGGAATAg                                                             <  3:762690/149‑1 (MQ=255)

                                                                                        aTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGc                                                            >  4:254796/1‑149 (MQ=255)

                                                                                        aTTACGAAAGCAAAATTAAGTAGTAATCGCGAATAGTTCAGTTGGTAGAGCACGACCTTGCCACGGTCGTGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAATTGCTGCAAAGTACTGACCACCCAAGCGGGAATAGc                                                            <  4:651985/149‑1 (MQ=17)

                                                                                         ttACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCt                                                           <  1:948419/149‑1 (MQ=255)

                                                                                           aCGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCa                                                         <  4:476580/149‑1 (MQ=255)

                                                                                             gAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGt                                                       <  1:370035/149‑1 (MQ=255)

                                                                                                aGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTgg                                                    >  3:937176/1‑149 (MQ=255)

                                                                                                 gCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGATGGt                                                   >  2:457515/1‑149 (MQ=255)

                                                                                                  cAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTa                                                  <  1:137886/149‑1 (MQ=255)

                                                                                                    aaaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTATTGAAAATCACCGACCACCCAAGCGGGAATAGCTGATTTgtttga                                                >  2:839113/1‑144 (MQ=17)

                                                                                                        tAAGTAGTGCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCGCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAACGGGAATAGCTCAGTTGGTAGAGCAc                                            <  3:652563/149‑1 (MQ=255)

                                                                                                            tagtaCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGAcc                                        <  3:918030/149‑1 (MQ=255)

                                                                                                               taCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTg                                     <  1:272591/149‑1 (MQ=255)

                                                                                                                aCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                    <  1:327728/149‑1 (MQ=255)

                                                                                                                     ggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAgg                               <  4:475996/149‑1 (MQ=39)

                                                                                                                      gAATAGCGCAGTAGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAACGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTGGAGCGCGACCTTGCCAAgga                              >  4:887103/1‑148 (MQ=11)

                                                                                                                             tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTc                       >  3:803460/1‑149 (MQ=40)

                                                                                                                             tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTc                       <  4:179820/149‑1 (MQ=39)

                                                                                                                             tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCCCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                    >  3:706910/1‑136 (MQ=255)

                                                                                                                             tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGCTCGCGAGTTCGAGTCTCGTTTCCCGCCCCAAAATTTGAAAAGTGCTGCAATGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                    <  4:706910/136‑1 (MQ=255)

                                                                                                                                 ttGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGGGCTGCAAAGCACAGACCACCCAAGCGGGAATAGTTCATTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                   >  3:468000/1‑149 (MQ=255)

                                                                                                                                    gTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGtt                <  1:351890/149‑1 (MQ=39)

                                                                                                                                     tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTc               <  4:168059/149‑1 (MQ=39)

                                                                                                                                     tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTACCGCTCCAAAATTTGAAAGGGGCTGAAAAGCACAGACCACCCAAGCGGGAATAGATAAGTGGTTAGAGCACGACCTTGCTAGGGACGGGGTCGCAAGTTc               >  2:858276/1‑149 (MQ=11)

                                                                                                                                     tAGAGCACGACCTTGCCAAAGTCGTGGTCGCGAGTTCGAGTCTCGATTCCCGCTCCAAAATTTGAAAATTGTTGCAAAGCCCAGACTACGCAATCTGGAATAGCTCATTTTGTAGGGCTCGACCTTGCCAAGTTCGGGGTCGCTAGTTc               >  4:357649/1‑149 (MQ=14)

                                                                                                                                       gagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa             <  1:388770/149‑1 (MQ=40)

                                                                                                                                        agCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGcaca                                                                                 >  1:845859/1‑80 (MQ=255)

                                                                                                                                        accccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg            >  4:595564/5‑149 (MQ=255)

                                                                                                                                         gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGAGCTGCAAAGCACAGGCGACCCAAGCGGGAATAGCTCGGTTGGTAGAGCACGACCATGCCAAGCTCGGGAGCGCTAGTTCGAgc           >  4:383867/1‑148 (MQ=14)

                                                                                                                                            ctaCCTTACCAAGGTCGGGGTCGCGAGTTCGAGTCCCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACTGCGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc        <  4:532341/147‑1 (MQ=14)

                                                                                                                                            cGACCTTGCCAAGGTCGGGGTCGCGAGTTCAAGTCTCCTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc        <  4:415766/149‑1 (MQ=255)

                                                                                                                                               ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTagag                                               <  1:150690/107‑1 (MQ=255)

                                                                                                                                               ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTagag                                               >  2:150690/1‑107 (MQ=255)

                                                                                                                                               ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCCGCAAAGCACAGACCACCCAAGCGGGAATAGCACAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGtt     >  2:872597/1‑149 (MQ=255)

                                                                                                                                                  tGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAATGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc  <  1:1052525/149‑1 (MQ=38)

AAAAACGCAAAAAAAACTTTTTTGGGGGGTTGCAGAGGGAAAGATTTCTCGTATAATGCGCCTCCCGTAACGACGCAGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCC  >  NC_000913/4392246‑4392540

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