BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	870,892	G→T	7.3%	S51R (AGC→AGA)	USA300HOU_RS04250 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	870,892	0	G	T	7.3%	484.2 / 3.8	153	S51R (AGC→AGA)	USA300HOU_RS04250	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (46/91); new base T (13/2); total (60/93)
Fisher's exact test for biased strand distribution p-value = 9.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

AGCAGTAATTTCCATTTTAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAGTGTACTTGCTAATAATATCCCCAT  >  USA300TCH1516_ALE/870754‑871032

aGCAGTAATTTCCATTTTAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTACCTTTTATATCCATA‑TTTTTTTTGCt                                                                                                                                              <  3:339811/141‑1 (MQ=255)

aGCAGTAATTTCCATTTTAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCt                                                                                                                                              >  2:397405/1‑141 (MQ=255)

     tAATTTCCATTTTAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTaa                                                                                                                                         <  3:363558/141‑1 (MQ=255)

      aaTTTCCATTTTAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTata                                                                                                                                        >  2:850097/1‑139 (MQ=255)

      aaTTTCCATTTTAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTaaa                                                                                                                                        <  1:1204349/141‑1 (MQ=255)

         tgtgaTTTTAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTc                                                                                                                                     <  1:1249072/137‑1 (MQ=255)

          tCCATTTTAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCt                                                                                                                                    >  2:549853/1‑141 (MQ=255)

          tCCATTTTAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCt                                                                                                                                    <  1:293642/141‑1 (MQ=255)

            cATTTTAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCttt                                                                                                                                  <  2:132325/141‑1 (MQ=255)

              ttttAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTaa                                                                                                                                <  2:159752/141‑1 (MQ=255)

              ttttAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTaa                                                                                                                                <  1:32829/141‑1 (MQ=255)

               tttAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTTATTTCTTTAAg                                                                                                                               >  2:563273/1‑141 (MQ=255)

                 tAAAAACGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGTCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGtt                                                                                                                             <  2:967599/141‑1 (MQ=255)

                  aaaaGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTAAATAAAACTCGCATAGACTCTAATTTTTATATTAATA‑TTTTTTTTGCTTTTTACTTCTTTAAGttt                                                                                                                            >  3:240032/1‑140 (MQ=255)

                    aaGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGttctt                                                                                                                          <  1:1019889/141‑1 (MQ=255)

                    aaGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGttctt                                                                                                                          <  2:656065/141‑1 (MQ=255)

                     aGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTTTTTTTTTTTTTTTTAAGtttttt                                                                                                                         >  2:101934/1‑140 (MQ=255)

                      gCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCa                                                                                                                        >  1:306398/1‑141 (MQ=255)

                       cGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTTTATTCTTTATGTTCTTCAt                                                                                                                       >  1:1185928/1‑141 (MQ=255)

                       cGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCAt                                                                                                                       <  2:495624/141‑1 (MQ=255)

                        gATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCAtt                                                                                                                      <  1:41371/141‑1 (MQ=255)

                         atatatTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCAttt                                                                                                                     <  4:407280/141‑1 (MQ=255)

                           atatTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTaa                                                                                                                   <  2:1437161/141‑1 (MQ=255)

                           atatTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTaa                                                                                                                   <  1:1081306/141‑1 (MQ=255)

                                   ttttCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCAt                                                                                                           <  1:1617889/141‑1 (MQ=255)

                                     ttCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTg                                                                                                         <  4:38207/141‑1 (MQ=255)

                                         cATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑‑TTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCAtt                                                                                                    >  4:637357/1‑141 (MQ=255)

                                           tCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCAttt                                                                                                   <  3:477319/141‑1 (MQ=255)

                                           tCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCAttt                                                                                                   <  2:627150/141‑1 (MQ=255)

                                             tttttCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATtttc                                                                                                 <  1:259441/141‑1 (MQ=255)

                                               tttCATTGATGATAAACGGTTCTTTTATTTTAGCTTTAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTACGTTTTTTTGTGGAGCCTTCATTGTCATtttctt                                                                                               <  2:836537/141‑1 (MQ=255)

                                                  cATTTATGATAAACGGTTTTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATATTTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATtttctt‑‑tc                                                                                           <  3:36389/141‑1 (MQ=255)

                                                  cATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTTTTTTTTTTTTCTTTAAGTTTTTCATTTTATTTTTTATTGTTGTTttctt‑‑tct                                                                                          >  2:697406/1‑141 (MQ=255)

                                                   aTTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATATTTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTttctt‑‑tct                                                                                          >  1:769890/1‑141 (MQ=255)

                                                   aTTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTTCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTtctt‑‑tctt                                                                                         >  2:99828/1‑141 (MQ=255)

                                                   aTTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTtctt‑‑tctt                                                                                         >  3:515362/1‑141 (MQ=255)

                                                    ttGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCCTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTctt‑‑tcttt                                                                                        >  1:412943/1‑141 (MQ=255)

                                                       atgatAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTACCAATAGTATCGAATTCCATTTTCATCCCTATT‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTtg                                                                                     <  1:1401911/141‑1 (MQ=255)

                                                       atgatAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTtg                                                                                     <  1:1019834/141‑1 (MQ=255)

                                                        tgatAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTtgt                                                                                    <  4:65397/141‑1 (MQ=255)

                                                        tgatAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTtgt                                                                                    <  4:65524/141‑1 (MQ=255)

                                                          atAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTTTTTTTTTTTTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTtgtgt                                                                                  >  2:1327619/1‑140 (MQ=255)

                                                            aaaCGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACt                                                                                <  3:416612/141‑1 (MQ=255)

                                                               cGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGt                                                                             >  1:1088465/1‑141 (MQ=255)

                                                                ggTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTTCTTTTTTTTTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGtt                                                                            >  4:208802/1‑140 (MQ=255)

                                                                 gTTCTTTTATTTTAGCTTCAAATTTGTCACTATCAATAGTATCGCCTTTAATTTTTATTTCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCTTTTTCTT‑‑TCTTTTTGTGACTAGTGc                                                                           <  4:709339/141‑1 (MQ=255)

                                                                  ttCTTTTATTTTAGCTTCTAATTTGTCACTTATAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTTCTTTTTAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGct                                                                          <  2:1197108/141‑1 (MQ=255)

                                                                  ttCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTTTTTTTTTTTTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTTTT‑‑TTTTTTTTTGTCTAGTGct                                                                          >  2:1362262/1‑141 (MQ=255)

                                                                   tCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGctc                                                                         <  1:506774/141‑1 (MQ=255)

                                                                    cTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGctct                                                                        <  1:1036067/141‑1 (MQ=255)

                                                                    cTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGctct                                                                        <  4:326765/141‑1 (MQ=255)

                                                                     ttttattttaGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTTCTTTTAATTTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGctctc                                                                       >  2:231456/1‑141 (MQ=255)

                                                                     ttttattttaGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGATTTGTGctttt                                                                       >  4:192152/1‑138 (MQ=255)

                                                                     ttttattttaGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGctctc                                                                       >  2:1642623/1‑141 (MQ=255)

                                                                      tttattttaGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTATATTCTTTAAGTTCTTCATTTAATTCTTCATTGTTATTTTCTT‑‑TCTTTTTGTGACTAGTGctttct                                                                      >  2:260723/1‑141 (MQ=255)

                                                                      tttattttaGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGctctct                                                                      >  1:1153478/1‑141 (MQ=255)

                                                                        tattttaGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCttt                                                                    <  1:1541417/141‑1 (MQ=255)

                                                                         attttaGATTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTACTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTGTT‑‑TCTTTTTGTGACTATTGCTTTCtttt                                                                   >  3:627068/1‑141 (MQ=255)

                                                                          ttttaGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCttttt                                                                  <  1:1513936/141‑1 (MQ=255)

                                                                            ttaGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCttttttt                                                                <  3:561936/141‑1 (MQ=255)

                                                                              aGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTTCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGc                                                              >  2:620792/1‑141 (MQ=255)

                                                                              aGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGc                                                              >  1:1513654/1‑141 (MQ=255)

                                                                              aGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGc                                                              >  1:1590592/1‑141 (MQ=255)

                                                                               gCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCa                                                             >  1:1545192/1‑141 (MQ=255)

                                                                                cTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCAc                                                            <  2:1201190/141‑1 (MQ=255)

                                                                                cTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCAc                                                            <  1:73529/141‑1 (MQ=255)

                                                                                cTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCAc                                                            <  2:238347/141‑1 (MQ=255)

                                                                                  tCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTa                                                          >  3:547926/1‑141 (MQ=255)

                                                                                  tCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTCAATTATTCATTGTCATTTTCTT‑‑TTTTTTTGTGACTAGTGCTCTCTTTTTTTGTACTa                                                          >  1:1397371/1‑141 (MQ=255)

                                                                                      aTTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCtt                                                      >  1:102496/1‑141 (MQ=255)

                                                                                         tGTCACTATCTATAGTTTCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTtgat                                                   <  2:924408/141‑1 (MQ=255)

                                                                                         tGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTtgat                                                   <  2:940059/141‑1 (MQ=255)

                                                                                          gTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTtgatg                                                  >  3:663068/1‑141 (MQ=255)

                                                                                           tCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTtgatga                                                 <  4:589086/141‑1 (MQ=255)

                                                                                           tCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTtgatga                                                 >  3:327000/1‑141 (MQ=255)

                                                                                           tCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTtgatga                                                 <  1:941801/141‑1 (MQ=255)

                                                                                             aCTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGatgatg                                               <  4:223978/141‑1 (MQ=255)

                                                                                               tAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAAttcttc                                                                                                             <  3:434459/79‑1 (MQ=255)

                                                                                                   aaTAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTTTCTCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCca                                           <  1:563012/141‑1 (MQ=255)

                                                                                                   aaTAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCcaca                                         <  1:855592/141‑1 (MQ=255)

                                                                                                    aTAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACaa                                        <  1:132359/141‑1 (MQ=255)

                                                                                                      aGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGc                                      <  1:628192/141‑1 (MQ=255)

                                                                                                       gTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTTTGACTAGTGCTCTCTTTTTTTGCACTATCTTTATGATGTCCACAAGCa                                     >  4:458965/1‑141 (MQ=255)

                                                                                                          tCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCt                                  <  1:467790/141‑1 (MQ=255)

                                                                                                                  aaTTTTTATATCCTTA‑CACTGTGTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAg                          <  1:924306/141‑1 (MQ=255)

                                                                                                                   aTTTTTATATCCATA‑‑TTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAgtg                        >  2:961672/1‑141 (MQ=255)

                                                                                                                    tttttATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAgtg                        <  4:501283/141‑1 (MQ=255)

                                                                                                                      tttATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAGTGTa                      >  3:588371/1‑141 (MQ=255)

                                                                                                                         atatCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAGTGTACtt                   <  2:106224/141‑1 (MQ=255)

                                                                                                                            tCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAGTGTACTTGCt                <  1:596437/141‑1 (MQ=255)

                                                                                                                            tCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAGTGTACTTGCt                <  1:549853/141‑1 (MQ=255)

                                                                                                                             ccATATTTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTTAGATATGTGTACTTGCt                <  2:295942/141‑1 (MQ=255)

                                                                                                                             ccATA‑TGTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTACTTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCCCCTAAGATAAGTGTACTTGCta               <  2:1581727/141‑1 (MQ=255)

                                                                                                                              cATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAGTGTACTTGCtaa              <  2:135050/141‑1 (MQ=255)

                                                                                                                                      tttttGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAGTGTACTTGCTAATAATATc       <  1:23520/141‑1 (MQ=255)

                                                                                                                                       ttttGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCCCAAGCACCTAAGATAAGTGTACTTGCTAATAATATcc      >  1:671505/1‑141 (MQ=255)

                                                                                                                                        tttGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAGTGTACTTGCTAATAATATccc     >  2:413222/1‑141 (MQ=255)

                                                                                                                                         ttGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAGTGTACTTGCTAATAATATcccc    <  4:365231/141‑1 (MQ=255)

                                                                                                                                           gCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCTCCTAAGATAAGTGTACTTGCTAATAATATCCCCAt  <  2:1404573/141‑1 (MQ=255)

                                                                                                                                           gCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAGTGTACTTGCTAATAATATCCCCAt  <  3:234428/141‑1 (MQ=255)

AGCAGTAATTTCCATTTTAAAAGCGATATATTTCTTTTTCTCATCTTTTTCATTGATGATAAACGGTTCTTTTATTTTAGCTTCAAATTTGTCACTAACAATAGTATCGCCTTTAATTTTTATATCCATA‑TTTTTTTTGCTTTTAAATTCTTTAAGTTCTTCATTTAATTCTTCATTGTCATTTTCTT‑‑TCTTTTTGTGACTAGTGCTCTCTTTTTTTGCACTATCTTGATGATGTCCACAAGCACCTAAGATAAGTGTACTTGCTAATAATATCCCCAT  >  USA300TCH1516_ALE/870754‑871032

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