D17 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007790	945	A→G	S222S (TCA→TCG)	SAUSA300_RS14705 →	replication protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007790	945	0	A	G	95.0%	428.5 / 4.7	160	S222S (TCA→TCG)	SAUSA300_RS14705	replication protein
Reads supporting (aligned to +/- strand): ref base A (3/5); new base G (146/6); total (149/11)
Fisher's exact test for biased strand distribution p-value = 2.86e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 161 total aligned reads displayed.

AAAACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAACAACAGACGAAGAAGAAAAAGCACATTCAATTACCGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAGACATTAGTGTTGAC  >  NC_007790/846‑1003

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

        aaaTTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGc                                                    <  1:67468/100‑1 (MQ=255)

        aaaTTTTAAACCTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGc                                                    <  1:3032440/100‑1 (MQ=255)

             tgcaaCTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAAttt                                               >  1:3420370/4‑100 (MQ=255)

             gtAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATAACAGCagatc                                               >  1:4357691/2‑96 (MQ=255)

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

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

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

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

                aatttAGATAATGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAACAACAGACGAAGAGGAAAAAGCACATTCAATTACAGCAATTGGGa                                            <  1:5403151/97‑1 (MQ=255)

                aaCTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACCTGGa                                            <  1:4016796/100‑1 (MQ=255)

                  tctAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACTTGGaag                                          >  1:3570474/3‑99 (MQ=255)

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

                   ctAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAAtt                                         >  1:4243546/2‑100 (MQ=255)

                    tcgACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAAttt                                        >  1:2959944/3‑100 (MQ=255)

                    tagataaTGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAACAACAGACGAAGAGGAAAAAGCACATTCAATTACAGCAATTTGGagaga                                        >  1:4256768/7‑96 (MQ=255)

                    tagataaTGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAACAACAGACGAAGAGGAAAAAGCACATTCAATTACAGCAATTTGGaaaga                                        >  1:854811/7‑97 (MQ=255)

                    tAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTCGGaaaga                                        >  1:3867342/1‑97 (MQ=255)

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

                    tAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATATGGaaaga                                        >  1:1696359/1‑97 (MQ=255)

                    tAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACTTGGaagat                                        >  1:4929895/1‑97 (MQ=255)

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

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

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

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

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

                     aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                       >  1:3955533/1‑100 (MQ=255)

                     aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                       >  1:1092286/1‑100 (MQ=255)

                     aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATGTGGaagatc                                       >  1:5246746/1‑96 (MQ=255)

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

                     aGACGATGCCGAAGATGGTAATTTGATTAATACAAGAGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGaagatc                                       >  1:6006033/1‑96 (MQ=255)

                      gACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCg                                      >  1:1419696/1‑100 (MQ=255)

                       tgtgtgccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGa                                     >  1:2807889/8‑100 (MQ=255)

                       ccGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGa                                     >  1:460131/2‑100 (MQ=255)

                       aCGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGaaagatcg                                     >  1:4368791/1‑94 (MQ=255)

                       aCGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGa                                     >  1:2518331/1‑100 (MQ=255)

                       aCGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGa                                     >  1:2903390/1‑100 (MQ=255)

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

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

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

                        tgATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATCTGGaagatcgga                                    >  1:665972/2‑93 (MQ=255)

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

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

                        ggATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                    >  1:2800431/2‑100 (MQ=255)

                        ctaTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATCTGGaaagatcgg                                    >  1:924747/3‑93 (MQ=255)

                        cgacgccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                    >  1:4501870/7‑100 (MQ=255)

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

                        cGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                    >  1:2790218/1‑100 (MQ=255)

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

                         taTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGatagatcgga                                   >  1:4178218/2‑91 (MQ=255)

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

                         gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                   >  1:1231704/1‑100 (MQ=255)

                         gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTCCGaaa                                   >  1:5503676/1‑100 (MQ=255)

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

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

                         gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACTTgggagatcggaa                                   >  1:516965/1‑90 (MQ=255)

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

                         cctgccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                   >  1:3511204/6‑100 (MQ=255)

                         atggccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                   >  1:1639762/6‑100 (MQ=255)

                         aaTGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAACAACAGACGAAGAGGAAAAAGCACATTCAATTACAGCAATTTGGaaagatcgga                                   >  1:3811973/2‑92 (MQ=255)

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

                         aaTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                   >  1:3608250/2‑100 (MQ=255)

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

                          aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGCTTACAGCAATTTGGAATTTCGaaaa                                  >  1:2462169/1‑100 (MQ=255)

                          aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGagatcggaaga                                  >  1:2707666/1‑90 (MQ=255)

                          aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGCATTTCGaaaa                                  >  1:1771898/1‑100 (MQ=255)

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

                          aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                  >  1:1314707/1‑100 (MQ=255)

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

                          aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                  >  1:5345404/1‑100 (MQ=255)

                          aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATGTGGaagatcggaag                                  >  1:6016257/1‑91 (MQ=255)

                          aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACTTGGaagatcggaga                                  >  1:3935401/1‑91 (MQ=255)

                           tgccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAg                                 >  1:2683243/4‑100 (MQ=255)

                           ggccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAg                                 >  1:1488398/4‑100 (MQ=255)

                            ggcGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                                >  1:3727135/3‑100 (MQ=255)

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

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

                            gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                                >  1:4069027/3‑100 (MQ=255)

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

                            gccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGc                                >  1:3721159/3‑100 (MQ=255)

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

                             ccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCa                               >  1:197020/2‑100 (MQ=255)

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

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

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

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

                                     ggTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattatt                       >  1:4709267/1‑100 (MQ=255)

                                     ggTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattatt                       >  1:1858775/1‑100 (MQ=255)

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

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

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

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

                                               ccaaTACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAACAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTaaaaaa             >  1:2235784/3‑96 (MQ=255)

                                                tAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTaaaaaat            >  1:5817114/1‑95 (MQ=255)

                                                 aaTACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAAAAAtt           >  1:1584261/1‑100 (MQ=255)

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

                                                     cAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAAAAAtttaag       >  1:3065342/1‑96 (MQ=255)

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

                                                          gatgaAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAAAAAtttaaggagtt  >  1:4465068/1‑91 (MQ=255)

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

AAAACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAACAACAGACGAAGAAGAAAAAGCACATTCAATTACCGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAGACATTAGTGTTGAC  >  NC_007790/846‑1003

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