D17 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007790	918	A→G	T213T (ACA→ACG)	SAUSA300_RS14705 →	replication protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007790	918	0	A	G	90.6%	471.1 / 23.9	215	T213T (ACA→ACG)	SAUSA300_RS14705	replication protein
Reads supporting (aligned to +/- strand): ref base A (15/5); new base G (184/10); total (200/15)
Fisher's exact test for biased strand distribution p-value = 7.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-08

Only 100 of 221 total aligned reads displayed.

TAAGTTATGGTGGATTGCTTAAACAAAAACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAAC‑AACAGACGAAGAAGAAAAAGCACATTCAATTACCGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAGACATTAGTGTTGAC  >  NC_007790/821‑1003

gcccttATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAaac‑gac                                                                                     >  1:3306926/5‑100 (MQ=255)

  aGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAg                                                                                   >  1:1886538/1‑100 (MQ=255)

    ttATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGAc                                                                                 <  1:4086432/100‑1 (MQ=255)

     tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAATGC‑GACAGatc                                                                                >  1:4142086/1‑98 (MQ=255)

     tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAGC‑GACAGatc                                                                                >  1:5361146/1‑98 (MQ=255)

     tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGAcc                                                                                >  1:3170368/1‑99 (MQ=255)

     tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGAca                                                                                >  1:3598581/1‑99 (MQ=255)

     tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACg                                                                                >  1:6065205/1‑100 (MQ=255)

     aaTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑AACAGaga                                                                                >  1:2335753/2‑98 (MQ=255)

      aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAATGC‑GACAGatcg                                                                               >  1:3510378/1‑97 (MQ=255)

      aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAac‑gacaatcga                                                                               >  1:2367318/1‑95 (MQ=255)

      aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAac‑aacaaagat                                                                               >  1:1035788/1‑95 (MQ=255)

      aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGAACAGatc                                                                                >  1:5842444/1‑98 (MQ=255)

      aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGatcg                                                                               >  1:1066781/1‑97 (MQ=255)

      aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACga                                                                               >  1:3045801/1‑100 (MQ=255)

      aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACga                                                                               >  1:3936170/1‑100 (MQ=255)

      aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAAACga                                                                               <  1:4457173/100‑1 (MQ=255)

      aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑AACAGatcg                                                                               >  1:5235762/1‑97 (MQ=255)

      aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑AACAGatcg                                                                               >  1:5818753/1‑97 (MQ=255)

        tgtggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACgaag                                                                             >  1:4024419/2‑100 (MQ=255)

        tgtggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑AACAGatcgga                                                                             >  1:5074517/2‑95 (MQ=255)

        cgtggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAGC‑GACAGatcgga                                                                             >  1:1643209/2‑95 (MQ=255)

         atggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGatcggaa                                                                            >  1:5168813/2‑94 (MQ=255)

           ggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACgaagaag                                                                          >  1:4572717/1‑100 (MQ=255)

           ggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GAAAGatcggaaga                                                                          >  1:5443939/1‑92 (MQ=255)

                      gccgtggcATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACAt                                                               >  1:4707198/8‑100 (MQ=255)

                              aTAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAGGAAAAAGCACATTCGATTAc                                                       <  1:5188321/100‑1 (MQ=255)

                              aTAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTAc                                                       <  1:4898876/100‑1 (MQ=255)

                                       tAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAAtaag                                              >  1:692560/1‑97 (MQ=255)

                                       tAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTg                                              >  1:1263688/1‑100 (MQ=255)

                                        aaaCTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTgg                                             >  1:3442104/1‑100 (MQ=255)

                                        aaaCTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTgg                                             >  1:2795202/1‑100 (MQ=255)

                                        aaaCTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTgg                                             >  1:1610953/1‑100 (MQ=255)

                                         aatttAGATAATGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAAC‑AACAGACGAAGAGGAAAAAGCACATTCAATTACAGCAATCTGGa                                            <  1:66264/97‑1 (MQ=255)

                                           cctAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAAt                                          >  1:2430823/3‑100 (MQ=255)

                                            atAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAGC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAAtt                                         >  1:362213/2‑100 (MQ=255)

                                             tAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATCTGGaaaga                                        >  1:2577506/1‑97 (MQ=255)

                                              aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAG‑GACAGACGAAGAAGAAAAAGCACATTCGATGACAGCAAtagatcggaag                                       >  1:1956062/1‑90 (MQ=255)

                                              aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGCAAAAGCACATTCGAGTACAGCAATTTGGaaagat                                       >  1:4810170/1‑96 (MQ=255)

                                              aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                       >  1:3794685/1‑100 (MQ=255)

                                              aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                       >  1:3907861/1‑100 (MQ=255)

                                              aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                       >  1:5020194/1‑100 (MQ=255)

                                              aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                       >  1:2803879/1‑100 (MQ=255)

                                              aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                       >  1:2522490/1‑100 (MQ=255)

                                              aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                       >  1:1596637/1‑100 (MQ=255)

                                              aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACTTGGaagatc                                       >  1:3317063/1‑96 (MQ=255)

                                              aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAACAAttagatcggaa                                       >  1:5975722/1‑91 (MQ=255)

                                                aCGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGa                                     >  1:1771274/1‑100 (MQ=255)

                                                aCGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGa                                     >  1:5267098/1‑100 (MQ=255)

                                                 tgATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGaagatcgga                                    >  1:4307438/2‑93 (MQ=255)

                                                 tgATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTAGGaagatcgga                                    >  1:5410427/2‑93 (MQ=255)

                                                 tcaTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                    >  1:3268443/3‑100 (MQ=255)

                                                 taaTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                    >  1:2633697/3‑100 (MQ=255)

                                                 gtgtgccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                    >  1:5744451/7‑100 (MQ=255)

                                                 caaTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                    >  1:6006805/3‑100 (MQ=255)

                                                 cGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                    >  1:4646701/1‑100 (MQ=255)

                                                 cGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                    >  1:154138/1‑100 (MQ=255)

                                                 acaTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACATGGaagatcgga                                    >  1:859209/3‑93 (MQ=255)

                                                  gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                   >  1:1193275/1‑100 (MQ=255)

                                                  gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                   >  1:3999545/1‑100 (MQ=255)

                                                  gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                   >  1:5513056/1‑100 (MQ=255)

                                                  gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                   >  1:3815417/1‑100 (MQ=255)

                                                  gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTGGGagatcggaag                                   >  1:812974/1‑91 (MQ=255)

                                                  gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTAGGaagatcggaa                                   >  1:3538761/1‑92 (MQ=255)

                                                  gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATCTGGaaagatcgga                                   >  1:1399803/1‑92 (MQ=255)

                                                  gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACTTGGaagatcggaa                                   >  1:2688305/1‑92 (MQ=255)

                                                  gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGACTACAGCAATTTGGAATTTCGaaa                                   >  1:2831110/1‑100 (MQ=255)

                                                  aaTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                   >  1:134982/2‑100 (MQ=255)

                                                  aaTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                   >  1:289060/2‑100 (MQ=255)

                                                   gtgccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                  >  1:4205688/5‑100 (MQ=255)

                                                   gtgccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                  >  1:5798203/5‑100 (MQ=255)

                                                   aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                  >  1:1159522/1‑100 (MQ=255)

                                                   aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                  >  1:5688862/1‑100 (MQ=255)

                                                    tgccGAAGATGGGAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAg                                 >  1:4362155/4‑100 (MQ=255)

                                                    agccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAg                                 >  1:3101964/4‑100 (MQ=255)

                                                     gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGGAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                                >  1:498090/3‑100 (MQ=255)

                                                     gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                                >  1:4508955/3‑100 (MQ=255)

                                                     gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                                >  1:4247670/3‑100 (MQ=255)

                                                     gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                                >  1:387108/3‑100 (MQ=255)

                                                     gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                                >  1:2980536/3‑100 (MQ=255)

                                                     gccGAAGATGGTAATTTGATTAATACAAGCGATGAAGATAAAAC‑GACAGACGAAGAAGCAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                                >  1:2268017/3‑100 (MQ=255)

                                                      ccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCa                               >  1:2498475/2‑100 (MQ=255)

                                                          aGATGGTAATTTGATTAGTACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAat                           >  1:5385183/1‑100 (MQ=255)

                                                           gATGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAatt                          >  1:739904/1‑100 (MQ=255)

                                                           cctGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAatt                          >  1:2216283/3‑100 (MQ=255)

                                                            ttGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAatta                         >  1:4933759/2‑100 (MQ=255)

                                                            gtGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAatta                         >  1:1476982/2‑100 (MQ=255)

                                                            aTGGTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAatta                         >  1:1522638/1‑100 (MQ=255)

                                                              ggTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattatt                       >  1:4534500/1‑100 (MQ=255)

                                                              ggTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattatt                       >  1:3130811/1‑100 (MQ=255)

                                                              ggTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattatt                       >  1:4028988/1‑100 (MQ=255)

                                                              ggTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattatt                       >  1:16931/1‑100 (MQ=255)

                                                              ggTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattatt                       >  1:2094288/1‑100 (MQ=255)

                                                              ggTAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattatt                       >  1:4731078/1‑100 (MQ=255)

                                                              gatAATTTGATTAATACAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattatt                       >  1:2919336/3‑100 (MQ=255)

                                                                          cgatctaGTGACGAAGATAAAAC‑AACAGACGAAGAGGAAAAAGCACATTCAATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAGACAtt           <  1:1981464/94‑1 (MQ=255)

                                                                            gccAAGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTaaaaaatcta         >  1:4185538/3‑92 (MQ=255)

                                                                               aaGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAAAAAtttaagg      >  1:4490604/1‑95 (MQ=255)

                                                                                aGTGATGAAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAAAAAttcaagga     >  1:666961/1‑94 (MQ=255)

                                                                                   gatgaAGATAAAAC‑GACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAAAAAtttaaggagca  <  1:3203668/100‑10 (MQ=255)

TAAGTTATGGTGGATTGCTTAAACAAAAACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAAC‑AACAGACGAAGAAGAAAAAGCACATTCAATTACCGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAGACATTAGTGTTGAC  >  NC_007790/821‑1003

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