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	87.4%	388.5 / 57.0	145	noncoding (71/88 nt)	leuV	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base A (0/2); major base G (89/36); minor base . (15/3); total (104/41)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-01

Only 100 of 154 total aligned reads displayed.

ATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTT  >  NC_000913/4603574‑4603753

atcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCt                                                                                         >  1:1517181/1‑93 (MQ=255)

atcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCt                                                                                         >  2:5066902/1‑93 (MQ=255)

 tcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCCAGG‑AAGAAAAAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGAGGGGGCtt                                                                                        >  1:5163893/1‑93 (MQ=255)

 tcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCtt                                                                                        >  2:6099224/1‑93 (MQ=255)

 tcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCtt                                                                                        >  1:5739766/1‑93 (MQ=255)

 tcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCtg                                                                                        >  1:5271760/1‑92 (MQ=255)

 tcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCtg                                                                                        >  1:5163892/1‑92 (MQ=255)

     gCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCttcagc                                                                                    >  1:3395753/1‑89 (MQ=255)

       tgatgaTAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAAccc                                                                                  >  2:5307323/1‑93 (MQ=255)

           gatAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAc                                                                              >  2:2407461/1‑93 (MQ=255)

           gatAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAc                                                                              >  2:4482136/1‑93 (MQ=255)

            atAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACg                                                                             >  1:4689014/1‑93 (MQ=255)

            atAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACg                                                                             >  1:3807869/1‑93 (MQ=255)

             tAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGt                                                                            >  1:5858955/1‑93 (MQ=255)

             tAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGt                                                                            >  2:2807243/1‑93 (MQ=255)

              aaTGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTc                                                                           <  1:3948664/93‑1 (MQ=255)

               aTGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCg                                                                         >  1:5565503/1‑93 (MQ=255)

                   aTTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTaa                                                                      >  1:2541247/1‑93 (MQ=255)

                    tttGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg                                                                     <  1:6126873/93‑1 (MQ=255)

                    tttGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg                                                                     <  1:6126872/93‑1 (MQ=255)

                    tttGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg                                                                     >  2:1057601/1‑93 (MQ=255)

                    tttGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAgg                                                                    >  1:5836006/1‑93 (MQ=255)

                       gAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGac                                                                  >  2:3061501/1‑93 (MQ=255)

                       gAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGac                                                                  <  2:2348963/93‑1 (MQ=255)

                       gAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGac                                                                  >  2:4649449/1‑93 (MQ=255)

                        aaTCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGaca                                                                 >  1:5002301/1‑93 (MQ=255)

                         aTCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGacac                                                                >  1:3789313/1‑93 (MQ=255)

                         aTCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGacac                                                                >  1:880827/1‑93 (MQ=255)

                          tCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACt                                                               <  2:2280188/93‑1 (MQ=255)

                           aaaggAAGAAAAAAAAAGGAAAAAAACAAGAAAGGAAGAAAAGATATTGGGGGGGA‑GGGGGGGGGTTTGACCCCCCACGTCCGTAAGGACACt                                                               <  2:1515655/90‑1 (MQ=12)

                            aGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTaa                                                             >  1:1167961/1‑93 (MQ=255)

                             gggATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAac                                                            >  1:5688961/1‑93 (MQ=255)

                              ggATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAaca                                                           >  1:2197044/1‑93 (MQ=255)

                               gATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACAcc                                                         >  2:2152070/1‑93 (MQ=255)

                                 tGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCt                                                        >  2:2568550/1‑93 (MQ=255)

                                 tGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCt                                                        <  2:3304055/93‑1 (MQ=255)

                                 tGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGAGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACCAACACCt                                                        >  1:4971260/1‑93 (MQ=255)

                                 aaaaaaaTAAAGA‑AAGAAAAAAAAAAGGAAGTAAAGATAATTGGGGCGA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCt                                                        <  2:2670994/91‑1 (MQ=18)

                                  gAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTg                                                       >  1:6119353/1‑93 (MQ=255)

                                  gAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTg                                                       >  2:4175276/1‑93 (MQ=255)

                                     aaaTGAAAG‑AAAAAAAAAGAAAAAAAGTAAAAATAATTGGTGCGG‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAg                                                    <  2:2886780/93‑1 (MQ=17)

                                     aaaTCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAg                                                    >  1:5314007/1‑93 (MQ=255)

                                      aaTCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATATTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGc                                                   <  2:4300933/93‑1 (MQ=255)

                                        tCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                >  2:1866375/1‑93 (MQ=255)

                                        aaaaGA‑AAAAAAAAAGAAAGGAAGTAAAGATAATTGGGGGGG‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTa                                                 <  2:1488424/91‑1 (MQ=18)

                                          aaGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgc                                               >  2:4955499/1‑93 (MQ=255)

                                           aGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGGAAACTAACACCTGAAGCTAgc                                               >  2:1688688/1‑93 (MQ=21)

                                           aGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                              <  1:5231019/93‑1 (MQ=255)

                                           aGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                              >  1:3977522/1‑93 (MQ=255)

                                           aGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcgc                                             >  2:1031579/1‑93 (MQ=255)

                                            gg‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcgc                                             >  1:6108195/1‑93 (MQ=255)

                                             g‑aagaaACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                            >  2:4404618/1‑93 (MQ=255)

                                               aagaaACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGt                                           >  2:4352800/1‑93 (MQ=255)

                                               aagaaACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGt                                           >  1:4238712/1‑93 (MQ=255)

                                                agaaACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                                          >  1:1371404/1‑93 (MQ=255)

                                                agaaACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                                          <  1:4373801/93‑1 (MQ=255)

                                                 gaaACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGAGCGTCTa                                        >  2:3728429/1‑93 (MQ=255)

                                                  aaaCAAGAAAGGAAGTAAATATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc                                       <  1:2025260/93‑1 (MQ=255)

                                                   aaCAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc                                       <  2:5322466/93‑1 (MQ=255)

                                                    aCAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAcc                                      >  2:3560368/1‑93 (MQ=255)

                                                     cAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTGTACCa                                     >  1:2746542/1‑93 (MQ=255)

                                                      aaGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa                                    <  2:3012165/93‑1 (MQ=255)

                                                      aaGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                   >  1:927283/1‑93 (MQ=255)

                                                       aGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                   <  1:3248956/93‑1 (MQ=255)

                                                       aGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                  >  1:2824960/1‑93 (MQ=255)

                                                       aGAAAGGAAGCAAAGATAAGTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGCCCGTAAGGACACGAACACCTGAAGCGAGCGCGTCTACCAAt                                   >  2:5493889/1‑93 (MQ=14)

                                                        gAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                  <  1:3335567/93‑1 (MQ=255)

                                                        gAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                 >  2:5919170/1‑93 (MQ=255)

                                                        aaaagaaagtaaaga‑aagtGGGGGGG‑GGGGGGGGGCTGGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                 <  2:1673062/92‑1 (MQ=11)

                                                        aaaaGGAAGAAAAGAAATTTGGTGCG‑‑GGGGGGGGGTTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                 <  1:5090911/92‑1 (MQ=14)

                                                         aaaGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                 >  1:1513428/1‑93 (MQ=255)

                                                          aaGGAAGTAGAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                >  1:4186060/1‑93 (MQ=255)

                                                          aaGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                >  2:5361671/1‑93 (MQ=255)

                                                            ggAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                              <  1:5381203/93‑1 (MQ=255)

                                                            ggAAGAAAAGATAATTGGGGCGA‑GGGGGGGGGTTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                              <  1:3046070/93‑1 (MQ=14)

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

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

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

                                                                 tAAAGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGTCACCt                         >  1:3064194/1‑93 (MQ=255)

                                                                  aaaGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                        <  2:2214168/93‑1 (MQ=39)

                                                                   aaGATAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCCCCACCTTc                       >  1:82340/1‑93 (MQ=255)

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

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

                                                                    aGATAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                     >  2:5897073/1‑93 (MQ=255)

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

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

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

                                                                       tAATTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcat                   >  1:2811177/1‑92 (MQ=17)

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

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

                                                                       tAATTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcaca                  >  2:4246861/1‑93 (MQ=255)

                                                                        aaTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATCCCGCCACCTTCGcacg                  >  1:2170971/1‑92 (MQ=255)

                                                                        aaTGGGTGCGG‑GGGGGGGGGCTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcaca                  <  1:5024280/93‑1 (MQ=11)

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

                                                                               gCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCtt          >  2:3750538/1‑93 (MQ=255)

                                                                                ggA‑GGGGGAGGACTCCCCCCCCCACGACCGTAAGGACACTAACCCCTGAAGCTAGCGCGTCTACCATATCCGCCACCTTCGCACAGTCATCtt          <  1:5488191/92‑1 (MQ=14)

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

                                                                                  a‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATcattc        >  1:1913433/1‑89 (MQ=255)

                                                                                   gggggggggGCTTGCCCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTAc        <  2:5639331/92‑1 (MQ=14)

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

ATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGG‑AAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTT  >  NC_000913/4603574‑4603753

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