BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,392,582	A→C	7.0%	intergenic (+35/‑1)	glyX → / → glyY	tRNA‑Gly/tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,582	0	A	C	7.0%	517.8 / 8.1	189	intergenic (+35/‑1)	glyX/glyY	tRNA‑Gly/tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base A (120/51); new base C (8/9); total (129/60)
Fisher's exact test for biased strand distribution p-value = 6.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

CAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTC  >  NC_000913/4392434‑4392731

cAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTG‑CAAGGGCTGTAAGGCACAGACCACCCCAg                                                                                                                                                      >  2:3881851/1‑150 (MQ=255)

 aaaaTTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGc                                                                                                                                                     >  2:1222923/1‑150 (MQ=255)

 aaaaTTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGc                                                                                                                                                     >  1:2731570/1‑150 (MQ=255)

 aaaaTTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGc                                                                                                                                                     >  2:2803377/1‑150 (MQ=255)

  aaaTTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAcgcg                                                                                                                                                    >  1:1061746/1‑150 (MQ=255)

   aaTTTGAAAAGTTCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCgg                                                                                                                                                   <  2:152262/150‑1 (MQ=255)

     tttGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCACGCGGGa                                                                                                                                                 >  1:1544086/1‑150 (MQ=255)

      ttGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGaa                                                                                                                                                >  2:3119664/1‑150 (MQ=255)

       tGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAAt                                                                                                                                               >  2:2731903/1‑150 (MQ=255)

         aaaaGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAATGCCCAGACCACCCAAGCGGGAATAg                                                                                                                                             >  2:1292253/1‑150 (MQ=255)

         aaaaGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAg                                                                                                                                             >  2:420820/1‑150 (MQ=255)

         aaaaGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTG‑AAAGTGATGTAATGAACAGACCACCCAAGAGGGAATAg                                                                                                                                             >  2:420838/1‑150 (MQ=18)

          aacgTGCTGCAACGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGc                                                                                                                                            <  2:2335154/147‑1 (MQ=255)

             gTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTG‑AAAGTGCTGTAAGGCACACACCACCCAAGCGGGAATAGCTcc                                                                                                                                         >  2:2233558/1‑149 (MQ=255)

                 tgcAAAGCACAGACCACCCAAGCGGGAATAGCTCAGGTGGTAGAGCACGACCTTGCCAAGGTCGGGGGCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCACGCGGGAATAGCTCAGTTg                                                                                                                                     >  1:2607568/1‑150 (MQ=255)

                   cAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGt                                                                                                                                   <  1:2149503/150‑1 (MQ=255)

                   cAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGt                                                                                                                                   <  1:2149758/150‑1 (MQ=255)

                    aaaGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTa                                                                                                                                  >  1:2411781/1‑150 (MQ=255)

                    aaaGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGAGGGAATAGCTCAGTTGGTa                                                                                                                                  >  1:1180836/1‑150 (MQ=255)

                                accCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCtt                                                                                                                      <  1:1028298/150‑1 (MQ=255)

                                 cccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTG‑AAAGTGCTGTAAGGCCCAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTg                                                                                                                     >  2:1826239/1‑150 (MQ=255)

                                  ccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCCCCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc                                                                                                                    >  1:1732020/1‑150 (MQ=255)

                                      gCGGTAATAGCTAATTTGGTGGTGCAGTACCTTGCCATGGTCGGGGTGGCTATTTCTATTCTCGTTTCCCGCTCCAAAATTTT‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAg                                                                                                                <  1:2299572/150‑1 (MQ=11)

                                      gCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCCGTTGGTAGAGCACGACCTTGCCAAg                                                                                                                >  1:1242815/1‑150 (MQ=38)

                                       cGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAgg                                                                                                               >  1:524721/1‑150 (MQ=38)

                                       cGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATCTG‑AAAGAGCTGAAAGCCACACACCACCCCCGCGCGAATAGCGCCGTTGGTAGAGCACCACCCAGCCACgg                                                                                                               >  2:4041347/1‑150 (MQ=11)

                                       cGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAgt                                                                                                               >  2:2163117/1‑149 (MQ=255)

                                            atcgCTCTGTTGGTTGTGAACTACCTGGCCAAGGTCGGGTTCGCGATTTCGAGTCTCGTTGCCCGCTCCAAAATTTT‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                          <  2:1789101/147‑1 (MQ=11)

                                            aTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                          <  2:1801791/150‑1 (MQ=38)

                                            aTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                          <  1:1260256/150‑1 (MQ=38)

                                            aTAGCCCAGCTGGTAGAGCACGACCTCGCCAAGGACGGGGCCGAGAGTGCGCGTCCCGTTTCCCGCCCCAAAATATG‑AAAGTGCTGTAAGGCACAGACCCACCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                          <  2:3659980/150‑1 (MQ=11)

                                                cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcg                                                                                                      >  2:2962394/1‑150 (MQ=38)

                                                cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcg                                                                                                      >  1:1334798/1‑150 (MQ=38)

                                                 tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                     >  1:3156002/1‑150 (MQ=38)

                                                 acAGGTGGCAGAGCACGACCTTGCCAAGGTGGGGGTCGCGAGTTCGAGTCTCGTTTCCCGATCCAAAATATG‑AAAGTGCTGTAAGGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                     <  2:1800061/149‑1 (MQ=11)

                                                  cAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgcg                                                                                                    >  1:4110235/1‑150 (MQ=38)

                                                   aGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTGGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                   >  1:3482915/1‑150 (MQ=38)

                                                   aGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGAACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACATTGCCAAGGTCGGGGTCGCGa                                                                                                   >  1:308742/1‑150 (MQ=255)

                                                     ttGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCGTGCCAAGGTCGGGGTCGCGAGt                                                                                                 >  2:2440369/1‑150 (MQ=38)

                                                       ggTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCAGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                                               >  1:703445/1‑150 (MQ=38)

                                                       ggTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCGCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCACGGTCGGGGTCGCGAGtgc                                                                                               >  1:1221288/1‑148 (MQ=255)

                                                       ggTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTAGAGTCTCGTTTCACGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTc                                                                                               >  1:825071/1‑150 (MQ=255)

                                                          agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAg                                                                                            >  1:3968240/1‑150 (MQ=38)

                                                              cACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                        >  2:1534308/1‑150 (MQ=38)

                                                               aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCg                                                                                       >  1:2393120/1‑150 (MQ=38)

                                                                       gCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGc                                                                               >  2:3963840/1‑150 (MQ=38)

                                                                         cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                             >  2:926702/1‑150 (MQ=38)

                                                                         cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                             >  2:2504182/1‑150 (MQ=38)

                                                                             gTCGGGGTCGCGAGTTCGAGTCTCGTTTACCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACACACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGGTCGAGTCGCGTTTCCCGCTCCaaa                                                                         >  1:3496395/1‑150 (MQ=255)

                                                                                      gcgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctct                                                                >  2:316778/1‑150 (MQ=255)

                                                                                      gcgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctct                                                                >  1:2964250/1‑150 (MQ=255)

                                                                                      cccacatCGAGTCACGTTACCCCCACCAAAATATC‑AAAGTCCTGAAACCCACACACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctct                                                                <  2:2729376/144‑1 (MQ=14)

                                                                                       cgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctctc                                                               >  1:1550923/1‑150 (MQ=255)

                                                                                       cgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctctc                                                               >  1:2755208/1‑150 (MQ=255)

                                                                                        gAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCa                                                              >  1:3303308/1‑150 (MQ=255)

                                                                                        gAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCa                                                              >  1:1686405/1‑150 (MQ=255)

                                                                                                 tctcGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAatat                                                     >  2:743438/1‑150 (MQ=255)

                                                                                                   tcGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaaaa                                                        <  2:1353447/146‑1 (MQ=255)

                                                                                                         cccGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCg                                             >  1:2243467/1‑150 (MQ=255)

                                                                                                           cGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGAc                                           >  1:1564789/1‑150 (MQ=255)

                                                                                                                cAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGAt                                      >  1:1875690/1‑150 (MQ=255)

                                                                                                                cAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGAt                                      >  1:2213052/1‑150 (MQ=255)

                                                                                                                 aaaaTTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATg                                     >  1:993185/1‑150 (MQ=255)

                                                                                                                 aaaaTTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATg                                     >  1:993121/1‑150 (MQ=255)

                                                                                                                  aaaTTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCACGGTGGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGc                                    <  1:3119664/150‑1 (MQ=255)

                                                                                                                     tttG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGtt                                 >  2:1335005/1‑150 (MQ=255)

                                                                                                                     tttG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGtt                                 >  2:58215/1‑150 (MQ=255)

                                                                                                                        g‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGGTAtt                              >  2:562287/1‑150 (MQ=255)

                                                                                                                           aaGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGc                            >  1:1907997/1‑150 (MQ=255)

                                                                                                                           aaGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGc                            >  1:835217/1‑150 (MQ=255)

                                                                                                                            cgTGCTGCAATGCACAGACCACGCAAGCGGGAATTGCCCAGTTGGTAGCGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCt                           <  2:2207459/149‑1 (MQ=18)

                                                                                                                              tGCTGTAAGGCACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTgg                         <  2:1665172/150‑1 (MQ=255)

                                                                                                                              tGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTgg                         >  1:2349903/1‑150 (MQ=255)

                                                                                                                              tGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTgg                         >  2:262643/1‑150 (MQ=255)

                                                                                                                              tGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTgg                         >  2:3075558/1‑150 (MQ=255)

                                                                                                                                cTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGtt                       >  1:3231953/1‑150 (MQ=255)

                                                                                                                                 tGTAAGTCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttt                      <  1:1965365/150‑1 (MQ=255)

                                                                                                                                 tGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttt                      >  1:3901933/1‑150 (MQ=255)

                                                                                                                                 tGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGgtt                      >  1:2847458/1‑150 (MQ=255)

                                                                                                                                   tAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                    >  1:599305/1‑150 (MQ=255)

                                                                                                                                    tccgCACTGACCACGCACGCGGGAATAGCTCAGTTGGTAGAGCACGCCGTTGCCACGGGCGGGGTCGCGTGTTGGAGTCGCGTTTCCGGCTCCAAATTCTTCTCTCACTAAAATATCCACCGCGACGCGATGCGTTATTGCTGGTTTttg                   <  2:11636/147‑1 (MQ=21)

                                                                                                                                    aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTttg                   <  2:264817/150‑1 (MQ=255)

                                                                                                                                      ggCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTttgtt                 <  1:1649759/150‑1 (MQ=255)

                                                                                                                                        cacaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTtgttgt               <  2:790762/150‑1 (MQ=255)

                                                                                                                                          caGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtct             <  2:304874/150‑1 (MQ=255)

                                                                                                                                          caGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtct             <  1:2233558/150‑1 (MQ=255)

                                                                                                                                          caGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtct             <  1:2680315/150‑1 (MQ=255)

                                                                                                                                           aGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtctc            <  1:339081/150‑1 (MQ=255)

                                                                                                                                              ccaccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGa         <  1:1811158/150‑1 (MQ=255)

                                                                                                                                              ccaccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGa         <  1:864460/150‑1 (MQ=255)

                                                                                                                                              ccaccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGa         >  2:490133/1‑150 (MQ=255)

                                                                                                                                               caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCAACAGCGACGCGATGCGTTATTGCTGGTGTTTGTTGTCTCTGAc        >  1:426076/1‑150 (MQ=255)

                                                                                                                                                 cccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaa      >  2:844668/1‑150 (MQ=255)

                                                                                                                                                   cAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAc    >  1:3764690/1‑150 (MQ=255)

                                                                                                                                                   cAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAc    >  2:3530315/1‑150 (MQ=255)

                                                                                                                                                   cAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAc    >  1:2411591/1‑150 (MQ=255)

                                                                                                                                                    aaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAct   <  1:3963840/150‑1 (MQ=255)

                                                                                                                                                    aaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAct   >  2:3504440/1‑150 (MQ=255)

                                                                                                                                                    aaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAct   >  2:2001696/1‑150 (MQ=255)

                                                                                                                                                     aGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAActc  <  1:436434/150‑1 (MQ=255)

CAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTC  >  NC_000913/4392434‑4392731

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`