BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,617	0	C	G	10.6%	319.1 / 26.1	103	noncoding (35/76 nt)	glyY	tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base C (45/47); new base G (6/5); total (51/52)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

Only 100 of 104 total aligned reads displayed.

GCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGT  >  NC_000913/4392472‑4392764

gCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCaa                                                                                                                                                  >  1:323969/1‑149 (MQ=38)

 cGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCaa                                                                                                                                                  <  4:748603/148‑1 (MQ=38)

  gggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGAAAGGCACAGACCACCGAAGCGGGAATAGCCCAGTTGGTAGAGCACGACCTAGCCAAgg                                                                                                                                                >  2:877768/1‑149 (MQ=255)

     aaTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTTTGAAAGTGCTGTAAGACACAGACCACCCAAGCGGGAATATATTAGTTGTTATAGCACGAGCTAACCAAGGTc                                                                                                                                              >  4:158429/1‑148 (MQ=11)

      aTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                                                            <  2:465180/149‑1 (MQ=38)

       tAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                           <  2:435914/149‑1 (MQ=38)

       tAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                           <  2:758974/149‑1 (MQ=38)

         gCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGt                                                                                                                                         <  1:894801/149‑1 (MQ=38)

           tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTc                                                                                                                                        >  3:118122/1‑148 (MQ=38)

              gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                           >  2:604171/1‑142 (MQ=38)

              gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                           <  1:604171/142‑1 (MQ=40)

              gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgcg                                                                                                                                     >  4:348512/1‑148 (MQ=38)

                    agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGCAAGTGCTGTAAGGCACAGACCACCCAAGCAGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGTTCGGGGTCGCGAGtttg                                                                                                                               <  4:82180/148‑3 (MQ=255)

                    agagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                              >  1:418959/1‑149 (MQ=39)

                     gagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                              >  3:2703/1‑148 (MQ=39)

                      agCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                                                                                                                            >  2:851017/1‑149 (MQ=38)

                       gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAACGTGCTATAAGGAACACACCACCCCAGCGGGAATAGCTCAGTTGGTAGAGCACGTCCTTGCCAAGGTCGGGGTAGCGAGTTCGAgac                                                                                                                           >  2:867934/1‑147 (MQ=11)

                       gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGt                                                                                                                            <  3:642609/148‑1 (MQ=40)

                          cGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                                                         >  4:482957/1‑148 (MQ=38)

                           gACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCg                                                                                                                        >  4:340959/1‑148 (MQ=38)

                           gACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGt                                                                                                                       >  2:125729/1‑149 (MQ=38)

                            aCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGtt                                                                                                                      <  2:418959/149‑1 (MQ=38)

                              cTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTc                                                                                                                    >  1:357267/1‑149 (MQ=39)

                              cTTACCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGCAAGTGCTGTAAGTCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGTGTTCGAGTCTCGttt                                                                                                                     >  4:253689/1‑148 (MQ=255)

                                   cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCt                                                                                                               >  1:233488/1‑149 (MQ=38)

                                        tCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgt                                                                                                                                            >  4:782635/1‑114 (MQ=255)

                                        tCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGGCAAGGTCgt                                                                                                                                            <  3:782635/115‑2 (MQ=255)

                                        tCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaa                                                                                                           >  3:592419/1‑148 (MQ=40)

                                               cgcgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATtcttct                                                                                                   <  1:196750/149‑1 (MQ=255)

                                                gcgcgctCGAGTCTCGTTTCCCGTTCCAAAATTTTACAGTGCTGTAAGGCACAGACCCCCCTAGCGGAAAGAGCTCAGTAGGTAGACCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGCCTGGTGTCCCGCTCCAAATTCTtctc                                                                                                  <  1:961267/143‑1 (MQ=14)

                                                  gAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACTCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTtctctc                                                                                                <  2:323969/149‑1 (MQ=255)

                                                     ttCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAAt                                                                                             >  2:902437/1‑149 (MQ=255)

                                                        gAGTCTCGTTTCCCACTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaa                                                                                           >  4:782844/1‑148 (MQ=255)

                                                         aGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATaaa                                                                                          <  4:876202/148‑1 (MQ=255)

                                                                  tCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCcaca                                                                                 <  3:940669/148‑1 (MQ=255)

                                                                   cccGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAg                                                                                <  3:340959/148‑1 (MQ=255)

                                                                     cGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCg                                                                              >  3:508329/1‑148 (MQ=255)

                                                                      gCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                    >  3:577624/1‑125 (MQ=255)

                                                                      gCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                    <  4:577624/125‑1 (MQ=255)

                                                                      gCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttctt                                                                                                     <  3:209906/124‑1 (MQ=255)

                                                                      gCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttctt                                                                                                     >  4:209906/1‑124 (MQ=255)

                                                                         ccAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                    >  3:506919/1‑122 (MQ=255)

                                                                         ccAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttcttc                                                                                                    <  4:506919/122‑1 (MQ=255)

                                                                          cAAAATTGGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGAcgcg                                                                         <  4:188088/148‑1 (MQ=255)

                                                                           aaaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTcc                                                                                                             <  1:939188/111‑1 (MQ=40)

                                                                           aaaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTcc                                                                                                             >  2:939188/1‑111 (MQ=255)

                                                                            aaaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGAt                                                                       >  3:471220/1‑148 (MQ=255)

                                                                             aaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATg                                                                      <  3:650657/148‑1 (MQ=255)

                                                                             aaTTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATg                                                                      <  4:592419/148‑1 (MQ=255)

                                                                                 tGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGtt                                                                  <  3:482957/148‑1 (MQ=255)

                                                                                 tGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTa                                                                 <  2:642373/149‑1 (MQ=255)

                                                                                    aaGTGCTGTAAGGCACAGTCCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTg                                                              <  1:267256/149‑1 (MQ=255)

                                                                                    aaGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCg                                                                                                                               <  3:286547/84‑1 (MQ=38)

                                                                                           gTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGtt                                                        >  4:933965/1‑148 (MQ=255)

                                                                                            tAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCg                                                                                                                               <  4:160353/76‑1 (MQ=35)

                                                                                             aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                     <  1:428637/149‑1 (MQ=255)

                                                                                             aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                     <  1:822861/149‑1 (MQ=255)

                                                                                             aaGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                     <  1:901018/149‑1 (MQ=255)

                                                                                              aGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGttttt                                                     >  3:1013946/1‑148 (MQ=255)

                                                                                               ggCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTttg                                                    >  3:895613/1‑148 (MQ=255)

                                                                                                    aGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtc                                               >  3:417299/1‑148 (MQ=255)

                                                                                                     gACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGtctc                                             <  2:546694/149‑1 (MQ=255)

                                                                                                      accaccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTctct                                            <  1:226046/149‑1 (MQ=255)

                                                                                                        caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTg                                           <  3:463436/148‑1 (MQ=255)

                                                                                                        caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGa                                          >  1:1026379/1‑149 (MQ=255)

                                                                                                        caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGa                                          >  1:302203/1‑149 (MQ=255)

                                                                                                        caccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGa                                          >  1:418940/1‑149 (MQ=255)

                                                                                                            cAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaa                                       <  4:471220/148‑1 (MQ=255)

                                                                                                            cAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaa                                       >  4:236645/1‑148 (MQ=255)

                                                                                                              aGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAc                                     >  3:868174/1‑148 (MQ=255)

                                                                                                              aGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAc                                     >  3:873951/1‑148 (MQ=255)

                                                                                                                   gAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGt                               <  2:418940/149‑1 (MQ=255)

                                                                                                                        gCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAAc                           <  3:236645/148‑1 (MQ=255)

                                                                                                                         cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTg                                >  3:492955/1‑142 (MQ=255)

                                                                                                                         cTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTg                                <  4:492955/142‑1 (MQ=255)

                                                                                                                          tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACaga                        <  1:846505/149‑1 (MQ=255)

                                                                                                                          tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACaga                        <  1:541573/149‑1 (MQ=255)

                                                                                                                          tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACag                         <  3:782844/148‑1 (MQ=255)

                                                                                                                          tCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACag                         <  4:663719/148‑1 (MQ=255)

                                                                                                                           cAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACagag                       <  1:125729/149‑1 (MQ=255)

                                                                                                                           cAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACaga                        >  3:979273/1‑148 (MQ=255)

                                                                                                                           cAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACaga                        <  4:873951/148‑1 (MQ=255)

                                                                                                                           cAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACaga                        >  3:905819/1‑148 (MQ=255)

                                                                                                                                  tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATcc                 <  3:578150/148‑1 (MQ=255)

                                                                                                                                     agCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAg             <  1:877768/149‑1 (MQ=255)

                                                                                                                                     accccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAg             >  2:688256/5‑149 (MQ=255)

                                                                                                                                      gCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAttc                                                                                                       <  2:131132/58‑1 (MQ=35)

                                                                                                                                       cccGACCTTGGCAAGGTCGTGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaaa                                      <  1:355552/120‑1 (MQ=255)

                                                                                                                                       cccGACCTTGGCAAGGTCGTGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACaaa                                      >  2:355552/3‑122 (MQ=255)

                                                                                                                                       cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTg                                                           >  1:20920/3‑101 (MQ=255)

                                                                                                                                       cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTg                                                           <  2:594291/99‑1 (MQ=255)

                                                                                                                                       cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTg                                                           >  1:594291/3‑101 (MQ=255)

                                                                                                                                       cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTg                                                           <  2:20920/99‑1 (MQ=255)

                                                                                                                                        aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCt          >  1:494697/1‑149 (MQ=255)

                                                                                                                                        aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCt          >  1:79099/1‑149 (MQ=255)

                                                                                                                                         cGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAct                                    <  3:301829/122‑1 (MQ=255)

                                                                                                                                         cGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAAct                                    >  4:301829/1‑122 (MQ=255)

                                                                                                                                         cGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCTGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTc         >  2:197185/1‑149 (MQ=255)

                                                                                                                                              ttGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCAACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCt    >  1:687985/1‑149 (MQ=255)

                                                                                                                                                 ccAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTGAGGCTGt  <  4:905819/148‑1 (MQ=255)

GCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAATTCTTCTCTCAATAAAATATCCACAGCGACGCGATGCGTTATTGCTGGTTTTTGTTGTCTCTGACAAACTCTTGTAAACAGAGTTATCCACAGCCTCAGGCTGT  >  NC_000913/4392472‑4392764

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