D17 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007790	837	G→A	L186L (TTG→TTA)	SAUSA300_RS14705 →	replication protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007790	837	0	G	A	89.6%	1487.5 / 108.6	695	L186L (TTG→TTA)	SAUSA300_RS14705	replication protein
Reads supporting (aligned to +/- strand): ref base G (68/4); new base A (586/37); total (654/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 695 total aligned reads displayed.

CTGATTTTTTAACTGATGATGATGAAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTGCTTAAACAAAAACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAACAACAGACGAAGAAG  >  NC_007790/740‑931

cTGACTTTTTAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACtt                                                                                             >  1:1959118/1‑100 (MQ=25)

  gactttttAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                                                                           >  1:2080621/4‑100 (MQ=25)

  gactttttAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                                                                           >  1:2754568/4‑100 (MQ=25)

  gactttttAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                                                                           >  1:2629105/4‑100 (MQ=25)

  gactttttAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                                                                           >  1:2588401/4‑100 (MQ=25)

  gactttttAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                                                                           >  1:740918/4‑100 (MQ=25)

  gactttttAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                                                                           >  1:5067127/4‑100 (MQ=25)

  gactttttAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                                                                           >  1:3831981/4‑100 (MQ=25)

  gactttttAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                                                                           >  1:1472029/4‑100 (MQ=25)

      agttAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGCGGATTACTTAAAGaa                                                                                       >  1:3338585/3‑100 (MQ=21)

         tAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaag                                                                                    >  1:282067/1‑99 (MQ=25)

         tAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaag                                                                                    >  1:5989366/1‑99 (MQ=25)

         tAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaag                                                                                    >  1:1066146/1‑99 (MQ=25)

          aaCTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaagc                                                                                   >  1:2521793/1‑98 (MQ=25)

          aaCTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaagc                                                                                   >  1:613247/1‑98 (MQ=25)

           aCTGATGACGACGAAAAAAATCAGGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCa                                                                                  >  1:1395574/1‑100 (MQ=21)

           aCTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCa                                                                                  >  1:4776408/1‑100 (MQ=25)

            cTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCAt                                                                                 >  1:2464909/1‑100 (MQ=25)

              tatagcgatgaAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATa                                                                                >  1:922305/7‑100 (MQ=25)

              gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:2040061/1‑100 (MQ=21)

              gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:5149247/1‑100 (MQ=21)

              gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:4407391/1‑100 (MQ=21)

              gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:5522259/1‑100 (MQ=21)

              gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:1229479/1‑100 (MQ=21)

              gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:803382/1‑100 (MQ=21)

              gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:2721919/1‑100 (MQ=21)

              gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGCTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:5578557/1‑100 (MQ=255)

                 gacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAt                                                                            >  1:4566596/7‑100 (MQ=21)

                 gacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAt                                                                            >  1:5769272/7‑100 (MQ=21)

                  acgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAtt                                                                           >  1:6050046/6‑100 (MQ=21)

                   cgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAttt                                                                          >  1:3715407/5‑100 (MQ=21)

                   agacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAttt                                                                          >  1:1898646/5‑100 (MQ=21)

                          aaaaaTCAAGGAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACtt                                                                   <  1:4844054/100‑1 (MQ=21)

                          aaaaaTCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACtt                                                                   <  1:847562/100‑1 (MQ=21)

                                           ctgTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATg                                                  >  1:3785068/3‑100 (MQ=21)

                                            cgTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCCAAGAtcg                                                 >  1:1542636/2‑98 (MQ=18)

                                             gTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGt                                                <  1:177920/100‑1 (MQ=21)

                                                gatctagaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTaa                                              >  1:3217636/8‑100 (MQ=21)

                                                    aGAAAAAGTTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAttta                                          >  1:1140344/1‑99 (MQ=21)

                                                    aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAttta                                          >  1:4511370/1‑99 (MQ=21)

                                                    aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAttta                                          >  1:3959395/1‑99 (MQ=21)

                                                    aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAttta                                          >  1:1429083/1‑99 (MQ=21)

                                                    aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAttta                                          >  1:408048/1‑99 (MQ=21)

                                                    aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAttta                                          >  1:3019358/1‑99 (MQ=21)

                                                    aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCCTAAACAAAAGCATAAAATTTTAAATTTAGATAACGTTGAAGAtcggaagagc                                          >  1:4584333/1‑91 (MQ=18)

                                                     agaagaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:318026/6‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGACAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTCGa                                         >  1:5323550/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGTTTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1192773/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:2190973/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:4655043/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:472271/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:2143399/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:2004081/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:5636371/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1623676/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:49264/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:4221578/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:4994327/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1658359/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1689225/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:5464905/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:5256710/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:5441198/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:3746536/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:293159/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1306875/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:692956/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:3385802/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:3451401/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1330988/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:100123/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:3899304/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:4019107/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:5861043/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:5825988/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:4351406/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:2222167/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATCTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:4462576/3‑100 (MQ=25)

                                                     agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGGGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGACGATGGTAATTTGa                                         >  1:3605730/3‑100 (MQ=21)

                                                     agaaaaGGCTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:2344443/3‑100 (MQ=25)

                                                      acaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGAt                                        >  1:2966335/3‑100 (MQ=25)

                                                        ctaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTa                                      >  1:5193177/3‑100 (MQ=25)

                                                        aaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTa                                      >  1:5921570/1‑100 (MQ=25)

                                                        aaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTa                                      >  1:307452/1‑100 (MQ=25)

                                                        aaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTa                                      >  1:2313814/1‑100 (MQ=25)

                                                         aaGGTTTGTATCGAAAACGTATGTTAAGCTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTaa                                     >  1:4971090/1‑100 (MQ=25)

                                                             ggtGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACa                                 >  1:5717583/3‑100 (MQ=25)

                                                             cttGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACa                                 >  1:3301031/2‑100 (MQ=25)

                                                               tGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAg                               >  1:4283474/1‑100 (MQ=37)

                                                                 tATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAATTTAATTAATACAAGtg                             >  1:892412/1‑100 (MQ=21)

                                                                         ccGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGAt                     >  1:4184402/2‑100 (MQ=255)

                                                                          tgTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATa                    >  1:5107982/2‑100 (MQ=255)

                                                                             aTGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATaaaa                 >  1:1860461/1‑100 (MQ=255)

                                                                                     tgaTGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAGACAACAGat         >  1:2367742/3‑99 (MQ=255)

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

                                                                                      tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACg        >  1:5365005/1‑100 (MQ=255)

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

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

                                                                                            ggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACgaagaag  >  1:4889856/1‑100 (MQ=255)

                                                                                            agATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACgaagaag  >  1:3138682/2‑100 (MQ=255)

CTGATTTTTTAACTGATGATGATGAAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTGCTTAAACAAAAACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAACAACAGACGAAGAAG  >  NC_007790/740‑931

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