BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	657,807	G→A	6.1%	G75S (GGT→AGT)	USA300HOU_RS03150 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	657,807	0	G	A	6.1%	343.4 / 8.6	117	G75S (GGT→AGT)	USA300HOU_RS03150	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (51/56); new base A (5/2); total (58/59)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 120 total aligned reads displayed.

TAATAAGCACTTGGATTAGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGCCTAATTAAAAAA  >  USA300TCH1516_ALE/657667‑657947

taataaGCACTTGGATTAGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAg                                                                                                                                              <  1:985606/141‑1 (MQ=255)

       cACTTGGATTAGTTTATTTTTTTTTTTTTTTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGTTAtctt                                                                                                                                       >  4:32727/1‑140 (MQ=255)

        aCTTGGATTAGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATctct                                                                                                                                      <  3:356393/141‑1 (MQ=255)

        aCTTGGATTAGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTACTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATctct                                                                                                                                      <  4:327482/141‑1 (MQ=255)

          ttGGATTAGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTTCAGGTATCTCTAc                                                                                                                                    >  4:107580/1‑141 (MQ=255)

           tGGATTAGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACt                                                                                                                                   <  4:332296/141‑1 (MQ=255)

           tGGATTAGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACt                                                                                                                                   >  3:281570/1‑141 (MQ=255)

           tGGATTAGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACt                                                                                                                                   <  1:807792/141‑1 (MQ=255)

           tGGATTAGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCAAAACGATAAGCCAAAATTAGTAGAGAAGCGTTTGCAAAATTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACt                                                                                                                                   <  1:648309/141‑1 (MQ=11)

               ttAGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGTTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTAt                                                                                                                               >  1:556212/1‑141 (MQ=255)

                tAGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTTAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCTAGTATCTCTACTTTAtt                                                                                                                              >  1:292290/1‑141 (MQ=255)

                 aGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTc                                                                                                                             >  4:330267/1‑141 (MQ=255)

                 aGTTTATTTATTCCTCTTATTAATTGGATCAACCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAAACTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTc                                                                                                                             <  2:542795/141‑1 (MQ=17)

                   tttATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGc                                                                                                                           <  3:77421/141‑1 (MQ=255)

                   tttATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGc                                                                                                                           >  1:557890/1‑141 (MQ=255)

                   tttATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGc                                                                                                                           >  1:557917/1‑141 (MQ=255)

                       ttttttCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGtt                                                                                                                       >  1:690786/1‑141 (MQ=255)

                        tttttCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTa                                                                                                                      >  1:1077029/1‑141 (MQ=255)

                         ttttCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTAc                                                                                                                     <  2:38626/141‑1 (MQ=255)

                              cTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTaa                                                                                                                <  2:984283/141‑1 (MQ=255)

                              cTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTaa                                                                                                                >  3:63309/1‑141 (MQ=255)

                                 tttttaATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTTGAGAAGAGTTTGAATATTTAAATATTTTGATTTTTTTTTTTTTTAATGTTTTTTGTCCATTTTCAGTTATCTCTACTTTATTCGCAGTTACATTAAGaa                                                                                                             >  4:365731/4‑141 (MQ=255)

                                  tattaATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTGTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGaaa                                                                                                            <  3:420528/141‑1 (MQ=255)

                                     taATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGta                                                                                                         >  3:173059/1‑141 (MQ=255)

                                        ttGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTTAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCTAGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATac                                                                                                      >  4:58288/1‑141 (MQ=255)

                                        ttGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGTCTTGTTGCAGGTATCTCTACTTTATTTGCAGTTACATTAAGAAAGTATac                                                                                                      >  2:826652/1‑141 (MQ=255)

                                        ttGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATac                                                                                                      >  1:1051422/1‑141 (MQ=255)

                                         tGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAg                                                                                                                         >  2:93137/1‑121 (MQ=255)

                                         tGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATaca                                                                                                     >  3:257338/1‑141 (MQ=255)

                                         tGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATaca                                                                                                     <  2:523797/141‑1 (MQ=255)

                                           gTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATAcaca                                                                                                   <  3:308607/141‑1 (MQ=255)

                                           gTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATAcaca                                                                                                   <  1:993422/141‑1 (MQ=255)

                                            ttCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTTAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCTAGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATacgctt                                                                                                  >  1:823826/1‑137 (MQ=255)

                                             tCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTTTTTTTTTTCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACAtt                                                                                                 >  2:303024/1‑141 (MQ=21)

                                               aTCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTa                                                                                               >  1:423001/1‑141 (MQ=255)

                                                tCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTTAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCTAGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACGCTTTTAc                                                                                              <  2:823826/141‑1 (MQ=18)

                                                tCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTAc                                                                                              >  2:158859/1‑141 (MQ=255)

                                                   ccAAAAGATACGTCAAAACTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTg                                                                                           <  4:120365/141‑1 (MQ=255)

                                                    cAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAATAAGGGTTTTGAGCCAGTTGCAGGTATCTCTACTTTAGTTGCAGTTACATTAAGAAAGTATACACATTTACTTGa                                                                                          >  1:377265/1‑141 (MQ=17)

                                                    cAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGa                                                                                          >  2:1063434/1‑141 (MQ=255)

                                                    cAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGa                                                                                          >  3:143329/1‑141 (MQ=255)

                                                    cAAAAGATACGTAAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGa                                                                                          <  2:780642/141‑1 (MQ=255)

                                                       aaGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGacac                                                                                       <  2:312410/141‑1 (MQ=255)

                                                       aaGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGacac                                                                                       <  2:312421/141‑1 (MQ=255)

                                                       aaGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGacac                                                                                       <  1:842400/141‑1 (MQ=255)

                                                       aaGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGacac                                                                                       <  1:842416/141‑1 (MQ=255)

                                                        aGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACt                                                                                      <  2:456143/141‑1 (MQ=255)

                                                        aGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACt                                                                                      <  2:615034/141‑1 (MQ=255)

                                                         gATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTc                                                                                     <  1:389463/141‑1 (MQ=255)

                                                          aTACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCa                                                                                    >  1:1086574/1‑141 (MQ=255)

                                                            aCGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAAc                                                                                  >  4:175196/1‑141 (MQ=255)

                                                             cGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTTTTTTTTTATTGGTATCTATACTTTATTCGCAGTTACATTATGAAAGTATACACATATACTTGACACTCACCt                                                                                 >  2:1020366/1‑141 (MQ=18)

                                                             cGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACt                                                                                 <  1:1067487/141‑1 (MQ=255)

                                                             cGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACt                                                                                 >  1:781031/1‑141 (MQ=255)

                                                              gTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACtt                                                                                <  1:129640/141‑1 (MQ=255)

                                                                cAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGa                                                                              <  2:991457/141‑1 (MQ=255)

                                                                 aaaaTTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGAt                                                                             <  4:55401/141‑1 (MQ=255)

                                                                 aaaaTTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGAt                                                                             <  2:839549/141‑1 (MQ=255)

                                                                       agtagAGAAGAGTTTGAAAATTTAAATATTGTTAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCTAGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACGCTTTTACTTGACACTCAACTTGATaaaaaa                                                                       <  3:318713/141‑1 (MQ=17)

                                                                        gtagAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTTTTTTTTTTTTTGTATCTCTACTTTATTCACATTTATATTAAGAAATTATACACATTTACTTGACATTCAACTTGATAGaaatt                                                                      >  2:980313/1‑139 (MQ=17)

                                                                          agagAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTa                                                                    >  1:830542/1‑141 (MQ=255)

                                                                           gagaAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAg                                                                   <  1:314675/141‑1 (MQ=255)

                                                                            agaagaGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGTCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGt                                                                  >  4:233967/1‑141 (MQ=255)

                                                                                gagTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTAtt                                                              <  3:327506/141‑1 (MQ=255)

                                                                                  gTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGc                                                            <  1:394627/141‑1 (MQ=255)

                                                                                   tttGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCt                                                           <  4:376990/141‑1 (MQ=255)

                                                                                    ttGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCta                                                          >  4:271772/1‑141 (MQ=255)

                                                                                       aaaaTTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTTTTTTTTTTTATGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTatat                                                       >  2:195748/1‑141 (MQ=21)

                                                                                            ttAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATAtgttgt                                                  >  4:348829/1‑141 (MQ=255)

                                                                                            ttAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATAtgttgt                                                  >  2:658147/1‑141 (MQ=255)

                                                                                             tAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTa                                                 <  2:557890/141‑1 (MQ=255)

                                                                                             tAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTa                                                 <  2:557917/141‑1 (MQ=255)

                                                                                              aaaTATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTAt                                                >  3:193493/1‑141 (MQ=255)

                                                                                               aaTATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATc                                               >  4:258595/1‑141 (MQ=255)

                                                                                                 tatTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATcac                                             >  1:280814/1‑141 (MQ=255)

                                                                                                    tgtgAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACAtt                                          <  3:7014/141‑1 (MQ=255)

                                                                                                     gttaaaCCTGCTAAAAAGAATGTTTTTTTTCTTTTTTTTAGTATCTATACTTTATTCGCAGTTACATTAAGAAAGTATACGCTTTTACTTGACACTCAACTTGATAAAAAATTAGTTATTGCCATATGTTGCACCACAttt                                         >  4:202074/4‑141 (MQ=14)

                                                                                                      tgAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTta                                        <  1:977810/141‑1 (MQ=255)

                                                                                                         aaCCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAg                                     <  4:148010/141‑1 (MQ=255)

                                                                                                          aCCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAgg                                    >  3:325294/1‑141 (MQ=255)

                                                                                                            cTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGa                                  <  4:43535/141‑1 (MQ=255)

                                                                                                             tGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGAt                                 >  1:64152/1‑141 (MQ=255)

                                                                                                             tGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGAt                                 >  1:64176/1‑141 (MQ=255)

                                                                                                               cTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTAATTTATTCGTAATTACACAAAGCAAGTACCAACATTTCCTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGAttt                               <  3:42371/141‑1 (MQ=18)

                                                                                                                 aaaaaGAATTTTTTTTGGGCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTa                             <  1:518492/141‑1 (MQ=255)

                                                                                                                  aaaaGAATGTTTTTTGGCCAGTTGCTAGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACGCTTTTACTTGACACTCAACTTGATAAAAAATTAGTTATTGCCATATGTTGCACCACATTTATAGGGATTTTaa                            >  1:297733/1‑141 (MQ=17)

                                                                                                                         tGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTa                     >  2:164849/1‑141 (MQ=255)

                                                                                                                             ttttGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTa                 >  1:192072/1‑141 (MQ=255)

                                                                                                                             ttttGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTa                 >  2:341054/1‑141 (MQ=255)

                                                                                                                                 ggCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGcc             <  3:242199/141‑1 (MQ=255)

                                                                                                                                  gCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGCCt            >  2:1077914/1‑141 (MQ=255)

                                                                                                                                      gTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGCCTAAtt        <  1:1072717/141‑1 (MQ=255)

                                                                                                                                        tGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGCCTAATTaa      <  2:129054/141‑1 (MQ=255)

                                                                                                                                          cAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGCCTAATTaaaa    >  2:767568/1‑141 (MQ=255)

                                                                                                                                          cAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGACTAAATaaaa    <  4:272895/141‑1 (MQ=255)

                                                                                                                                           aGGTATCTCTACTTTATTCGCAGTTACATTAAGAGAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGCCTAATTaaaaa   >  2:667901/1‑141 (MQ=255)

                                                                                                                                           aGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGCCTAATTaaaaa   >  1:622572/1‑141 (MQ=255)

                                                                                                                                           aGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGCCTAATTaaaaa   >  2:32116/1‑141 (MQ=255)

                                                                                                                                            ggTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGCCTAATTaaaaaa  <  4:267440/141‑1 (MQ=255)

                                                                                                                                            ggTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGCCTAATTaaaaaa  <  2:502044/141‑1 (MQ=255)

TAATAAGCACTTGGATTAGTTTATTTTTTCCTTTTATTAATTGGTTCATCCCCAAAAGATACGTCAAAATTAGTAGAGAAGAGTTTGAAAATTTAAATATTGTGAAACCTGCTAAAAAGAATGTTTTTTGGCCAGTTGCAGGTATCTCTACTTTATTCGCAGTTACATTAAGAAAGTATACACATTTACTTGACACTCAACTTGATAGAAAATTAGTTATTGCTATATGTTGTATCACATTTATAGGGATTTTAACATTTTATGTACGCCTAATTAAAAAA  >  USA300TCH1516_ALE/657667‑657947

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