D17 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007790	781	G→T	V168L (GTA→TTA)	SAUSA300_RS14705 →	replication protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007790	781	0	G	T	83.0%	632.9 / 119.2	336	V168L (GTA→TTA)	SAUSA300_RS14705	replication protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base T (250/28); minor base A (48/9); total (299/37)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 336 total aligned reads displayed.

CGGCAATCAAAGAGACCTCAAAATATTCGGTTAAATCATCTGATTTTTTAACTGATGATGATG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTGCTTAAACAAAAACATAAAATTTTAAACTTAGACGATGTCGA  >  NC_007790/701‑878

taacAATCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGaaaaggt                                                                                >  1:2750975/4‑97 (MQ=255)

gaccAATCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGaaaaggt                                                                                >  1:3737961/4‑97 (MQ=255)

gaacAATCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGaaaaggt                                                                                >  1:5576103/4‑97 (MQ=255)

ctgctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAgg                                                                                >  1:2943434/6‑100 (MQ=255)

ctgctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAgg                                                                                >  1:1713272/6‑100 (MQ=255)

 aacAATCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGaaaaggtt                                                                               >  1:2484140/3‑96 (MQ=255)

  gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                                                              >  1:5603125/4‑100 (MQ=255)

  gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                                                              >  1:4370099/4‑100 (MQ=255)

   gggggagggagaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGttt                                                                             >  1:6015640/9‑100 (MQ=255)

     ctCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGGGATCTAGAAAAAGGTTTGt                                                                           >  1:776831/2‑100 (MQ=255)

     cccAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTa                                                                          >  1:5887239/3‑100 (MQ=255)

     aTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATGAAAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTg                                                                            >  1:2803011/1‑100 (MQ=255)

     aTCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTCTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTa                                                                          >  1:1087240/1‑100 (MQ=255)

      acagaaggaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTAt                                                                         >  1:4196737/8‑100 (MQ=255)

         aaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACA‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGa                                                                      >  1:2701390/1‑100 (MQ=255)

            tgaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGaaaa                                                                   >  1:4534128/2‑100 (MQ=255)

            agaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGaaa                                                                    >  1:489545/1‑100 (MQ=255)

                    aaaTATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGtt                                                           <  1:1312032/100‑1 (MQ=37)

                    aaaTATTCCGTTAAAGCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGtt                                                           <  1:1792724/100‑1 (MQ=37)

                       ttctccgTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTaa                                                         >  1:1355200/7‑100 (MQ=25)

                        atTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGt                                                       >  1:1238017/1‑100 (MQ=25)

                          tccgTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTa                                                     >  1:1145648/4‑100 (MQ=25)

                               tAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt                                                 <  1:5911955/100‑1 (MQ=25)

                               tAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt                                                 <  1:730791/100‑1 (MQ=25)

                               tAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt                                                 <  1:5398156/100‑1 (MQ=25)

                                   ccatcTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGAtt                                            >  1:2285589/2‑100 (MQ=25)

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

                                         gactttttAACTGATGACGACG‑AAAAAAATCACGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGGTATGGTGGATTACTTaa                                      >  1:3712791/4‑100 (MQ=25)

                                         gactttttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                      >  1:2821448/4‑100 (MQ=25)

                                         gactttttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                      >  1:4485980/4‑100 (MQ=25)

                                         gactttttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                      >  1:3331396/4‑100 (MQ=25)

                                         gactttttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                      >  1:3517658/4‑100 (MQ=25)

                                         gactttttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                      >  1:254728/4‑100 (MQ=25)

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

                                         gactttttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                      >  1:2105153/4‑100 (MQ=25)

                                         gactttttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                      >  1:5761198/4‑100 (MQ=25)

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

                                         gactttttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaa                                      >  1:1066067/4‑100 (MQ=25)

                                         gactttttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGACTACTTaa                                      >  1:2620762/4‑100 (MQ=25)

                                         gactttttAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTa                                       >  1:67812/4‑100 (MQ=25)

                                         gactttttAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTa                                       >  1:897014/4‑100 (MQ=25)

                                             ttttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaa                                  >  1:4333931/1‑100 (MQ=25)

                                             gctccccTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaa                                  >  1:471303/7‑100 (MQ=25)

                                             ctttAACTGATGACGACGAAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTaaaga                                   >  1:5591185/2‑98 (MQ=25)

                                              cttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTACGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGAGGATTACTTAAAGaaa                                 >  1:875542/2‑100 (MQ=21)

                                               ttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaa                                >  1:5125405/1‑100 (MQ=25)

                                               ttAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaa                                >  1:5458921/1‑100 (MQ=25)

                                                tAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaag                               >  1:5530561/1‑99 (MQ=25)

                                                tAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGCGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaag                               >  1:6027167/1‑99 (MQ=25)

                                                caaCTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaag                               >  1:5188375/2‑99 (MQ=25)

                                                 agcTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaagc                              >  1:2977837/3‑98 (MQ=25)

                                                 aaCTGATGACGACGAAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGCTATGGCGGATTACTTAAAGaaaag                               >  1:3085151/1‑99 (MQ=21)

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

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

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

                                                   cTGTTGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCAt                            >  1:1010559/1‑100 (MQ=21)

                                                   cTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTGAAGAAAAGCAt                            >  1:3646464/1‑100 (MQ=25)

                                                   cTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCAt                            >  1:4910257/1‑100 (MQ=25)

                                                   cTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCAt                            >  1:3809079/1‑100 (MQ=25)

                                                   cTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCAt                            >  1:458975/1‑100 (MQ=25)

                                                    tgatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATa                           >  1:1873540/1‑100 (MQ=25)

                                                    tcccgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGAATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATa                           >  1:5486142/11‑100 (MQ=255)

                                                    agatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATa                           >  1:1623974/2‑100 (MQ=21)

                                                     ggtgacgacg‑aTAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:2999097/3‑100 (MQ=21)

                                                     gctgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:4012274/3‑100 (MQ=21)

                                                     gatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAACGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:2831850/1‑100 (MQ=21)

                                                     gatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAGAGAAAAGCATaa                          >  1:3240024/1‑100 (MQ=21)

                                                     gatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:764587/1‑100 (MQ=21)

                                                     gatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:808473/1‑100 (MQ=21)

                                                     gatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:5831688/1‑100 (MQ=21)

                                                     gatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:3352774/1‑100 (MQ=21)

                                                     gatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:1384121/1‑100 (MQ=21)

                                                     gatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:5009703/1‑100 (MQ=21)

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

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

                                                     gatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:4861944/1‑100 (MQ=21)

                                                     gatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTAGTATCGAAAACGAATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:3928237/1‑100 (MQ=255)

                                                     gatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGCTTGTATCGAAAACGACTGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:4266193/1‑100 (MQ=255)

                                                     gatgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGATTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                          >  1:2377609/1‑100 (MQ=21)

                                                      ctgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaaa                         >  1:2044752/2‑100 (MQ=21)

                                                      atgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaaa                         >  1:4224464/1‑100 (MQ=21)

                                                       tgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaaaa                        >  1:5302826/1‑100 (MQ=21)

                                                       tgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaaaa                        >  1:3428990/1‑100 (MQ=21)

                                                       tgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaaaa                        >  1:2141456/1‑100 (MQ=21)

                                                       cgacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaaaa                        >  1:4976006/8‑100 (MQ=21)

                                                        gacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAt                       >  1:1573335/7‑100 (MQ=21)

                                                        gacgacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAt                       >  1:1453126/7‑100 (MQ=21)

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

                                                         ttcacgaaaaaaaaTCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAt                       >  1:2192524/8‑100 (MQ=21)

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

                                                         accacg‑aAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAtt                      >  1:5228031/6‑100 (MQ=21)

                                                                      ctCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGa           >  1:2957041/2‑100 (MQ=25)

                                                                        cAAGAAATTTTAAGTGATCTAG‑GAAAGGTTTGGATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGa        >  1:3205416/1‑100 (MQ=21)

                                                                        cAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGa        >  1:4121396/1‑100 (MQ=21)

                                                                        cAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGa        >  1:4022919/1‑100 (MQ=21)

                                                                        cAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGata         >  1:377496/1‑98 (MQ=21)

                                                                             aaTTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACTTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCg   <  1:3398770/100‑1 (MQ=255)

                                                                             aaTTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATgccc   <  1:1218237/100‑4 (MQ=21)

                                                                              aTTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGa  <  1:2356436/100‑1 (MQ=21)

                                                                              aTTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGa  <  1:3746474/100‑1 (MQ=21)

CGGCAATCAAAGAGACCTCAAAATATTCGGTTAAATCATCTGATTTTTTAACTGATGATGATG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTGCTTAAACAAAAACATAAAATTTTAAACTTAGACGATGTCGA  >  NC_007790/701‑878

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: `‑`