BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,603,664	A→G	noncoding (71/88 nt)	leuV ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,603,664	0	A	G	97.9%	452.4 / 2.6	142	noncoding (71/88 nt)	leuV	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base A (0/1); major base G (100/37); minor base . (2/1); total (102/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.67e-01

Only 100 of 155 total aligned reads displayed.

ATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTT  >  NC_000913/4603574‑4603752

atcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCt                                                                                        >  1:5867387/1‑93 (MQ=255)

 tcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCtt                                                                                       >  2:5057114/1‑93 (MQ=255)

 tcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCtt                                                                                       <  2:5087621/93‑1 (MQ=255)

   atGCTGATGATAATGATTTTGAATCAGGGATGAAAAATCAAGGCAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGa                                                                                     >  2:2541402/1‑93 (MQ=255)

      cTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAAcc                                                                                  <  1:4362411/93‑1 (MQ=255)

       tgatgaTAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAAccc                                                                                 >  2:5459287/1‑93 (MQ=255)

       tgatgaTAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAAccc                                                                                 >  1:2433797/1‑93 (MQ=255)

       tgatgaTAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAAccc                                                                                 >  1:2316802/1‑93 (MQ=255)

          tgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCa                                                                              <  1:2743060/93‑1 (MQ=255)

           gatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAc                                                                             >  1:5617498/1‑93 (MQ=255)

            atAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACg                                                                            >  1:1019656/1‑93 (MQ=255)

             tAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACgc                                                                           >  2:2731970/1‑92 (MQ=255)

             tAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGt                                                                           <  1:1923696/93‑1 (MQ=255)

             tAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGt                                                                           >  2:6437995/1‑93 (MQ=255)

               aggAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTcc                                                                         >  1:4188445/3‑93 (MQ=255)

                tGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCg                                                                        >  1:133785/1‑93 (MQ=255)

                    tttGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg                                                                    >  1:5835895/1‑93 (MQ=255)

                     ttGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAgg                                                                   >  2:5166303/1‑93 (MQ=255)

                     ttGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCACCACGTCCGTAAGGa                                                                  <  2:5418228/93‑1 (MQ=255)

                     tggAATGAGGGAT‑AAAAAAGAAGGAAAAAAAAAAAAAGGGAAAAAAAGATAATTGGGGCGA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAgg                                                                   <  2:2699323/91‑1 (MQ=12)

                      tGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGa                                                                  >  1:2516450/1‑93 (MQ=255)

                         aTCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGacac                                                               >  1:1551035/1‑93 (MQ=255)

                         aTCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGacac                                                               >  1:3536329/1‑93 (MQ=255)

                         aTCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAAATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGacac                                                               <  2:6605919/93‑1 (MQ=255)

                            aGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTaa                                                            >  1:2139747/1‑93 (MQ=255)

                            aGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTaa                                                            >  2:1889984/1‑93 (MQ=255)

                            aGGGATAAAAAATCAAGAAAAAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTaa                                                            <  2:3020152/93‑1 (MQ=25)

                             gggATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAac                                                           >  1:712960/1‑93 (MQ=255)

                             gggATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAac                                                           >  2:2173736/1‑93 (MQ=255)

                                ataaaaaaaaaAGAAAGAAAAAGTAAA‑GGAATTAAGGATAATTGGTCCGA‑GGGGGGGGGTTTAACCCCCCACGTCCGTAAGGACACTAACAcc                                                        <  1:5263703/90‑1 (MQ=12)

                                aTGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACAcc                                                        <  2:2015078/93‑1 (MQ=255)

                                aTGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACAcc                                                        <  1:3130524/93‑1 (MQ=255)

                                aTGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACAcc                                                        <  1:1123419/93‑1 (MQ=255)

                                 tGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCt                                                       >  1:4187995/1‑93 (MQ=255)

                                       ataaaGAAAAAAACAAGAAA‑GGAAGAAAAGATAATTGGTGGGG‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCt                                                 <  2:575554/90‑1 (MQ=18)

                                        tCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTa                                                >  2:3893244/1‑93 (MQ=255)

                                        tCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAAGTCCGTAAGGACACTAAAACCTGAAGCTa                                                >  2:5570555/1‑93 (MQ=255)

                                         cAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                               >  2:4956273/1‑93 (MQ=255)

                                         cAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                               >  1:5322941/1‑93 (MQ=255)

                                         cAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                               >  1:413829/1‑93 (MQ=255)

                                         aaagaaagaAACAAGAAA‑GAAAGAAAAGATAGTGGGTGCGG‑GGGGGGGGGTTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                               <  2:5984365/92‑1 (MQ=17)

                                         aaagaaagaAAAAAAAAA‑GGAATAAAAGATAATTGGTGCGG‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                               <  1:2871743/92‑1 (MQ=17)

                                         aaaGGAAAAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                               <  2:4063941/92‑1 (MQ=21)

                                           aGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                             >  2:4170509/1‑93 (MQ=255)

                                           aGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                             >  2:6539171/1‑93 (MQ=255)

                                           aGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                             >  2:2231947/1‑93 (MQ=255)

                                           aGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                             >  2:1083189/1‑93 (MQ=255)

                                           aGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcgc                                            >  1:6439741/1‑93 (MQ=255)

                                             gaagaaACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                           >  2:6352809/1‑93 (MQ=255)

                                              aaaaaaCAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGt                                          <  2:6523509/90‑1 (MQ=255)

                                               agaaACAAGAAA‑GGAAGTAAAGATAATTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGt                                          >  1:3770626/1‑93 (MQ=255)

                                               agaaACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                                         <  2:4085147/93‑1 (MQ=255)

                                               agaaACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                                         >  2:2045864/1‑93 (MQ=255)

                                               aaaaaaaaaaaa‑GAAAGAAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                                         <  1:5044366/91‑1 (MQ=14)

                                                 aaaCAAGAAA‑GGAAGTAAAGATAATTGGCGCGG‑GGGGGGGGTCTTGAACCCCCACGTCCGTAAGGACACTAACCCCTGAAGCTAGCGCGTCTa                                       <  2:4472395/93‑1 (MQ=18)

                                                  aaCAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc                                      >  1:3915389/1‑93 (MQ=255)

                                                     aaGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa                                   >  1:5769909/1‑93 (MQ=255)

                                                       gAAA‑GGAAGTAAAGATAATTGGTGCTA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                 <  1:3295016/93‑1 (MQ=255)

                                                          a‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAAGACCTGAAGCTAGCGCGTCTACCAATTCCg                              >  2:6064874/1‑93 (MQ=255)

                                                            ggAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                             >  1:3213580/1‑93 (MQ=255)

                                                             gAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                            >  1:3617664/1‑93 (MQ=255)

                                                             gAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                            >  2:4221905/1‑93 (MQ=255)

                                                               aGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                          >  2:3848255/1‑93 (MQ=255)

                                                                 tAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                        >  1:2723824/1‑93 (MQ=255)

                                                                  aaaGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                       <  1:6464216/93‑1 (MQ=255)

                                                                  aaaGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                       >  1:1435420/1‑93 (MQ=255)

                                                                  aaaGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                       >  2:523235/1‑93 (MQ=255)

                                                                   aaGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                      >  1:2047660/1‑93 (MQ=255)

                                                                   aaGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                      >  1:4714620/1‑93 (MQ=255)

                                                                   aaGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGATAGCGCGTCTACCAATTCCGCCACCTTc                      >  1:2223332/1‑93 (MQ=255)

                                                                    aGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                     >  1:6260032/1‑93 (MQ=37)

                                                                     gATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                    >  1:2452419/1‑93 (MQ=25)

                                                                     gATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                    >  1:809074/1‑93 (MQ=25)

                                                                      aTAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGca                   >  1:3967232/1‑93 (MQ=21)

                                                                       tAATTGGTGCGA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcac                  <  2:657067/93‑1 (MQ=255)

                                                                       tAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcac                  >  1:1020378/1‑93 (MQ=21)

                                                                        aaTTGGTGGGA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcaca                 <  2:6212898/93‑1 (MQ=255)

                                                                        aaTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcaca                 >  1:4652078/1‑93 (MQ=21)

                                                                         aTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAg                >  2:2115008/1‑93 (MQ=21)

                                                                         aTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTATACCAATTCCGCCACCTTCGCACAg                >  1:720858/1‑93 (MQ=255)

                                                                          ttGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcac                  >  2:4283381/1‑90 (MQ=14)

                                                                           tGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcacacca              >  1:1607307/1‑90 (MQ=0)

                                                                           gggTGCGA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTc              <  2:6127587/92‑1 (MQ=255)

                                                                            ggTGGGA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCa             <  2:5184435/93‑1 (MQ=255)

                                                                              ggggA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATc           <  2:1992866/90‑1 (MQ=255)

                                                                               gCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCt          <  2:890898/93‑1 (MQ=255)

                                                                               gCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCt          <  1:3099888/93‑1 (MQ=255)

                                                                                ggggggggggggGTTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCt          <  1:2638737/89‑1 (MQ=14)

                                                                                 gggggggggggGCTGAACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCtt         <  1:4920522/90‑1 (MQ=14)

                                                                                 gA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTa        <  2:689777/93‑1 (MQ=255)

                                                                                 gA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTa        >  2:5136788/1‑93 (MQ=255)

                                                                                   gggggggggCCTCCACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTAc       <  1:4577394/92‑1 (MQ=14)

                                                                                    ggggtgggGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACt      <  2:1094164/93‑1 (MQ=255)

                                                                                    ggggggggGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACt      >  2:4502406/1‑93 (MQ=255)

                                                                                    ggggggggGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACACTCATCTTACt      >  2:4969407/1‑93 (MQ=255)

                                                                                     ggggggggCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCCTCTTACtt     >  2:3165606/1‑93 (MQ=255)

                                                                                      ggggggGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttt    <  1:1143812/93‑1 (MQ=255)

                                                                                       ggggggCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACtttt   <  2:698504/93‑1 (MQ=255)

                                                                                        ggggGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttttt  <  2:1845860/93‑1 (MQ=255)

                                                                                        ggggGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttttt  >  1:3404350/1‑93 (MQ=255)

ATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAA‑GGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTT  >  NC_000913/4603574‑4603752

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`