D22 :: 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	96.5%	305.2 / ‑1.3	115	S222S (TCA→TCG)	SAUSA300_RS14705	replication protein
Reads supporting (aligned to +/- strand): ref base A (3/1); new base G (108/3); total (111/4)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01

Only 100 of 115 total aligned reads displayed.

AACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAACAACAGACGAAGAAGAAAAAGCACATTCAATTACCGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAGACATTAGTGTTGA  >  NC_007790/848‑1002

agcATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAACAACAGACGAAGAGGAAAAAGCACATTCAAt                                                         >  1:2622231/3‑100 (MQ=255)

agcATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGAt                                                         >  1:2344990/3‑100 (MQ=255)

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

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

      aaaTTTTAAATTTAGATAATGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAACAACAGACGAAGAGGAAAAAGCACATTCAATTACAGc                                                   <  1:4782792/100‑1 (MQ=255)

       ccttttAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCa                                                  >  1:931462/3‑100 (MQ=255)

          tgtAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCaaga                                               >  1:113435/3‑98 (MQ=255)

            tAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTAcgcagatcgg                                             >  1:2037708/1‑91 (MQ=255)

             aaaCTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTgg                                            <  1:4096352/100‑1 (MQ=255)

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

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

             aaaCTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTAgg                                            >  1:2194246/1‑100 (MQ=255)

             aaaCTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATATgg                                            >  1:1933850/1‑100 (MQ=255)

             aaaCCTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCaaagatc                                            >  1:1016498/1‑95 (MQ=255)

               gcctAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGaa                                          >  1:4849241/4‑100 (MQ=255)

                  tagataaTGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAACAACAGACGAAGAGGAAAAAGCACATTCAATTACAGCAATCTGGaaaga                                       >  1:926568/7‑97 (MQ=255)

                  tAGACNATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAAttt                                       >  1:1290492/1‑100 (MQ=255)

                  tAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCaagatcggaaga                                       >  1:3381636/1‑90 (MQ=255)

                  tAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAAttt                                       >  1:1190852/1‑100 (MQ=255)

                  tAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAAttt                                       >  1:3488724/1‑100 (MQ=255)

                  tAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAAttt                                       >  1:4582547/1‑100 (MQ=255)

                   aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGGCAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTc                                      >  1:3075960/1‑100 (MQ=255)

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

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

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

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

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

                   aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACTTgggagatc                                      >  1:3912352/1‑94 (MQ=255)

                   aGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACTTgggagatc                                      >  1:610929/1‑94 (MQ=255)

                   aGACGATGCCGAAGATGGCAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGaaagat                                      >  1:1647712/1‑96 (MQ=255)

                     aCGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGGCGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAAtaagat                                    >  1:2391716/1‑95 (MQ=255)

                     aCGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTCGGagagatcg                                    >  1:5224368/1‑93 (MQ=255)

                     aCGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTAGGAAtagatc                                    >  1:3849257/1‑95 (MQ=255)

                      tgATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                   >  1:4126783/2‑100 (MQ=255)

                      tgATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                   >  1:4947558/2‑100 (MQ=255)

                      tgATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                   >  1:1829784/2‑100 (MQ=255)

                      tgATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaa                                   >  1:3567782/2‑100 (MQ=255)

                      tgATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTCGGaagatcgga                                   >  1:615312/2‑93 (MQ=255)

                      ccaTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATCTGGaagatcgga                                   >  1:1615269/3‑93 (MQ=255)

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

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

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

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

                      cGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGACTACAGCAATTTGGAATTTCGaa                                   >  1:1626372/1‑100 (MQ=255)

                       naTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                  >  1:3431217/2‑100 (MQ=255)

                       gcggccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                  >  1:4639889/6‑100 (MQ=255)

                       gATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGaaagatcgga                                  >  1:1384277/1‑92 (MQ=255)

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

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

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

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

                       gATGCCGAAAATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGaagatcggaa                                  >  1:179832/1‑92 (MQ=255)

                       caTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                  >  1:537447/2‑100 (MQ=255)

                       caTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                  >  1:877606/2‑100 (MQ=255)

                       actgccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaa                                  >  1:856309/6‑100 (MQ=255)

                       aaTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACTTGGaagatcggaa                                  >  1:3821780/2‑92 (MQ=255)

                        ttgccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                 >  1:3999501/5‑100 (MQ=255)

                        cctccGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGaaaa                                 >  1:323055/5‑100 (MQ=255)

                        aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCGCATTCGATTACAGCAATTAGGaagatcggaag                                 >  1:4503372/1‑91 (MQ=255)

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

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

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

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

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

                        aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAAGTTGGaagatcggaag                                 >  1:2531838/1‑91 (MQ=255)

                        aTGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAACTTGGaagatcggaag                                 >  1:3766346/1‑91 (MQ=255)

                         tgccGAAGGTGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAg                                >  1:1557283/4‑100 (MQ=255)

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

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

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

                         tgccGAAGATGGAAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAg                                >  1:776867/4‑100 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                               ccaTGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAat                          >  1:1333145/3‑100 (MQ=255)

                                tgtGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAatt                         >  1:4540774/3‑100 (MQ=255)

                                actGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAatt                         >  1:496171/3‑100 (MQ=255)

                                 ttGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAatta                        >  1:4682473/2‑100 (MQ=255)

                                 aTGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAGGCAAAatta                        >  1:1370015/1‑100 (MQ=255)

                                  tGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattat                       >  1:1964286/1‑100 (MQ=255)

                                  aggTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattat                       >  1:3608802/2‑100 (MQ=255)

                                   ggTAATTTGATTAATACAAGTGATGAAGATAAAACGACATACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAattatt                      >  1:1220072/1‑100 (MQ=255)

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

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

                                           gATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTaaaa              >  1:2262278/1‑98 (MQ=255)

                                               aaTACAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTaaaaaatg          >  1:48856/1‑94 (MQ=255)

                                                atacAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAAAAAttg         >  1:3332414/4‑99 (MQ=255)

                                                  acAAGTGATGAAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAAAAAtttaa       >  1:4788123/2‑97 (MQ=255)

                                                   cAAGTGACGAAGATAAAACAACAGACGAAGAGGAAAAAGCACATTCAATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAGTCATCAgat      >  1:2399702/1‑98 (MQ=255)

                                                       tgatgaAGATAAAACGACAGACGAAGAAGAAAAAGCACATTCGATTACAGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAAAAAtttaaggagt  >  1:4715017/1‑92 (MQ=255)

AACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAACAACAGACGAAGAAGAAAAAGCACATTCAATTACCGCAATTTGGAATTTCGAAAAGCAAAATTATTATTTAAGACATTAGTGTTGA  >  NC_007790/848‑1002

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