BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,691,962	G→A	5.3%	F225F (TTC→TTT)	ydfJ ←	Membrane protein YdfJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,691,962	0	G	A	5.3%	278.1 / 3.5	95	F225F (TTC→TTT)	ydfJ	Membrane protein YdfJ
Reads supporting (aligned to +/- strand): ref base G (35/54); new base A (3/2); total (38/57)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

ATTGCAAGACCAATGCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTACAACAAATGCAACAATAATACCAACT  >  USA300TCH1516_ALE/2691822‑2692097

aTTTCACGACCAATCCCTTTAATGCAAACTGTTTCAGTTTTCATAACTTGTCGATATCTGAATGCAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCg                                                                                                                                         <  1:738920/141‑1 (MQ=255)

aTTGCAAGACCAATGCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCg                                                                                                                                         <  2:414207/141‑1 (MQ=255)

aTTGCAAGACCAATGCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCg                                                                                                                                         >  1:513396/1‑141 (MQ=255)

  tGCAAGACCAATGCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGaa                                                                                                                                       >  3:56364/1‑141 (MQ=255)

     aataCCAATTCCTCTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTTAATAAAATAAATAATTCATAATCGAAACCAACTACTAGTCCAATCATTACGGCTAATGTCCGTGTGACATTTGGTATATCGAATGc                                                                                                                                    <  1:526431/138‑1 (MQ=255)

     aaGACCAATGCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGc                                                                                                                                    >  1:658591/1‑141 (MQ=255)

      aGACCAATGCCTTTAATGTAATCTGTTTCAGTTTTCATAACTGGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCa                                                                                                                                   >  1:935697/1‑141 (MQ=255)

        aCCAATGCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATa                                                                                                                                 >  3:265547/1‑141 (MQ=255)

         ccAATGCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATaa                                                                                                                                >  1:841165/1‑141 (MQ=255)

          cgaTGCTTTTAATTTGATCTTTTTCGGTTTTCATAACTGGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGGCCAATGATTACGGCTAATGTCAGTGTTACATTTGGTATATCGAATGCATAAg                                                                                                                               <  2:893186/139‑1 (MQ=255)

          cAATTCCTTTAATGTAATCTTTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAg                                                                                                                               <  2:893112/141‑1 (MQ=255)

          cAATGCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACCGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAg                                                                                                                               <  2:72300/141‑1 (MQ=255)

           aaTGCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGt                                                                                                                              <  4:229716/141‑1 (MQ=255)

            aTGCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGtt                                                                                                                             >  2:976854/1‑141 (MQ=255)

              gCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTaa                                                                                                                           <  2:783796/141‑1 (MQ=255)

              gCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTaa                                                                                                                           <  2:446221/141‑1 (MQ=255)

                cTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCATCTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACa                                                                                                                         <  2:90350/141‑1 (MQ=255)

                cTTTAAAGTAATCTGTTTCAGTCTTCATAACTTCTCGTTCTCTGAATAAAATGAATAATGCATAATCGATACTATCTGGTAGTACAATCCTTACGCCTCCTGTCACTGTGACATTTGGTATATCCAATGCCTAAGTTACCa                                                                                                                         <  4:206718/141‑1 (MQ=255)

                  ttAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTaaaaaa                                                                                                                       >  1:23438/1‑141 (MQ=255)

                  tgaaTGTAATCTGGTTCAGGTTTCATAATTTGTCGATAGCTGAATAAAATAAATAGTGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACaaa                                                                                                                       <  2:582766/139‑1 (MQ=255)

                   tAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAAc                                                                                                                      >  3:303590/1‑141 (MQ=255)

                    aaTGTAATATGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACt                                                                                                                     <  4:333480/141‑1 (MQ=255)

                        tAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGata                                                                                                                 <  2:96900/141‑1 (MQ=255)

                         aaTCTGTTTCAGTTTTCATAACTTGTCGATACCTGAATAAAATAAACAAAGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAATGTGACATTTGGTATATCAAATGCATAAGTTAACAAACTGataa                                                                                                                >  2:996718/1‑141 (MQ=255)

                          aTCTGTTTCAGTTTTCATAACTTGTCGATACCTGAATAAAATAAACAAAGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAATGTGACATTTGGTATATCAAATGCATAAGTTAACAAACTGataat                                                                                                               <  2:59856/141‑1 (MQ=255)

                           tCTGTTTCAGTTTTTATGACCTGTAAATATCTGAATAATATAAATAATTCATAAATGATAACTACTGTCAATCCATACAATACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAACCTAACAAGCTGacgcta                                                                                                              <  2:844008/141‑6 (MQ=255)

                                ttCAGTTTTCATAGCTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATActacctac                                                                                                         <  2:885507/141‑1 (MQ=255)

                                  cAGTTTTNATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATAAATacac                                                                                                       <  2:816492/141‑1 (MQ=255)

                                  cAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGAGACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTacac                                                                                                       <  2:572661/141‑1 (MQ=255)

                                   aGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACAcc                                                                                                      <  2:118895/141‑1 (MQ=255)

                                   aGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGATAACAAACTGATAATACCTACAcc                                                                                                      >  1:808408/1‑141 (MQ=255)

                                     ttttCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCag                                                                                                    <  4:258836/141‑1 (MQ=255)

                                     ttttCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCag                                                                                                    <  4:439513/141‑1 (MQ=255)

                                      tttCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCaga                                                                                                   >  3:457180/1‑141 (MQ=255)

                                        tCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAgg                                                                                                 >  1:73698/1‑141 (MQ=255)

                                          aTAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCt                                                                                               >  1:771224/1‑141 (MQ=255)

                                          aTAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCt                                                                                               >  1:865588/1‑141 (MQ=255)

                                          aTAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCt                                                                                               <  4:94301/141‑1 (MQ=255)

                                             acgtGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGa                                                                                            <  2:177661/138‑1 (MQ=255)

                                              cTTGTCAATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCGACACCTGAGGCGAGAc                                                                                           <  2:846982/141‑1 (MQ=255)

                                               ttGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGAcc                                                                                          <  2:303869/141‑1 (MQ=255)

                                               ttGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCAAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGAcc                                                                                          >  1:548658/1‑141 (MQ=255)

                                                 gTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACcaa                                                                                        <  2:448453/141‑1 (MQ=255)

                                                   cGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATACTACCTACACCAGCGGCTAGACcaatc                                                                                      >  1:804787/1‑141 (MQ=255)

                                                   cGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACcaatc                                                                                      >  3:149572/1‑141 (MQ=255)

                                                          tGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAGTCAATGCAc                                                                               <  2:160559/141‑1 (MQ=255)

                                                          tGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCAc                                                                               <  4:361217/141‑1 (MQ=255)

                                                          tGAACAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGTCTAATGTCATTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCAc                                                                               <  3:38335/141‑1 (MQ=255)

                                                           gAATAAAATACAGAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCTGAGGCTAGACCAATCAATGCACt                                                                              <  2:2844/141‑1 (MQ=255)

                                                           gAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACt                                                                              <  2:903530/141‑1 (MQ=255)

                                                             ataaaataaaCAAAGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAATGTGACATTTGGTATATCAAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTta                                                                            >  1:310205/1‑141 (MQ=255)

                                                                  ataaataaTGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAAtt                                                                       <  4:298353/141‑1 (MQ=255)

                                                                     aataatGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGt                                                                    <  2:651410/141‑1 (MQ=255)

                                                                     aataatGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGt                                                                    <  4:218307/141‑1 (MQ=255)

                                                                     aataatGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTAGATCGAATGCATGGTTTAACGAAGTGTTAATTGCTATATCAGAGGCTAGACCATTCAATGCAGTTATATTTGGt                                                                    >  3:340633/1‑141 (MQ=255)

                                                                     aataatGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGGCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGt                                                                    >  3:57344/1‑141 (MQ=255)

                                                                     aataaaGCATAATCGATACCGACTGCTAGTCCAAGCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTTCACCAGCGGCTAGACCAATCAATGCTCTTATAATTGGt                                                                    <  2:770088/141‑1 (MQ=255)

                                                                       taatGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTaa                                                                  >  1:477497/1‑141 (MQ=255)

                                                                       taatGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTaa                                                                  >  2:167731/1‑141 (MQ=255)

                                                                          tGCATATTCGATACCAACTGCTAGTCCAATCATTACGGCCAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATcc                                                               <  1:799875/141‑1 (MQ=255)

                                                                            cATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTg                                                             <  1:164464/141‑1 (MQ=255)

                                                                            cATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTg                                                             <  2:515590/141‑1 (MQ=255)

                                                                             aTAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTCTATCGAATGCATAAGTTAAAAAACTGATAATACCTACAACAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTgc                                                            >  3:148536/1‑141 (MQ=255)

                                                                             aTAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTgc                                                            >  3:382711/1‑141 (MQ=255)

                                                                                 tCGATACNAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTtcagca                                                        <  2:813720/141‑1 (MQ=255)

                                                                                 tCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTccaaga                                                        >  1:547086/1‑138 (MQ=255)

                                                                                       ccAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACt                                                  >  1:542624/1‑141 (MQ=255)

                                                                                        cagcTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTg                                                 <  2:404586/138‑1 (MQ=255)

                                                                                         aaCTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGa                                                <  2:1047623/141‑1 (MQ=255)

                                                                                          aCTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGa                                                                                            >  1:895908/1‑96 (MQ=255)

                                                                                          aCTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGa                                                                                            <  2:895908/96‑1 (MQ=255)

                                                                                           cTGCTAGTCCGTTCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGTATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAACGCACTTATAATTGGTAATCCTGCAGCAATGACTGAAc                                              <  1:289083/141‑1 (MQ=255)

                                                                                           cTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTGACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAAc                                              <  2:843660/141‑1 (MQ=255)

                                                                                           cTGCCAGTCCAATCATTACGGATAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTAAACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAAc                                              <  1:661959/141‑1 (MQ=255)

                                                                                            tGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAAcc                                             >  1:520979/1‑141 (MQ=255)

                                                                                             gCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATGAGTGGTAATCCTGCAGCAATGACTGAACCg                                            <  2:389232/141‑1 (MQ=255)

                                                                                              cTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGa                                           <  2:561218/141‑1 (MQ=255)

                                                                                              cTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGa                                           >  3:421991/1‑141 (MQ=255)

                                                                                              cTAGTCCAATCATTACGGCTAATGTCAATGTGACATTTGGTATATCAAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGa                                           >  2:277636/1‑141 (MQ=255)

                                                                                               tAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGaa                                          <  2:710921/141‑1 (MQ=255)

                                                                                                   cgaaacATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGACCCGAAtgtg                                      <  4:368369/136‑1 (MQ=255)

                                                                                                    cAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGa                                     <  2:236930/141‑1 (MQ=255)

                                                                                                          ttACGGCTAATGTCAGTGTAACATTTAGTATATCGAAAGCATAAGTTAACAACATGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAAACCTGCAGAAATGACTGACCCGAATGTGATTAACa                               >  2:702117/1‑141 (MQ=255)

                                                                                                                cTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTacaa                         <  2:505427/141‑1 (MQ=255)

                                                                                                                cTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTacaa                         <  4:123736/141‑1 (MQ=255)

                                                                                                                   aTGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTAcaacaa                      <  2:615553/141‑1 (MQ=255)

                                                                                                                       cacactgACATTTGGTATATAGAATGCATAAGTTAACAAACTGATAATACCTACACCAGCGGCTAGACCAATCAATGCACTTATAATTGGTCAACCTGCAGCAATGACTGACACGAATCTGATTACCAGTGCAACAAATGc                  >  4:110930/6‑141 (MQ=255)

                                                                                                                       cAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTACAACAAATGc                  >  1:253595/1‑141 (MQ=255)

                                                                                                                       cAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTACAACAAATGc                  >  1:386927/1‑141 (MQ=255)

                                                                                                                       cAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCATGCAGCAATGACTGAACCGAATGTGATTAACAGTACAACAAATGc                  >  3:19728/1‑141 (MQ=255)

                                                                                                                              aCATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGAACTGACTGAAACGAATGTGGTGCACATTAAAACAAATGGGACaata           >  1:39436/1‑141 (MQ=255)

                                                                                                                                aTTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTACAACAAATGCAACaataat         >  1:152366/1‑141 (MQ=255)

                                                                                                                                   tGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTACAACAAATGCAACAATAATAcc      <  2:865588/141‑1 (MQ=255)

                                                                                                                                    ggTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTACAACAAATGCAACAATAATACCa     <  2:87691/141‑1 (MQ=255)

                                                                                                                                    ggTATATCGAATGCAGATGTTAACAAGCTAATAAGAACTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGAAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTACAACAAATGCAACAATAATACCa     <  2:486027/141‑1 (MQ=255)

                                                                                                                                     gTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTACAACAAATGCAACAATAATACCaa    >  1:752297/1‑141 (MQ=255)

                                                                                                                                      tgtatCAAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTACAACAAATGCAACAATAATACCAAc   <  3:262749/139‑1 (MQ=255)

                                                                                                                                       atatCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTACAACAAATGCAACAATAATACCAACt  <  4:270448/141‑1 (MQ=255)

ATTGCAAGACCAATGCCTTTAATGTAATCTGTTTCAGTTTTCATAACTTGTCGATATCTGAATAAAATAAATAATGCATAATCGATACCAACTGCTAGTCCAATCATTACGGCTAATGTCAGTGTGACATTTGGTATATCGAATGCATAAGTTAACAAACTGATAATACCTACACCAGAGGCTAGACCAATCAATGCACTTATAATTGGTAATCCTGCAGCAATGACTGAACCGAATGTGATTAACAGTACAACAAATGCAACAATAATACCAACT  >  USA300TCH1516_ALE/2691822‑2692097

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