D22 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007790	844	C→G	Q189E (CAA→GAA)	SAUSA300_RS14705 →	replication protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007790	844	0	C	G	90.3%	1562.7 / 159.8	630	Q189E (CAA→GAA)	SAUSA300_RS14705	replication protein
Reads supporting (aligned to +/- strand): ref base C (57/4); new base G (543/26); total (600/30)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-04

Only 100 of 631 total aligned reads displayed.

TTTTAACTGATGATGATGAAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTGCTTAAACAAAAACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAACAACAGACGAAGAAGAAAA  >  NC_007790/746‑935

 actAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaa                                                                                          >  1:3365107/3‑100 (MQ=25)

  ttAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaa                                                                                         >  1:1739411/1‑100 (MQ=25)

  gtAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaa                                                                                         >  1:4401680/2‑100 (MQ=25)

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

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

    taCTGATGACGACGAAAAAAATCAAGAAATTTTAAGCGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaagc                                                                                       >  1:1265784/2‑98 (MQ=25)

    gaCTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaagc                                                                                       >  1:4324269/2‑98 (MQ=25)

    aaCTGATGACGACGAAAAAAATCAAGAAACTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaagc                                                                                       >  1:4432002/1‑98 (MQ=21)

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

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

      cctctgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCAt                                                                                     >  1:4905269/5‑100 (MQ=21)

       tgatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATa                                                                                    >  1:3872055/1‑100 (MQ=25)

        gatgacgacgaAAAAAATCAAGAAATTTTAAGTGCTCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                                   >  1:831788/1‑100 (MQ=21)

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

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

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

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

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

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

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

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

        gatgacgacaaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                                   >  1:1571105/1‑100 (MQ=21)

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

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

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

           cacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAt                                                                                >  1:2065812/7‑100 (MQ=21)

           acgacgaaaaaaaaTCAAGAAATTTTAAGTGATCTAG‑AAAAGGGTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAt                                                                                >  1:4633195/8‑100 (MQ=21)

            ccgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAtt                                                                               >  1:5104890/6‑100 (MQ=21)

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

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

                        aTCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGAc                                                                   >  1:3281807/1‑100 (MQ=21)

                                aTTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATgcgga                                                           <  1:3019208/100‑5 (MQ=21)

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

                                 ttatcaGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGa                                                           <  1:4133746/95‑1 (MQ=18)

                                     caGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGAt                                                       >  1:2416918/2‑100 (MQ=25)

                                     aaGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAGGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATg                                                      >  1:5207277/1‑100 (MQ=21)

                                      ggTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATgg                                                     >  1:3563295/2‑100 (MQ=21)

                                      aGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATgg                                                     >  1:3974866/1‑100 (MQ=25)

                                             ctagaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAAttt                                               >  1:1293365/5‑100 (MQ=25)

                                              aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGAtcggaagaca                                              >  1:3270673/1‑91 (MQ=21)

                                              aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAttta                                              >  1:4363551/1‑99 (MQ=25)

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

                                              aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAttta                                              >  1:1152662/1‑99 (MQ=25)

                                              aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAttta                                              >  1:1353266/1‑99 (MQ=25)

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

                                              aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATATTGAAGAtcggaagagc                                              >  1:2523621/1‑91 (MQ=18)

                                               agaaaaGGTTTGTATCGGAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                             >  1:4499830/3‑100 (MQ=25)

                                               agaaaaGGTTTGTATCGAAAACGTATGTTATGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGCAATTTGa                                             >  1:4575011/3‑100 (MQ=25)

                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGATGATGGACATTTGa                                             >  1:5174688/3‑100 (MQ=21)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGAAATTTGa                                             >  1:1865235/3‑100 (MQ=25)

                                                 aaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGAtt                                           >  1:3674312/1‑100 (MQ=37)

                                                             ccGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGtgat                               >  1:5031323/2‑100 (MQ=37)

                                                              ngAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGtgatg                              >  1:2563510/2‑100 (MQ=37)

                                                              cGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGACGGTAATTTGATTAATACAAGtgatg                              <  1:1709081/100‑1 (MQ=37)

                                                                  ctcGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGa                          >  1:2815527/3‑100 (MQ=255)

                                                                  aaCGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGa                          >  1:4267485/1‑100 (MQ=38)

                                                                   aCGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGAt                         >  1:629049/1‑100 (MQ=38)

                                                                    cGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATa                        >  1:4465781/1‑100 (MQ=25)

                                                                       aTGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATaaga                     >  1:3399164/1‑98 (MQ=255)

                                                                         gTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAaaca                   >  1:1390916/1‑100 (MQ=255)

                                                                           tAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAgacgac                 <  1:949061/100‑1 (MQ=255)

                                                                            aaGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAacaaca                >  1:1053711/1‑100 (MQ=255)

                                                                            aaGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAATTTAATTAATACAAGTGACGAAGATAAAacaaca                >  1:4830429/1‑100 (MQ=255)

                                                                             ggTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACTATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAg               >  1:1642547/2‑100 (MQ=255)

                                                                                gaTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACg            >  1:3224719/2‑100 (MQ=255)

                                                                                caTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAGGACAGatc            >  1:3120740/2‑98 (MQ=255)

                                                                                 aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAacgacaagaga           >  1:2886136/1‑95 (MQ=255)

                                                                                 aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACga           >  1:4789539/1‑100 (MQ=255)

                                                                                 aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACAACAGatcg           >  1:3471113/1‑97 (MQ=255)

                                                                                     tggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACTGACgaagat       >  1:2326862/1‑99 (MQ=255)

                                                                                          tacTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACGAAGAAGaaaa  <  1:2225754/98‑1 (MQ=255)

TTTTAACTGATGATGATGAAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTGCTTAAACAAAAACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAACAACAGACGAAGAAGAAAA  >  NC_007790/746‑935

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