BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,603,785	0	A	G	96.8%	822.1 / 18.5	285	noncoding (71/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (158/117); minor base . (6/3); total (164/121)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

Only 100 of 319 total aligned reads displayed.

CACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAA  >  NC_000913/4603695‑4603874

  ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTg                                                                                       >  1:8316124/1‑93 (MQ=255)

   cTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGa                                                                                      <  2:3266873/93‑1 (MQ=255)

    tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGGCTCGTCTGGCGC‑GA‑‑‑GG‑GGGGGGGCTTGaa                                                                                     >  2:7520504/1‑93 (MQ=37)

      aaGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAAcc                                                                                   <  1:867394/93‑1 (MQ=255)

      aaGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAAcc                                                                                   <  1:9386136/93‑1 (MQ=255)

       aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCTTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAAccc                                                                                  <  2:5854613/93‑1 (MQ=255)

          tAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGGACCCCCa                                                                               >  2:513152/1‑93 (MQ=255)

          tAGCGCCTCTACCAATTCCGCCACCTTCGCAAAGTCATCTTCCTTTTTTTGATATCGCCACGTTTGGGGC‑GAG‑‑GG‑GGGGGGGCTTGAAccccc                                                                                <  2:3085158/93‑1 (MQ=25)

          tAGCCCTTCTACCAATTCCGCCTTCTTCGCTCAGTTATCTTACTTTTTTTGATATCGCCTCGTTGGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAAccccc                                                                                <  2:10037984/93‑1 (MQ=255)

           aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCa                                                                               <  2:2080343/93‑1 (MQ=38)

           aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCa                                                                               <  1:7781177/93‑1 (MQ=38)

            gcgcGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCAc                                                                              <  2:1544245/93‑1 (MQ=255)

            gcgcGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACg                                                                             >  1:6507477/1‑93 (MQ=255)

              gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAGG‑GG‑GGGGGGGCTTGAACCCCCAcc                                                                             >  1:7426915/1‑92 (MQ=38)

              gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTc                                                                           >  1:6288005/1‑93 (MQ=255)

               cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTc                                                                           >  1:5955160/1‑93 (MQ=255)

               cgTCTACCAATCCCGCCCCCTTCGCACAGTCATCTTATTTTTTTTGATATCCCCTCTTTTGGTGC‑GAGG‑GG‑GGGGGGGCTTGAACCCCCACGt                                                                            <  2:596696/93‑1 (MQ=25)

                   tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTa                                                                       >  2:7020763/1‑93 (MQ=255)

                    aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTaa                                                                      >  2:8274189/1‑93 (MQ=255)

                    aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTaa                                                                      >  2:22210/1‑93 (MQ=255)

                    aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAg                                                                     <  1:6018160/93‑1 (MQ=255)

                    aCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAg                                                                     >  1:4802323/1‑93 (MQ=255)

                     ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGa                                                                    >  1:9091618/1‑93 (MQ=255)

                      cAACTCCCCCCCCTTCCCACAGTCATTTTATTTTTTTTGATTTCGCCTCGTTTGGGGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGa                                                                    <  2:7914546/93‑1 (MQ=21)

                        aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTTC‑GA‑‑‑GGAGGGGGGGGTTGAACCCCGAGGTCCGAAAAAac                                                                  >  1:8056189/1‑93 (MQ=255)

                        aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAaca                                                                 >  2:2593291/1‑93 (MQ=255)

                        aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAaca                                                                 <  1:5704581/93‑1 (MQ=255)

                          tCCGCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAacac                                                                <  1:9254447/93‑1 (MQ=37)

                          tCCGCCCCCCTCGCACAGGCATTTTTCTTTTTTTGATATCCCCTCGTTTGGGGCGGG‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAacac                                                                <  1:1177257/93‑1 (MQ=21)

                             gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGG‑CTTGAACCCCCACGTCCGTAAGAACACTAaca                                                           >  1:704566/1‑93 (MQ=255)

                             cccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAGG‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTa                                                              <  1:1141803/92‑1 (MQ=37)

                              ccACCTTCGCACAGTCATCTTCCTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAac                                                            <  1:9165113/93‑1 (MQ=37)

                               cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAacac                                                          <  1:2329099/93‑1 (MQ=255)

                                 ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCCCGTCCGTAAGAACACTAACACCt                                                        >  2:7225319/1‑93 (MQ=255)

                                 ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                        >  1:4513795/1‑93 (MQ=255)

                                 ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                        >  2:6405509/1‑93 (MQ=255)

                                 ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGG‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTg                                                       >  2:4536469/1‑93 (MQ=255)

                                 ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGG‑CTTGAACCCCCACGTCCGTAAGAACACTAACACATg                                                       >  2:5244745/1‑93 (MQ=255)

                                  cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAGGGGG‑GGGGGGGCTTGAACCCCCCCCGCCGAAAaaacacaaacac                                                          >  2:6600216/1‑93 (MQ=255)

                                  cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTg                                                       >  1:9487792/1‑93 (MQ=255)

                                    tCGCACAGTCATCTTACTTTTATTGACATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGTGCCTGAACTCCCACGTCCGTAAGAACATCTACACCTGaa                                                     >  2:8958301/1‑93 (MQ=21)

                                     cGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAg                                                    <  2:8019170/93‑1 (MQ=255)

                                      gCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAg                                                    <  1:3427799/93‑1 (MQ=255)

                                      gCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGc                                                   <  1:7517645/93‑1 (MQ=255)

                                        acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                 >  1:2010554/1‑93 (MQ=255)

                                          aGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                               >  2:4289944/1‑93 (MQ=255)

                                           gTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAAACCCCACGTCCGTAAGAACACTAACACCTGACGCTAgcg                                              >  1:9426424/1‑93 (MQ=255)

                                               tCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                          >  2:1774370/1‑93 (MQ=255)

                                                cTTACTTTTTTTGATATCGCCTCGTTTGGTGCGGG‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                          <  1:6745143/93‑1 (MQ=21)

                                                cTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCCCGTCCGTAAGAAACCTAACACCTGAAGCTACGGCGTc                                          >  2:9438005/1‑93 (MQ=18)

                                                    cTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                     >  1:2803551/1‑93 (MQ=255)

                                                     tttttttGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCGCGTCCGTAAGAACACTAACACCTGAACCTAGCGCGTCTCCCa                                     >  2:5552718/1‑93 (MQ=18)

                                                     tttttttGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                     <  1:4367565/93‑1 (MQ=255)

                                                     tttttttGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                    >  2:1300390/1‑93 (MQ=255)

                                                     tttgtttgATATCGCCTGGTTGGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                    <  2:77640/93‑1 (MQ=17)

                                                      ttttttGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGG‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                  >  1:3184185/1‑93 (MQ=255)

                                                       tttttGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTTAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCTTCTACCAAt                                   >  1:6940954/1‑93 (MQ=17)

                                                         tttGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                 <  1:8123162/93‑1 (MQ=255)

                                                          ttGATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                >  2:5060981/1‑93 (MQ=255)

                                                            gATATCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAAAACCTGAAGCTAGCGCGTCTACCAATTCCGc                              >  2:3122357/1‑93 (MQ=17)

                                                             atatCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                             >  1:3027734/1‑93 (MQ=255)

                                                              tatCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGCCTACCAATTCCGccc                            >  1:3135218/1‑92 (MQ=14)

                                                              tatCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAAACCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAATCCGCCa                            >  1:4600011/1‑93 (MQ=14)

                                                               ttCGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                           <  1:4111930/92‑1 (MQ=14)

                                                                tCGCCTCGTTTGGTGC‑GAGG‑GG‑GGGGGGGCTTGAACCCCCCCGTCCGTAAGAACACTAACACCTGAAGCTAGCGAGTCTACCAATTCCCCCAc                           >  2:3492928/1‑93 (MQ=14)

                                                                tCGCCTCGTTTGGTGC‑GAGG‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACAATAACACCTGAAGCTAGCGCGGCTACCAATTCCGCCAc                           >  1:5604897/1‑93 (MQ=14)

                                                                tCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGG‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                        >  2:5745905/1‑93 (MQ=255)

                                                                tCGCCTCGTTTGGGGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                          <  2:4771833/93‑1 (MQ=14)

                                                                 cGCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                         >  1:9730170/1‑93 (MQ=255)

                                                                  gCCTCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                        >  1:99313/1‑93 (MQ=255)

                                                                  gCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                       <  1:5328783/93‑1 (MQ=255)

                                                                  gCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                       <  1:6214847/93‑1 (MQ=255)

                                                                     tCGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                     >  1:1457410/1‑93 (MQ=255)

                                                                      cGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                    >  2:8756650/1‑93 (MQ=255)

                                                                      cGTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                    >  1:3758272/1‑93 (MQ=255)

                                                                       gTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                   >  2:3812217/1‑93 (MQ=255)

                                                                       gTTTGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                   >  1:7605792/1‑93 (MQ=255)

                                                                       gTTTGGTGC‑GA‑‑‑GG‑GGGGGG‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg                 >  2:986092/1‑93 (MQ=255)

                                                                        tttGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                  <  2:1573278/93‑1 (MQ=255)

                                                                        tgtGGTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                  <  1:287590/91‑1 (MQ=11)

                                                                          tGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg               <  1:5045531/93‑1 (MQ=255)

                                                                          tGGTGC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg               <  1:2488989/93‑1 (MQ=255)

                                                                            gTGC‑GAG‑‑GG‑GGGGGGGCTTGAACCCCCACGCCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc              >  2:2910434/1‑93 (MQ=14)

                                                                              gC‑GA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg           <  1:8564670/93‑1 (MQ=255)

                                                                              gC‑GA‑‑‑GG‑GGGGGG‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGc          <  1:5770547/93‑1 (MQ=255)

                                                                                 gA‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAAAACTAACACCTGAAGCTAGCGCGTCTTCCAATTCCGCCACCTTCGCCCACTGCggaca         >  2:2322516/1‑89 (MQ=14)

                                                                                  gagggg‑ggggggCCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg           <  2:5261617/89‑1 (MQ=255)

                                                                                  a‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                              >  1:3565123/1‑55 (MQ=255)

                                                                                  a‑‑‑GG‑GGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCaa        <  1:6672232/93‑1 (MQ=255)

                                                                                     ggg‑ggggggGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCaa        <  1:4645063/92‑1 (MQ=255)

                                                                                     ggg‑ggggggGCTTGAACCCCAACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCaa        <  1:1865332/92‑1 (MQ=255)

                                                                                      gg‑ggggggGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                              >  2:6852779/1‑70 (MQ=25)

                                                                                      gg‑ggggggGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat       >  2:610000/1‑93 (MQ=255)

                                                                                       g‑gggggggCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCCCCTTCGCCCAGTGCGAGCAata      <  2:8834508/93‑1 (MQ=255)

                                                                                       g‑gggggggCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata      <  2:2825511/93‑1 (MQ=255)

                                                                                         ggggggGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatat     >  1:3863490/1‑93 (MQ=255)

                                                                                         ggggggGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatat     <  2:4901219/93‑1 (MQ=255)

                                                                                           ggggGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCa   >  2:1683632/1‑93 (MQ=255)

                                                                                           ggggGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCa   <  2:6744535/93‑1 (MQ=255)

                                                                                           ggggGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCa   <  1:9850838/93‑1 (MQ=255)

CACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAA  >  NC_000913/4603695‑4603874

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