BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,608,861	C→T	15.7%	A72A (GCC→GCT)	yjjG →	dUMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,608,861	0	C	T	15.7%	2004.9 / 306.4	779	A72A (GCC→GCT)	yjjG	dUMP phosphatase
Reads supporting (aligned to +/- strand): ref base C (328/328); new base T (67/56); total (395/384)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

Only 100 of 797 total aligned reads displayed.

TGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTG  >  NC_000913/4608713‑4609009

       gATGTTTCTTGATTACAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGc                                                                                                                                               <  2:5086576/149‑1 (MQ=255)

        aTGTTTCTTGATTACAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCt                                                                                                                                              >  2:3999486/1‑149 (MQ=255)

           tttCTTGATTACAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCTGAACGGCTGaa                                                                                                                                           >  4:827547/1‑149 (MQ=255)

            ttCTTGATTACAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCTGTTAACAAGCCACTGGGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAAc                                                                                                                                          >  4:3195749/1‑149 (MQ=255)

             tCTTGATTACAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACTGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACg                                                                                                                                         >  2:4191514/1‑149 (MQ=255)

               ttGATTACAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTc                                                                                                                                       <  3:3474531/149‑1 (MQ=255)

               ttGATTACAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTc                                                                                                                                       >  2:6884404/1‑149 (MQ=255)

               ttGATTACAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTc                                                                                                                                       <  4:3608278/149‑1 (MQ=255)

                   ttACAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGc                                                                                                                                   >  2:4901157/1‑149 (MQ=255)

                   ttACAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGc                                                                                                                                   >  3:37180/1‑149 (MQ=255)

                     aCAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCa                                                                                                                                 >  3:767723/1‑149 (MQ=255)

                      cAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAg                                                                                                                                <  3:5638879/149‑1 (MQ=255)

                      cAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAg                                                                                                                                <  2:5849945/149‑1 (MQ=255)

                            cACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTATCAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAAc                                                                                                                          <  4:5602876/149‑1 (MQ=255)

                            cACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAAc                                                                                                                          <  1:7081728/149‑1 (MQ=255)

                                tttACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCaa                                                                                                                      >  1:4020511/1‑149 (MQ=255)

                                tttACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCaa                                                                                                                      >  2:2952752/1‑149 (MQ=255)

                                 ttACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAAc                                                                                                                     <  2:3547974/149‑1 (MQ=255)

                                 ttACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAAc                                                                                                                     <  1:6967844/149‑1 (MQ=255)

                                          aaGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCttt                                                                                                            <  2:1096531/149‑1 (MQ=255)

                                          aaGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCttt                                                                                                            <  1:6395142/149‑1 (MQ=255)

                                           aGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTtta                                                                                                           <  4:2815165/149‑1 (MQ=255)

                                           aGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTtta                                                                                                           <  4:6383321/149‑1 (MQ=255)

                                            gATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTttat                                                                                                          <  2:1296137/149‑1 (MQ=255)

                                            gATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTAGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTttat                                                                                                          <  3:4645064/149‑1 (MQ=255)

                                              ttttCAGGACTATCAGGCCGTTAACAAGCCGCTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTttatta                                                                                                        >  1:3462396/1‑149 (MQ=255)

                                              ttttCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTttatta                                                                                                        >  1:1387381/1‑149 (MQ=255)

                                               tttCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTaa                                                                                                       <  3:970024/149‑1 (MQ=255)

                                               tttCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGACGAGCCAGGTAAACTCAACGAAGCCTTTATTaa                                                                                                       <  3:1572953/149‑1 (MQ=255)

                                               tttCAGGACTATCAGGCCGTTAACAAGCCACGGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTaa                                                                                                       <  4:5375608/149‑1 (MQ=255)

                                                ttCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAAt                                                                                                      >  2:658549/1‑149 (MQ=255)

                                                ttCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAAt                                                                                                      <  1:674109/149‑1 (MQ=255)

                                                   aGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCg                                                                                                   >  1:2961809/1‑149 (MQ=255)

                                                    ggACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGa                                                                                                  >  4:1169308/1‑149 (MQ=255)

                                                    ggACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGa                                                                                                  <  3:1064428/149‑1 (MQ=255)

                                                      aCTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATg                                                                                                <  2:3692728/149‑1 (MQ=255)

                                                        tATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGc                                                                                              <  2:606366/149‑1 (MQ=255)

                                                        tATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGa                                                                                                                                     <  1:5275103/110‑1 (MQ=255)

                                                             ggCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGaaa                                                                                         >  4:6868144/1‑149 (MQ=255)

                                                              gCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAAt                                                                                        <  3:6841294/149‑1 (MQ=255)

                                                               ccGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATc                                                                                       <  2:1932073/149‑1 (MQ=255)

                                                                 gTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCTCTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTg                                                                                     >  3:4368153/1‑149 (MQ=255)

                                                                 gTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACCCAACGAAGCCTttat                                                                                                          <  1:5549594/128‑1 (MQ=255)

                                                                       aaGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATg                                                                                                >  4:4282260/1‑132 (MQ=255)

                                                                       aaGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGc                                                                                              <  3:3911428/134‑1 (MQ=255)

                                                                       aaGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCAcgcc                                                                               >  3:6253191/1‑149 (MQ=255)

                                                                        aGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTtta                                                                                                           <  3:2994224/120‑1 (MQ=255)

                                                                         gCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCAcgccgc                                                                             >  4:1316214/1‑149 (MQ=255)

                                                                         gCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGcc                                                                                                               >  3:1236329/1‑115 (MQ=255)

                                                                          ccACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAAc                                                                                                                     >  2:4741499/1‑108 (MQ=255)

                                                                            aCTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGc                                                                          >  1:1856435/1‑149 (MQ=255)

                                                                              tgtgGGTGGATGATCAAAACGAAGAGATCACTTGATTACAGCTTTAGCACGGGCGGTTTGAGAGCTAGGCCGAACGGCTGAACGTCGAGCGAAGTAAACTCAACGAATCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCg                                                                        >  4:6537556/1‑149 (MQ=255)

                                                                               gtgggtggATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCgg                                                                       <  2:7165962/149‑1 (MQ=255)

                                                                                   gtggATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGgcgc                                                                   >  4:6146134/1‑149 (MQ=255)

                                                                                    tggATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGcc                                                                         >  4:496535/1‑142 (MQ=255)

                                                                                    tggATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGcc                                                                         >  2:4425527/1‑142 (MQ=255)

                                                                                    tggATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGcc                                                                         <  1:4425527/142‑1 (MQ=255)

                                                                                     ggATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTttat                                                                                                          <  4:1680088/108‑1 (MQ=255)

                                                                                      gATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGt                                                                >  3:4216519/1‑149 (MQ=255)

                                                                                          aTCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTtct                                                            <  3:5168217/149‑1 (MQ=255)

                                                                                           tCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTtctc                                                           <  4:5844831/149‑1 (MQ=255)

                                                                                           tCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCtt                                                                                                             >  3:5487008/1‑99 (MQ=255)

                                                                                           tCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTtctc                                                           >  4:379741/1‑149 (MQ=255)

                                                                                           tCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTtctc                                                           >  1:3626686/1‑149 (MQ=255)

                                                                                            cAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTctct                                                          <  3:3130650/149‑1 (MQ=255)

                                                                                            cAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTctct                                                          <  3:2669053/149‑1 (MQ=255)

                                                                                               aaCGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCt                                                       >  3:6535632/1‑149 (MQ=255)

                                                                                               aaCGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCt                                                       <  2:5643747/149‑1 (MQ=255)

                                                                                                 cGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTa                                                     <  4:4776513/149‑1 (MQ=255)

                                                                                                   gcgcGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCt                                                                                      >  3:1875435/1‑114 (MQ=255)

                                                                                                            cTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAATGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGt                                          >  2:6291162/1‑149 (MQ=255)

                                                                                                            cTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTctct                                                          <  3:2089515/133‑1 (MQ=255)

                                                                                                            cTTCATTACAGCTTCAGCACGGGCGGTTAGAGAGCTGGGCCGAACGGCTCAACGTCGAGCCAGGTAAACTCAACGCAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCTGTTTCTCTGCTTAACGCCATTCGt                                          <  3:6017123/149‑1 (MQ=255)

                                                                                                             ttCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTg                                                         <  4:25519/133‑1 (MQ=255)

                                                                                                             ttCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTg                                         <  1:4187029/149‑1 (MQ=255)

                                                                                                             ttCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTg                                         >  4:2085171/1‑149 (MQ=255)

                                                                                                             ttCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTg                                         >  1:2593918/1‑149 (MQ=255)

                                                                                                              tCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTgg                                        <  1:1033731/149‑1 (MQ=255)

                                                                                                                 ttACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCaa                                     <  1:816834/149‑1 (MQ=255)

                                                                                                                       cTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCaa                               >  1:2255291/1‑149 (MQ=255)

                                                                                                                         tCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCaaaa                             <  4:4040850/149‑1 (MQ=255)

                                                                                                                             cACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCgg                         >  1:2132336/1‑149 (MQ=255)

                                                                                                                              aCGGGCGGTTTGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGc                        >  1:1636368/1‑149 (MQ=255)

                                                                                                                              aCGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGc                        <  4:2286086/149‑1 (MQ=255)

                                                                                                                              aCGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGc                        <  2:2163256/149‑1 (MQ=255)

                                                                                                                              aCGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGc                        <  2:2731916/149‑1 (MQ=255)

                                                                                                                               cGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGca                       <  2:873428/149‑1 (MQ=255)

                                                                                                                                gggCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGcat                      >  4:5137977/1‑149 (MQ=255)

                                                                                                                                  gCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGcc                                                                                                               >  1:4864612/1‑58 (MQ=255)

                                                                                                                                    ggTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGCatcatc                  <  1:376748/149‑1 (MQ=255)

                                                                                                                                     gTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGCAtcatca                 <  3:4870456/149‑1 (MQ=255)

                                                                                                                                      tttGAGAGCTGGGCTGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGCATCATCAc                <  1:578528/149‑1 (MQ=255)

                                                                                                                                      tttGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTgg                                        <  1:2824398/125‑1 (MQ=255)

                                                                                                                                      tttGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGCATCATCAc                <  4:3003349/149‑1 (MQ=255)

                                                                                                                                         gagaGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGCATCATCACCaa             <  3:2320798/149‑1 (MQ=255)

                                                                                                                                              cTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCt        >  2:2474271/1‑149 (MQ=255)

                                                                                                                                                  gCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAg    <  1:6045618/149‑1 (MQ=255)

                                                                                                                                                   ccGCACGGCTGAACGCCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCCGCTTAACGCCATTCGCGGCAACGCCAAAATCGGCACCATCACCAACGGCTTTAGt   <  2:3128928/149‑1 (MQ=255)

                                                                                                                                                   ccGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGt   <  1:471906/149‑1 (MQ=255)

                                                                                                                                                   ccGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGt   <  2:1836506/149‑1 (MQ=255)

TGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGAAATCTGCACGCCGCTGCCGGGCGCGGTTTCTCTGCTTAACGCCATTCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTG  >  NC_000913/4608713‑4609009

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